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PREFACE 

This  Memorandum,  part  of  Rand's  research  on  the  use  of  airpower  in  support 
of  ground  operations,  is  one  of  a  series  of  studies  concerned  with  the  predic¬ 
tion  of  damage  to  a  single  or  multiple  targets.  Aspects  of  the  target-weapon 
relationship,  target-weapon  errors,  and  the  coverage  problem  were  considered  in 
RM-4566-PR.  FAST-VAL:  A  Theoretical  Approach  to  Some  General  Target  Coverage 
Problem ,  March  1966  (For  Official  Use  Only);  the  computer  program,  originally 
designed  to  compute  inputs  to  the  Rand  FAST-VAL  (Forward  Air  Strike  Evaluation) 
simulation  model,  was  described  in  a  companion  study,  RM-456  7-PR,  FAST-VAL:  Target 
Coverage  Model ,  March  1966  (For  Official  Use  Only). 

From  the  point  of  view  of  production  computations,  however,  the  target  cover¬ 
age  program  had  two  serious  limitations:  the  length  of  computer  time  required  to 
make  ths  computations,  and  dependence  of  precision  on  the.  size  of  the  integration 
cell,  which  in  turn  depends  on  the  machine  capacity  available.  To  alleviate  these 
problems,  a  model  was  designed  that  replaced  the  empirical  damage  function  used 
in  the  general  model  with  a  simpler  and  far  leas  time-consuming  analytic  expression. 
The  results  of  this  work  were  described  in  an  interim  reference  report,  RM-5152-PR, 

A  Simplified  Target  Coverage  Model ,  Novenier  1967  (For  Official  Use  Only).  As  these 
results  have  been  incorporated  in  the  present  study,  RM-5152-PR  has  been  withdrawn 
and  replaced  by  this  Memorandum. 

Thfe  work  of  simplification  continues  with  the  present  study,  which  describes 
a  model  that  is  sufficiently  broad  to  cover  almost  all  problems  that  arise  in  the 
field  of  nonnuclear  weapons  evaluation  and  a  machine  program  that  is  designed  to 
make  coverage  computations  for  some  of  the  problems  considered  in  RM-4567-PR  in  a 
fraction  of  the  time  previously  needed.  Without  materially  reducing  the  accuracy 
of  the  answers,  the  computer  time  required  for  some  problems  has  been  reduced  bv 
as  much  an  a  factor  of  100. 

After  continued  use  of  the  model  at  Rand  and  at  the  Air  Force  Armament  Labo¬ 
ratory,  further  refinements  to  the  program  to  make  it  more  efficient  both  as  to 
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time  and  accuracy  were  deemed  desirable.  An  improved  method  for  determining  the 
integration  step  size  was  implemented  with  a  considerable  increase  in  performance. 
The  method  of  fitting  an  empirical  bomblet  distribution  to  an  analytical  distri¬ 
bution  used  in  the  program  has  been  changed  so  that  any  smoothing  necessary  in 
irregular  data  is  done  by  the  machine  rather  than  by  the  user,  thus  avoiding  possi¬ 
ble  input  errors.  The  process  of  using  a  multiple  wing  station  configuration  was 
also  simplified.  A  direct  input  is  now  provided  for  the  numbers  of  wing  stations 
and  the  numbers  of  weapons  per  impulse  per  station. 
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SUMMARY 


Two  restrictions  permit  a  simplification  of  the  target  coverage  model  and  com¬ 
puter  program  described  in  FAST-VAL;  A  Theoretical  Approach  to  Some  General  Target 
Coverage  Problem  LU  and  m  FAST-VAL:  Target  Coverage  Model  [21.  The  problems 
are  restricted  to  the  case  of  a  gaussian  aiming  error  distribution  and  a  rectangu¬ 
lar  target  area  (uniform  distribution  of  the  target  elements  in  the  area).  As  a 
result  it  is  possible  to  reduce  the  coverage  computations  to  two  stages,  each  in¬ 
volving  a  double  integration,  in  contrast  to  the  three  stages  required  in  the 
original  model.  The  second  restriction  concerns  both  the  assumed  form  of  the 
damage  function  and  the  ballistic  error  distribution.  It  is  necessary  that  (1) 
the  damage  function  be  an  analytic  function,  rather  than  an  empirical  function; 

(2)  the  ballistic  error  distribution  be  one  of  three  types:  gaussian,  uniform, 
or  stick  type;  and  (3)  the  damage  function  b»  lntegrable  in  a  closed  form  with  re¬ 
spect  to  the  ballistic  error  distribution.  Under  these  restrictions,  the  coverage 
computations  are  reduced  to  a  single  stage,  involving  only  one  double  integration. 
For  some  cases,  it  is  possible  to  reduce  the  problem  to  a  summation  of  functions 
in  close-  form  with  no  integration  necessary. 

Two  types  of  damage  functions  are  considered,  corresponding  to  two  different 
types  of  weapon-target  effects:  a  fragment-sensitive  target,  or  one  in  which  the 
major  damage  mechanism  is  due  to  fragments  rather  than  to  a  direct  Impact  by  the 
weapon;  and  an  impact-sensitive  target,  or  one  for  which  there  is  a  definite  geo¬ 
metric  figure  that  must  be  Impacted  by  the  weapon.  For  the  fragment-sensitive 
target,  the  empirical  function  that  i >  usually  obtained  from  a  computer  program 
using  fragmentation  data  is  replaced  by  an  analytic  function,  the  "gaussian  damage 
function,"  fitted  through  the  choice  of  three  parameters.  For  an  impact-sensitive 
target,  it  is  assumed  that  the  target  element  is  e  rectangle,  and  that  there  '3  a 
fixed  probability  of  damage,  given  a  hit  on  the  element.  Under  these  assumptions, 
the  damage  function  is  exact. 

Stick  distribution  la  defined  in  App*  B. 
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The  two  types  of  damage  function,  the  aiming  and  ballistic  error  distributions, 
and  the  basic  coverage  relation  are  considered  in  turn,  and  the  expected  coverage 
K(X,Y)  is  expressed  as  one  den*’-  '-.tegration  in  terms  of  the  damage  function  and 
the  various  forms  of  the  aiming  and  ballistic  error  distribution.  For  various 
coverage  problems  that  occur  in  practice,  an  explicit  expression  for  the  coverage 
is  derived  in  terms  of  the  pertinent  parameters.  The  set  of  formulas  developed 
for  the  coverage  function-  provides  an  answer  to  the  weapon-target  effectiveness 
problem  that  corresponds  to  most  of  the  current  weapon  delivery  systems. 

The  FORTRAN  program  computes  numerical  answers  for  each  of  the  target  coverage 
problems.  The  output  of  the  program  is  the  particular  value  of  K  depending  on  the 
parameters  considered,  such  as  aiming  errors  ballistic  errors,  spacing,  and  weapon- 
target  effectiveness  indices. 
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1 .  INTRODUCTION 


A  number  of  general  target  coverage  problems  were  considered  in  FAST-VAL:  A 
Theoretioal  Approach  [1J.  The  companion  report,  FAST-VAL:  Target  Coverage  Model 
[2]  presents  a  machine  program  designed  to  make  coverage  computations  for  some  of 
the  problems  formulated  in  Ref.  1.  However,  from  the  point  of  view  of  production 
computations,  the  target  coverage  program  described  in  Ref.  1  and  Ref.  2  suffers 
from  two  serious  drawbacks:  lengthy  computation  time  and  dependence  of  precision 
on  integration  step  cell  size.  In  some  cases  a  single  problem  may  require  20  min¬ 
utes  of  computer  time,  and  the  cell  size  used  is  usually  limited  by  the  memory 
capacity  of  the  machine  as  well  as  the  time  element  in  computation.  To  alleviate 
these  problems,  a  model  was  designed  that  replaced  the  empirical  damage  function 
used  in  the  general  model  with  a  simpler  and  far  less  time-consuming  analytic  ex¬ 
pression.  The  results  of  this  work  are  described  in  Ref.  3,  A  Simplified  Target 
Coverage  Model,  and  have  been  incorporated  in  the  present  study.* 

Continuing  the  work  of  simplification,  the  present  study  describes  a  machine 
program  designed  to  make  coverage  computations  for  some  of  the  problems  considered 
in  Ref.  2  in  a  fraction  of  the  time  previously  needed.  The  program  in  Ref.  2  ac¬ 
complishes  the  coverage  computations  in  three  stages,  each  involving  a  double  in¬ 
tegration.  In  contrast,  this  study  restricts  the  problems  to  the  case  of  a  gaussian 
aiming  error  distribution  and  a  rectangular  target  area  (uniform  distribution  of 
the  target  elements  in  the  area).  As  a  result,  it  is  possible  to  reduce  the  cover¬ 
age  computations  to  two  stages,  each  involving  a  double  integration.  For  most 
problems,  these  restrictions  are  acceptable,  and  assumptions  of  a  gaussian  aiming 
error  distribution  and  of  a  rectangular  target  area  are  usually  valid.  The  second 
simplification  in  this  study  concerns  both  the  assumed  foim  of  the  damage  function 
and  the  ballistic  error  distribution.  It  is  necessary  that  (1)  the  damage  function 
be  an  analytic  function,  rather  than  an  empirical  function;  (2)  the  ballistic  error 

£  ^ 

Reference  3  will  be  withdrawn  and  replaced  by  this  Memorandum. 


distribution  be  one  ‘.nree  .types :  gaussian,  uniform,  or  stick  type;  and  (3)  the 
damage  function  be  integrable  in  closed  form  with  respect  to  the  ballistic  error 
distribution,  Under  these  restrictions,  the  coverage  computations  are  reduced  to 
s  single  stage,  involving  only  one  double  lntegretlon.  For  some  cases  it  is  pos¬ 
sible  to  reduce  the  problem  to  a  summation  of  functions  in  closed  form  with  no 
integration  necessary. 

Two  types  of  damage  functions  are  considered,  corresponding  to  two  different 
types  of  weapon-target  effects:  a  fragment-sensitive  target  and  an  impact-sensitive 
target.  For  the  fragment-sensitive  target,  the  damage  function  used  in  the  program 
in  Ref.  2  is  an  empirical  function  that  is  usually  obtained  from  a  computer  pro¬ 
gram  using  the  fragmentation  data  (size,  velocity,  and  distribution  of  fragments) 
from  arena  tests.  The  simplification  used  here  is  to  replace  this  matrix  function 
with  an  analytic  function,  the  "gausslan  damage  function"  [1],  fitted  through  the 
choice  of  three  parameters.  The  form  of  this  function  fits  most  of  the  empirical 
nonnuclear  weapon-effects  tables  very  well.  The  use  of  this  function  is  not  new. 

It  has  been  used  in  work  at  Sandia  Corporation,  where  it  is  called  the  "casualty 
function."  A  Ballistic  Research  Laboratories  report  [4]  considered  the  function 
for  the  special  case  of  an  elliptical  target.  Current  work  both  at  the  Strategic 
Air  Command  and  the  Naval  Weapons  Laboratory  has  involved  use  of  this  function. 

It  also  occurs  naturally  as  an  approximation  to  the  circular  coverage  function, 
where  it  has  been  referred  to  as  the  "Carleton"  approximation.  In  other  uses,  it 
has  been  called  the  "Carleton"  function. 

For  an  impact-sensitive  target,  i.c.,  one  for  which  there  is  a  definite  geo¬ 
metric  figure  that  must  be  impacted  by  the  weapon,  it  is  assumed  that  the  target 
element  is  a  rectangle,  and  that  there  is  a  fixed  probability  of  damage,  given  a 
hit  on  the  element.  Under  these  assumptions,  the  damage  function  is  exact,  l.e., 
the  same  as  would  be  used  in  the  program  of  Ref.  2. 

Section  2,  "Basic  Functions,"  considers  in  turn  the  two  types  of  damage  func¬ 
tions,  the  aiming  and  ballistic  error  distributions  to  be  used  and  the  basic  cover¬ 
age  relation,  which  expresses  the  expected  coverage  K(X,Y)  as  one  double  integration 
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in  terms  of  the  damage  function  and  the  various  forms  of  the  aiming  and  ballistic 
error  distributions. 

Sections  3  and  4  consider  various  coverage  problems  that  occur  in  practice. 

For  each  of  these,  there  is  derived  an  explicit  expression  for  the  coverage  K(X,Y) 
in  terms  of  the  pertinent  parameters.  Section  3  examines  the  ripple  delivery  of 
weapons  with  gaussian  distribution  of  both  the  ballistic  and  aiming  errors  (air 
delivered  bombs).  Included  is  one  case  of  an  elliptical  target  area  for  a  fragment' 

sensitive  target,  since  the  methods  used  here  are  applicable  to  this  case.  Various 

types  of  dispenser  delivery  of  weapons  are  considered  in  Sec.  4.  Section  4.1  dis¬ 
cusses  deliveries  for  fixed  dispensers  such  as  the  SUU-7  and  SUU-14  (CBU-1,  2,  3, 

7,  etc.).  For  these  cases,  the  ballistic  distribution  in  range  is  usually  avail¬ 
able  as  a  table  obtained  from  test  data  and  is  fed  into  the  model  of  Ref.  2  in 
this  form.  For  use  in  the  model  in  chJs  Memorandum,  this  empirical  distribution 
is  fitted  by  a  stick  distribution  through  the  use  of  two  parameters.  Section  4.2 
considers  the  case  of  a  ripple  of  dispensers,  released  from  the  carrier,  each  of 

which  spreads  its  subweapons  uniformly  over  a  rectangular  pattern  (Hayes  dispenser) 

A  ripple  of  dispensers,  released  from  the  carrier,  each  of  which  spraads  its  sub¬ 
weapons  uniformly  over  an  ellipse  (Rock-eye),  is  described  in  Sec.  4.3.  Section 
4.4  examines  the  same  case  for  dispensers  that  spread  their  subweapons  uniformly 
over  an  elliptical  annulus  (self-dispersing  subweapons,  such  as  the  CBU-24). 

Sections  3  and  4  thus  derive  a  set  of  formulas  for  the  coverage  function 
K(X,Y)  that  provide  the  answer  to  the  weapon-target  effectiveness  problem  corre¬ 
sponding  to  most  of  the  current  weapon  delivery  systems.  Section  5  describes  a 
machine  program  that  computes  numerical  answers  for  each  of  these  problems.  The 
output  of  the  machine  program  is  the  particular  value  of  K  depending  on  the  param¬ 
eters  considered  (aiming  errors,  ballistic  errors,  spacing,  weapon-target  effec¬ 
tiveness  indices,  etc.)  and  is  designated  by  X(j),  where  j  is  an  index  that  desig¬ 
nates  the  particular  type  of  problem:  X(l),  X(ll),  X ( 5 ) ,  and  X(15)  are  answers 
to  the  problems  considered  in  Sec.  3,  air-delivered  bombs;  X(2) ,  tf(12),  XO), 
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AT ( 13) ,  Y(4),  and  Y(14)  are  answers  to  the  problems  in  Sec.  4.1,  deliveries  from 
fixed  dispensers;  Y(103),  Y(104),  Y(105),  and  Y(106)  are  answers  to  problems  in 
Sec.  4.2,  dispensers  with  a  uniform  pattern  over  a  rectangle;  Y(107),  Y(108),  and 
^ (109)  are  answers  to  problems  in  Sec.  4.3,  dispensers  with  a  uniform  pattern  over 
an  ellipse;  ,Y(100),  A'(101),  Y(102),  Y(110),  Y(lll),  and  Y(112)  are  answers  to  prob¬ 
lems  considered  in  Sec.  4.4,  dispensers  with  a  uniform  pattern  over  an  annulus. 

Each  pertinent  formula  in  Secs.  3  and  4  is  designated  by  the  proper  X  designation 
to  correlate  it  with  the  corresponding  computer  output. 

Section  5,  a  description  of  the  machine  program,  was  planned  for  the  person 
interested  only  in  running  the  program  and  can  be  used  without  necessarily  refer¬ 
ring  to  the  previous  sections.  Thus,  some  description  of  the  various  problems 
considered  in  Secs.  3  and  4  is  repeated. 
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2,  BASIC  FUNCTIONS 

2.1.  DAMAGE  FUNCTIONS 

We  will  be  concerned  with  two  types  of  damage  functions,  depending  on  the 
weapon-target  damage  mechanism:  (1)  a  fragment-sensitive  target  and  (2)  an  impact- 
sensitive  target  with  a  defined  area  In  which  a  hit  is  necessary  for  damage,  i.e. , 
the  probability  of  damage  if  hit,  p^,  may  be  less  than  1. 

2.1.1.  Fragment-sensitive  Gaussian  Damage  Functicn 

As  the  name  implies,  a  fragment-sensitive  weapon-target  case  is  one  in  which 
the  major  damage  mechanism  is  due  to  fragments  rather  than  to  a  direct  impact  by 
the  weapon.  In  general,  the  damage  function  will  be  computed  in  matrix  form.  We 
will  assume  the  form  of  the  damage  function  to  be 

(2.1)  0(*. *)-£><,  e*p(-Z)(I"'-2-~+-^- — *1 V 

\  U2(l>  *2( 2>J/ 

where  D(x,y)  la  the  probability  that  a  target  element  at  ( x,y )  is  "killed"  (damaged 
at  least  to  a  specified  degree),  and  J?(l),  R( 2),  and  DQ  are  parameters  to  be  ob¬ 
tained  from  empirical  data. 

In  Ref.  1,  Sec.  1.6.3,  certain  parametera  are  discussed  that  may  be  used  to 
characterize  a  damage  function.  The  weapon  radius  R  is  defined  as 


Using  (2.1),  ve  find  that 

(2.2)  R2-R(l)R(2). 

Along  any  ray  0,  the  mean  square  radial  damage  distance  r2(6)  is,  if  (2.1)  is  used, 

_  ,• 

(2.3)  r2(6)«2j  rD( r  cos9,r  ain9)cfr 


1 

x2(e) 


where 


*2(0)«£°dV-^. 

ff2(l)  R7  (2) 


Letting  0=0°  and  0=90°,  we  find  that 


fl2(l)=[  2xD(x ,0)dx*r2 (0°) , 
J0 

(2.4) 


fl2(2)=J  2yD(Oty)dy~r2(W°). 


Thus,  f?(l)  and  R{ 2)  in  (2.1)  are  simply  the  mean  square  radial  damage  distances 
along  the  x  and  y  axes,  respectively. 

Considering  the  mean  radial  damage  distance  r(6)  and  using  (2.1),  we  obtain 


(2.5) 


i*(6)«  [  Dir  cos8,i*  sin0)dr 

Jo 

m2XW' 


The  relative  damage  variance  E2(9)  is 


r2(e)  4 


(2.6) 


Do 


,4[i-s2(e)] 


Thus,  we  can  view  the  three  parameters  in  (2.1)  as  the  triplet  C)f(l) ,R\2) , E] , 
where  tf(l)  and  /?( 2)  are  the  mean  square  radial  distances  along  the  x  and  y  axes, 
and  E2  is  the  relative  damage  variance.  Instead  of  1?(1)  and  R( 2),  we  can  use  R 
and  p,  considering  the  triplet  (ff,p,E),  where 


2.2.  DELIVERY  ERRORS 

In  general,  we  shall  consider  the  case  of  multiple  delivery  of  weapons.  Each 
weapon  is  assumed  to  be  subject  to  a  ballistic  error  in  range  and  deflection,  ex¬ 
pressed  by  a  distribution  function  T(x,y) ,  independent  of  the  other  weapons.  The 
whole  delivery  is  considered  to  be  subject  to  aiming  errors  in  range  and  deflection, 
<  expressed  by  another  distribution  function. 
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2.2.1.  Aiming  Error  Distributions 

We  shall  confine  our  attention  to  the  case  of  the  aiming  error  distribution 
having  the  form  of  a  noncorrelated  bivariate  gaussian,  i.e., 


(2.10) 


Prob(*S)f,i/sy)  =  r  f  exp^-4|— - - -|\ 

J.J.oo  \2U2(1)  t2 (2)J/ 


dxdit 


2nt(l)t(2> 


where  t(l)  and  fc( 2)  are  the  standard  deviations  in  the  x  and  y  directions,  respec¬ 
tively.  For  convenience,  we  define  the  linear  gaussian  distribution  G(X)  and  its 
density  function  g(X)=G' (X)  as  follows: 

(2.11) 

(x)-.e.^.(-£2/2) 

9  Jli 

Thus,  Eq.  (2.10)  can  be  written  as 

(2.12)  Prob  (*S*.»sr>"<;(tfr))<,(tf2)) ' 

If  we  wish  to  use  REP  (range  probable  error)  and  DEP  (deflection  probable  error), 
the  relations  between  them  and  the  standard  deviations  are 


(2.13) 


REP-.6744t(l) , 
DEP- ,6744£ (2) . 


Further,  if  CEP  (circular  probable  error)  is  given  by 

(2,14)  CEP-1. 1774t-1.7459REP. 

2.2,2.  Ballistic  Error  Distributions 

We  shall  assume  that  the  ballistic  errors  are  independent  in  the  x  and  y 
directions.  Therefore,  the  ballistic  error  distribution  T(X,Y)  is  of  the  form 
r(.Y,lr)-7’j  (X)72(7) .  We  shall  consider  three  types  of  linear  ballistic  error  dis¬ 
tributions:  gaussian,  uniform,  and  stick. 
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a.  Gaussian  Ballistic  Error.  As  for  the  aiming  error,  we  assume  for  the 
first  ballistic  error  a  linear  gaussian  error  distribution,  i.e., 


(2.15) 


7 1 (X) -6 1 


dx 

/jiiied)  ’ 


If  the  errors  are  gaussian  in  both  directions,  we  have 


(2.16) 


7«-»-<5TrT>(sy- 


b.  Uniform  Ballistic  Error.  For  the  second  ballistic  error,  we  assume  a 
uniform  distribution,  i.e., 

(2.17)  TiUO-J  a(x,L)^, 


where  a(x,L)  is  defined  in  (2.9).  A  similar  expression  holds  in  the  y  direction. 

c.  Stick  Ballistic  Error.  A  third  type  of  ballistic  error  is  the  so-called 
stick  type,  which  is  a  combination  of  a  uniform  and  a  gaussian  distribution,  dis¬ 
cussed  in  App.  B.  As  given  in  (BA),  we  have 


(2.18) 


wherj 

J 

R{XtL)=\  Mx,L)dx, 

J  —00 

(2.19) 

l  rL 

7i(x,L)“2^J  g(x*y)dy. 


If  Lemma  3  in  App.  D  is  used,  the  function  H(.X,L)  can  be  evaluated,  for  L* 0,  as 

(2.20)  H{XtL)-(.X+L)G{X+L)~(X-L)C(.X-L)+g(X+l)-g<.X-L) , 

where  G(X)  and  g(x)  are  defined  in  (2.11). 
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We  note  the  following  limits  for  h  and  H: 


(2.21) 


1.h(x/e,L/e)  ^(x/b) 

e*  rt  * 


L-+  0 


,,_h{xla,Ll8)  a (x,L) 
~ . . 2L  * 

8"*0 


lirnWl 

L*  0 


mo. 


limtfl 

8-»0 


2.3.  THE  BASIC  COVERAGE  RELATIONS 

Define  X(X,Y)  as  the  fractional  coverage  to  a  target  area  with  center  at  (A',y) 
inflicted  by  a  pattern  of  weapons  subject  to  a  guussian  aiming  error  with  standard 
deviations  t(l)  and  t( 2).  The  pattern  damage  function  D^(.u,v)  is  the  kill  prob¬ 
ability  to  a  target  element  at  ( u,v )  if  the  center  of  the  pattern  is  at  (0,0). 

The  function  D^(.u,v)  depends  on  the  specific  problem  and  will  be  determined  for 
each  case. 

Define  K^(.u,v)  as  the  point  target  kill  probability  to  a  target  element  at 
( u,v )  if  the  aiming  point  is  at  (0,0).  We  obtain 


(2,22) 


Kp(u,v) 


(u-x,V-y)g  (^7^)  ^(tfer) 


dxdy 

*(l)t(2)‘ 


The  fractional  coverage  K(X,Y)  for  a  target  area  A  is  thus 


(2.23) 


W,y) 


(u,u) 


dudv 

A 


Dp(u-x,v-y)g(^y-%~) 


dxdydudv 


Consider  a  rectangular  area  of  dimensions  I  2/1(1)  ,2/1(2)  j. 


Using  (2.23),  we 


obtain 
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( l )  (2)  „  .  /  «  \  /  y  \  dxdududv 

(2.24)  K(*,Y)-[  J  I  ^p^-**y“y)ff\.T(T)'^U(2)/4t(l)t(2M(lM(2)' 


dxdududv 


'x-Aixyy-A'i) 


Tnrough  appropriate  changes  in  variable  and  reversing  the  order  of  integration, 


we  obtain 


(2.25) 


•  -  r  ,  Jx+AO)]/t(i) 

X(X,ry  1  1  D  (3*X,y+Y)  tott  ,  g(u)du 


x  Sy+A (2)  ]/t (2) 


_ir^+' 

24(2)Jr 

L  [y-< 


•X(2)]/t(2) 


g(.v)dv  \dxdy. 


Using  the  definition  of  h(x,L)  from  (2.19)  in  (2.25),  we  get 


(2.26)  *a’Mj 


If  Op(a :,t/)  can  be  expressed  without  integrals,  we  have  thus  reduced  the  problem 
of  obtaining  K  to  a  double  integration.  An  alternate  form  for  K  used  in  the  nu¬ 
merical  procedure  is 


(2.27)  *(*,*) -J  J  D pC«t(l)+Y,^(2)+r]7i(*,|^^(y,^)d*di/. 


For  a  gaussian  target,  i.e.,  one  whose  location  ( u,v )  is  subject  to  a  gaus- 
sian  distribution  with  variances  t2( 3)  and  t2(41  and  center  at  (Y,Y),  the  prob¬ 
ability  of  kill  K(X,Y)  is 


(2.28) 


K(X,Y) 


dudv 

ti 3)t(4)’ 


where  Ap(w,v)  is  given  in  (2.22).  Substituting  Eq.  (2.22)  in  (2.28)  and  using 
Lemma  1  in  App.  D,  we  obtain 
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(2.29)  •  uxt  r)-|  J_  j>p 


where 


(2.30) 


t2(5)-t2(l)+t2(3), 

t2(6)-t2(2)+i2(4). 


Now  from  (2.21),  we  get 


Thus  we  can  put  X  in  (2.29)  In  the  form 

(2.32)  x(*,*)»J  J  J>p (^ 'y+y)h{m)>Q)h($6)' °)t(WiW 

When  we  compare  this  equation  with  Eq.  (2.26),  we  see  that  the  problem  is  equiva¬ 
lent  to  that  of  a  point  target  with  the  aiming  error  variances  equal  to  the  sum 
of  the  actual  aiming  error  variances  and  the  target  location  variances. 
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3.  RIPPLE  DELIVERY  WITH  GAUSSIAN  DELIVERY  ERRORS 


3.1.  RECTANGULAR  TARGET  AREA:  FRAGMENT-SENSITIVE  T/.RC-ET 

Let  us  consider  first  a  single  weapon  against  a  target  element  at  (u,v)  sub¬ 
ject  to  ballistic  errors  according  to  a  gaussian  distribution  with  standard  devi¬ 
ations  a(l)  and  8(2).  If  the  ballistic  center  of  impact  is  at  (0,0),  the  kill 
probability  D^(u,v)  is  given  by 

( 3 . 1 )  Db(u,v)~^J-D(i*-xti>-y)  ?(ife)8(nf(2)’ 


where  D(x,y)  is  the  damage  function  in  (2.1).  Substituting  for  D  from  (2.1)  and 
integrating,  using  Lemma  1  of  App.  D,  we  obtain 


(3.2) 


Db(u,v) 


.JSDM21 

2q(l)q(2) 


where 


q2(l)"^2p~+«2(1)  > 

q2(2)-^li)+82(2). 


Consider  the  delivery  of  N(l)  weapons  against  a  target  element  at  (u,v) . 
Each  weapon  is  subject  to  ballistic  errors  as  above.  The  #(1)  weapons  will  be 
delivered  in  a  ripple  mode.  Examine  first  a  zero  aiming  error.  Then  the  tth 
weapon  will  have  a  center  of  impact  (Cl)  at  (£. ,n«).  The  array  of  Cl's  form  the 
ripple  pattern.  We  will  use  the  center  of  gravity  of  the  pattern  as  the  pattern 
reference  point;  i.e.,  if  the  reference  point  is  at  (x,y),  then 


(3.3) 


We  will  measure  a.  and  b.  from  the  reference  point,  i.e., 

w*. 
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Thus,  we  obtain 


(3.4) 


/Ml) 

y  b.*o. 


Consider  the  pattern  reference  point  at  ( x.y )  and  a  single  weapon  subject 
only  to  ballistic  error.  The  center  of  impact  of  the  ith  weapon  is  thus  at 
(aH-a tffb .) .  Let  p  be  the  weapon  reliability  factor,  i.e.,  the  probability  of 
the  weapon  working  properly.  If  the  center  of  the  pattern  is  at  (0,0),  the  kill 
probability  D^{u,v)  due  to  the  pattern  is  given  by 


N(\) 

(3.5)  £p(u, v)«l-  JJ'  [1  -pZ^u-a^.y-b.)  ], 


where  D ^  is  given  by  (4.2).  Explicitly,  we  have 


(3.6) 


tf(l) 

V“.y)- 1-  fl1 

i- 1 


1  Pi?0)g(2) 
l'2q(l)q{2) 


Using  this  expression  for  D ^  in  (2.27),  we  obtain  the  coverage  X(X,Y)'. 


(3.7) 


v  t-i 


\2?L  <?(1)  . 


Equation  (3.7)  is  the  basis  for  obtaining  ( 1 1 )  in  the  machine  computation  program 
discussed  in  Sec,  5. 

Expanding  the  product  in  (3.6),  we  obtain  as  a  sum,  which,  when  substituted 
in  (2,24),  gives  an  alternate  expression  for  K(XtY)  that  can  be  integrated  term  by 
term.  Thus,  after  expansion,  we  have 
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KU.Y)-”)"  (~l)k~  f  °~fef  1  f1'-*2  N p> 

*  U<i)^<2)J  frl  Lllfl  -  J.  , 


2-1  !+l  VW 


.X+MD  » 


”*-/!(  1) 


f  exP(^  f — jL — Y/_g__\_..4Kft>  pM(I)r 


2,2(1) 


/  k^.ry)A 

x  expl-  |  - — -£ —  Uj'-i — \ — mmki. 

\A  2,2(2)  fwvhtwAur 

Using  Lemma  5  in  Apa.  D  to  evaluate  the  inner  integrals  in  (3.8),  we  obtain 

N(  1 ) 

(3.9)  K{X,Y).  I  (- 1 ) 1  fg/?Q)/?(2)] *fa(l)q(2)  ~| 

g i  L2?(1),(2)J  LW)«(2)J 


(WU>\ 

V  k  4*M(i) 

/  tu-a 

•  I J 

expl - 

Jx-A(l) 

\ 

/ 

[y-b(t1,t; 

j  e*p- 

Jr->4(2)  \ 

2  Q‘ 

'fc[oy(i1,i2,...,i  . t.)]) 


’•y  “ 
2q2(I) 


'  W(ii,i2,...t^)-g2tflti2 . ju) 

2,2(2)  > 


:ii.i2....,ifc)4l  a|  .  5(ii,£2 . £*)4  \  a.  , 

,•»  3  ,7-1  j 


fc2(il^2 . V"i,£‘ bi.*  i  b.  , 

j  *  vi,  V 


where 


(3.10) 


<?2(2)=4^+fl2<2). 


Q2(  2)  =  >-  (2)-i'I-rP  =t  ‘  (2)4/--(-2-)-^|-^--6'.  •^■, 


wu)\  /y(D ! 

y,  k  )  *k\  ( /V(l)^£]T  * 


N(l)-k+l  N(l)-k+2 


•  1  l 

tj*!  l2=ti+l 


N(.  1) 

VV 1+1 


(T})  //v(  i  >\ 

The  sum  £  is  the  sum  over  the  \  k  )  combinations  of  the  N(l)  integers 
1,2,3,. ..  ,N(l)  taken  k  at  a  time. 

If  we  have  a  point  target,  i.e.,  4(1)*j4(2)*0,  Eq.  (3.9)  for  /(becomes 


(3.11) 


Nil) 

w,y>=  l 

k- 1 


(-1) 


/c-l 


fp/?(l)/?(2)1  V  g(l)q(2)~| 

|_2<7 ( 1  )<y ( 2)J  Lfe«(l)«(2)J 


In  the  special  case  of  a.Mb.= 0,  i.e.,  salvo  firs  against  a  point  target,  the 
(N(  1)\  t  1 

V  k  )  (H  lj)  _  _  _2 

sum  1  in  (3.11)  is  simply  \  k  /  since  a2ma2mb2"b  *0.  Thus,  we  obtain 


JV(D  » 
(3.12)  /((A-,7)=  7  (-1) 

Jc=i 


y?  > 

2Q?(2); 
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For  an  area  target,  if  the  definition  of  h{xtL)  given  in  (2.19)  is  used,  Eq. 
(3.9)  for  K  becomes 


(3.13) 


K(X  ><7(2)1 

(  *  kil  (  °  L*^CI><f(2)J  r[<?(l)C(2)J 

(»<») 

*  y  h(£*~  diil)h(M-  A<2A 

L  n\Q(l)  •Q(iyn\Qt2VQi2)) 


(k[az(t  i,i2.  •  •  •  ,ik)-a2U- \,iz . ife)  A 

2??(1)  / 


-  ,  - 

(k(.fcZ(t- 1  ,t-2 » •  •  •  (^lf  l2»  •  ♦  •  »^)  A 

2e?2(2)  / 


where  the  symbols  are  defined  in  (3.10). 

In  the  special  case  of  salvo  fire,  i.e..  0.“fc .*0,  the  expression  for  K(.X,Y) 

f<  I* 

in  (3.13)  simplifies  to 

,,,  „,v  .,fc-l/ff(lArpm)/?(2)lWa(l)<7(2)1,V  X  A(l)\,./  y  4(2)\ 

(3.H)  X(y.y).^  (-1)  ^  ;L^(T)?(2)J 


Equations  (3.13)  and  (3.14)  form  the  basis  for  obtaining  X(l)  in  the  machine  pro¬ 
gram.  Thus,  X(l)  and  X(ll)  are  equivalent  expressions  in  which  different  compu¬ 
tation  techniques  are  used. 


3.2.  RECTANGULAR  TARGET  AREA:  IMPACT-SENSITIVE  TARGET 


Consider  a  single  weapon  subject  only  to  ballistic  errors.  Substituting 
the  impact  damage  function  D  from  (2.8)  in  the  expression  for  D ^  in  (3.1),  we  ob¬ 
tain  for  D^(u,v)  the  equation 


(3.15) 


.u+3(l)  .V+B(2)  ,  _  .  .  .  ,  , 

Db(u’v)mpdH[u_B(l)  lv_B(2Xt mX^Tcrrfh)' 


If  we  define  f(xty)  as 


(3.16)  f(.xty)"G(x*-y)-G(x-y) , 

from  (3.15)  becomes 

(3<l7)  V  “  • w)  . 

Proceeding  exactly  as  in  the  preceding  section  on  a  fragment-sensitive  target 
for  a  ripple  of  N( 1)  weapons,  we  obtain  an  expression  for  D^(u,v)  corresponding 
to  (3.5): 


HD 

(3.18)  £>  (tt,t>)- 1-  Jj  U-p‘ Db(u-a^tV-bj)) , 

t-1 


where  D b  is  given  in  (3.17)  and 


(3.19) 


P'-PPj*. 


i.e.,  p'  is  the  product  of  the  reliability  factor  for  the  weapon  and  the  probability 
of  damage  if  hit.  Using  (3.17)  in  (3.18)  and  substituting  in  (2.27),  we  obtain  for 
K(X,Y) 


(3.20) 


<•  SS11 

’8(1) 


'yt(2)+Y-b 


8(2) 


t  ssm 


Eduction  (3.20)  is  the  hssis 


for  obtaining  X(15)  in  the  machine  program. 


Expanding  the  product  in  (3.2Q1,  we  obtain  an  alternate  expression  for  K,  l.e.. 
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*CP  k.  x  .. 

(3.21)  K(XtT)-  l  (-1)*  l 
h*\ 


^(D\ 

t(Dj 


dx 


<yt{2)+Y-bi 

L  8(2) 


im. 

’em, 


We  note  that:  we  have  a  sum  of  products  involving  only  a  single  integration  rather 
than  a  double  integration  as  in  (3.20).  For  the  salvo  case,  we  obtain  from  (3.21) 


(3.22)  /WO-T5  <-!>*"  W}) 
k- l  v  *  ' 


,  f  Mm m  mj\hL  f  Mumhi  mi)hL  AmU 

J_/  \  •<«  ’«0>/  V  ’*(!>/  J _/  V  8^2)  ’e(2))  \y,tmJay’ 


Equations  (3.21)  and  (3.22)  form  the  basis  for  obtaining  .Y(5)  in  the  machine  pro¬ 
gram.  Again  *(5)  and  Jf(  1 5)  are  equivalent  expressions  for  K,  but  differ  in  the 
method  of  computation. 


3.3.  ELLIPTICAL  TARGET  AREA:  FRAGMENT-SENSITIVE  TARGET 

Consider  an  elliptical  target  area  with  axes  4(3)  and  4(4)  in  the  x  and  y 
directions.  The  fraction  coverage  K(.X,X)  for  this  target  area  is  given  by  (2.23). 
Specifically,  we  have 


(3.23)  wr,«. 


./  *  \ 

J  y  \ 

\£(1)/ 

Atm) 

t(2)/rt(l )fc (2)4 (3)4(4) ’ 


where  4  is  defined  by 


For  a  fragment-sensitive  target,  the  pattern  damage  function  D^{u,v)  is  given  by 
Eq.  (3.6),  the  same  as  was  used  in  Sec.  3.!.  Expanding  the  product  in  (3.6), 


-20- 


substituting  in  (3.3),  and  integrating  as  in  F.q.  (3.9),  we  obtain  an  expression 
for  K  that  is  similar  to  (3.9): 


(r) 


,.,r  n  NV1) .  ..kritmimnVr  a(l)a(2)lV  £  '  (  k(a2-d2)\  ,  (  Hb2-b?)\ 

J,  l‘1>  LWnj  btrnfeyj  z  T^f]  w»T/' 


(3.24) 


|,J  “p  * j)  > 


dudv 


f  Wi4  (3)4  (4)  ’ 

where  i4  is  given  in  (3.23)  and  the  other  symbols  are  as  defined  in  (3.10). 


Consider  the  integral  expression  I in  (3.24);  i.e., 


(3.23) 


where 


If 


dudv 

x/l(3)i4(4) 


2g(l)g(2) 

4(3)4(4) 


<■*•*>)#". 


C  .  Uu-(^-a)J/g(l)}2|U»-Q,-fe)]/0(2)}2.:lj 

r  U(3)/g(l)j2  r/l(4)/g(2)12 

Define  the  offset  ellipse  function  P{A,B’tx,y)  as  the  integral  of  a  gaussian  dis¬ 
tribution  of  unit  standard  deviation  over  an  offset  ellipse,  center  at  ( x,y ),  with 
axes  A  and  B  in  the  x  and  y  directions,  respectively;  i.e., 

(3.26)  P(A  ,B;x,i/)  =  JJ  exp(-|(i;2+n2))^~, 


A 2  S2 


where 


where  -Tq (*)  ia  the  zero-order  Beaael  function  of  the  first  kind  with  imaginary 
argument. 

Conaider  the  apecial  case 

«<l)-<?(2)  an<J  4(3)-4(4), 
i.e.,  the  target  is  a  circle  of  radiua  4(3). 
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(3.29) 


H 1) 

M,Y)~  l  (-1) 
A>  l 


k- 1 


exp  -2 


A:  a2-a2.b2-b 

■  ' 


p4o)  ^X-*n2+(r-2>)2\ 

50)’  «(l) 


where  P(R,r )  is  defined  in  (3.28).  For  the  salvo  case,  this  expreseion  reduces  to 


4,  DISPENSER- TYPE  DELIVERIES  AGAINST  RECTANGULAR  TARGET  AREAS 

Let  us  consider  a  dispenser-type  delivery.  This  type  includes  both  the  mass 
delivery  of  a  number  of  weapons  and  the  delivery  of  a  sequence  of  weapons  at  timed 
intervals.  The  latter  type  is  similar  to  the  ripple  delivery  of  the  previous  sec¬ 
tion  except  that  the  ballistic  dispersion  of  the  individual  weapon  is  not  neces¬ 
sarily  gaus8ian,  and  the  number  delivered  may  be  much  larger  than  in  the  case  of 
Sec.  3. 

The  dispensers  themselves  may  stay  with  the  delivery  vehicle  (no  ballistic 
errors  for  the  dispenser)  or  may  be  released  in  the  ripple  mode  and  thus  be  sub¬ 
ject  to  dispenser  ballistic  errors.  In  either  case,  the  subweapons  are  distributed 
according  to  some  type  of  pattern  distribution.  In  addition  to  the  ballistic  dis¬ 
tributions  of  Sec.  3,  we  will  also  consider  cases  in  which  the  subweapons  are  uni¬ 
formly  distributed  over  a  pattern  (rectangle,  ellipse,  annulus).  We  will  view 
these  pattern  distributions  as  ballistic  error  distributions  for  the  individual 
subweapons.  Thus,  we  have  cases  where  there  are  two  ballistic  error  distributions, 
one  for  the  dispenser  and  one  for  the  individual  subweapons. 

We  will  consider,  in  turn,  the  following  types  of  subweapon  distribution 
functions: 

1.  Stick  distribution  in  range;  gausslan  distribution  in  deflection. 

2.  Uniform  distribution  in  range  and  deflection  (over  a  rectangle). 

3.  Uniform  distribution  over  an  ellipse. 

4.  Uniform  distribution  over  an  elliptic  annulus. 

In  this  section  we  restrict  ourselves  to  the  consideration  of  either  point 
targets  or  rectangular  target  areas  of  dimensions  [24(0,24(2) J. 

4.1 .  RIPPLE-TYPE  DISPENSER  DELIVERIES 

We  consider  first  a  sequence  of  #( 2)  weapons  delivered  at  timed  intervals, 

i.e.,  with  centers  of  impact  forming  a  pattern  (a.  ,fc.),t* 1 ,2,3, ... ,17(2) .  This 

* 

,  type  of  delivery  is  similar  to  the  type  of  delivery  of  Sec.  3  except  that  several 


patterns  may  be  delivered  at  once,  either  with  the  same  aiming  points  or  with 
offset  aiming  points,  which  involve  the  delivery  of  several  tubes  in  a  container 
or  from  several  containers.  Further,  the  ballistic  data  for  the  subweapon  are 
usually  not  obtained  for  each  weapon  individually.  The  ballistic  errors  in  de¬ 
flection  are  assumed  to  be  gaussian,  but  the  range  errors  are  governed  by  an  em¬ 
pirical  distribution. 

We  assume  each  of  the  tf(2)  weapons  to  be  aimed  individually  at  the  center 
point  of  the  pattern.  We  fit  the  empirical  data  for  range  dispersion  by  assuming 
that  its  distribution  K\(.X)  is  the  stick  distribution: 


(4.1) 


K^X) 


r  h(-z~  _ 

%(0)  *q(0)J  J_JV°>  ’q{0))q(0)  * 


dx 


where  L(l)  is  the  "length"  of  the  stick,  and  q( 0)  is  a  parameter  of  the  distri¬ 
bution.  Appendix  B  contains  a  discussion  of  this  pattern. 

The  ballistic  distribution  in  deflection  is  gaussian  with  standard  deviation 
s(2).  Consider  a  single  weapon  against  a  target  element  at  (u,v) .  The  kill  prob¬ 
ability  D^{u,V)  is  given  by 


(4.2) 


\  dxdu 


fi(2)/e(2)q(0)  * 


where  D(u-x,v-y )  is  the  damage  function.  If  N(2)  weapons  are  released  from  the 
dispenser  with  pattern  center  at  (0,0),  the  kill  probability  D^(u,v)  to  a  target 
element  at  (u,V)  due  to  a  single  dispenser  is 


(4.3) 


0p(u,tf)-l-[l-p0fc(u,y)/(2), 


where  D^{u,V)  is  as  given  in  (4,2),  and  p  is  the  subweapon  reliability  factor.  We 
define  v  as  the  probability  of  a  dispenser  operating  properly,  i.e.,  the  dispenser 
reliability  factor.  If  there  are  N(l)  dispensers  with  pattern  centers  at  ( a-,b .), 
a  pattern  of  dispensers  is  formed.  Thus,  the  kill  probability  D^(.u,v)  due  to  the 
whole  pattern  is 
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#(1) 

(4.4)  Cp(M,t>)-l-  [1  -rDjjiu-a^tV-b^)] 

i"  l 

#U> 

-1-  rf  Cl-i»{l-[l-pDi(u-ai,y-bi)ru,>]. 

t-l 


4.1.1.  Fragment-sensitive  Targets 

For  a  fragment-sensitive  target,  substituting  the  damage  function  (2.1)  in 
(4.2),  we  obtain  for  D^(u,v)  the  expression 


(4.5) 


The  integrals  are  evaluated  using  App.  D;  Lemma  6  is  used  for  the  first  integral 
and  Lemma  1  for  the  second  Integral.  We  obtain 


(4.6) 


where 


<?2(2)-^P+82(2), 

(4.7) 

?2(3)-^~+<72(0). 

Using  (4.6)  in  (4.4)  to  get  D  (u,v)  and  substituting  in  (3.27)  to  obtain  the 

P 

coverage  K(X,i)t  we  have 
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(A.8) 


If  r«l,  we  obtain 


(A. 9) 


K(.x,y )- 


1- 


>^ff(i)./gt(1HX~at  loA 

■?(3)  \  <?(3>  V3)/ 


X 


«(2) 


The  expression  for  K(X,Y)  in  (A. 9)  forms  the  basis  for  X(12)  in  the  machine 
program. 

For  a  single  dispenser,  N(  1)-1,  we  expand  the  expression  in  (A. 9)  to  obtain 


The  second  integral  may  be  evaluated  using  Lemma  6  of  App.  D,  so  chat 


<,  u,  w.v-Y  <-»*-(?’)  mfflM&MfJ MH9*- 


where  q(2)  and  <7(3)  are  given  in  (A. 7)  and 

,2 


Q2W-3-j~+t2{2). 
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Equation  (A. 11)  forma  the  basis  for  £(2)  in  the  machine  program,  and  X(2)  and  £(12) 
are  thus  also  equivalent  except  for  the  difference  in  computational  methods. 


4.1.2.  Impact-sensitive  Targets 

For  a  hard  target-type  subweapon,  substituting  the  damage  function  (3.1)  in 
(A. 2),  we  obtain  for  D^(u,v) 


Thus,  we  can  express  D^(u,v)  as 


(A. 12) 


where  the  function  F{x,y,L)  is  defined  as 

F(x,y,L)~H(x*y,L)-H(.x-y,L) , 

and  tf(x,L)  ia  defined  in  (2.19),  and  f(.xty )  is  defined  in  (3.16).  Using  (A. 12) 
in  (A. A)  to  obtain  D^u,V)  and  substituting  in  (3.20)  to  get  the  coverage  K(X,Y), 
we  obtain 


(A. 13) 


where  p'  is  as  defined  in  (3.19).  The  expression  for  K(.X,Y)  in  (A.  13)  forms  the 
basis  for  £(14)  in  the  machine  program. 

For  a  single  dispenser,  #(!)*!,  we  expand  the  expression  in  (A. 13)  to 


obtain 
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(4.14) 


N'2) 

m.Y)~  l 
k~  1 


The  expression  for  K(X,Y)  in  (4.14)  forms  the  basis  for  X(4)  in  the  machine  pro¬ 
gram.  Again  X(4)  and  X(14)  are  equivalent  but  are  computed  differently. 

4.1.3,  Approximate  Method:  Fragment-sensitive  Target 

In  the  problem  considered  in  Sec.  4.1.1  for  a  fragmentation-type  weapon,  we 
can  simplify  the  expression  for  X(X,Y)  if  the  range  offsets  a.m 0  by  an  approxima- 
tion  for  the  range  ballistic  distribution.  Referring  to  the  expression  for  D^(u,v) 
in  (4.6),  we  have 


«  L(  1)\ 
\g(3) *q(3)  / 
<?(3) 


which  is  a  density  function  for  a  distribution  function.  If  we  replace  this  func¬ 
tion  by  a  uniform  distribution  o(u,£(5))/2L(5) ,  where  of  is  defined  in  (2.9),  we 
obtain  the  approximate  expression  for  D^{utV)i 


(4.15) 


D^(.U,V) 


2£(5) 


ar(u,M5)) 


'  m) 


We  will  demand  that  the  functions  h(u/q(.3)  ,L(l)/<f(3))/q(3)  and  a(u,t(5))/2£(5) 
have  the  same  variance.  Then  L{ 5)  is  determined  by 


(4.16) 


i2(5)-£2(l)+3*2<3). 


Using  (4.15)  in  (4.4)  to  determine  D^(.u,v),  we  obtain 


0£j(u,y)-a(u,£,{5)) 


Nil) 

‘-It 


L  i«l 


2L{5)  7?5T2) 


exp 
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Substituting  for  D^{u,x'j  in  (2.27),  we  obtain  the  coverage  K(X,Y)t 


(4.17) 


K(X,X) 


00  00 

-J  <*(*t(l)+X,£(5))>i(x,^)<£cJ 


N(l)( 

l-  rr  1 


i-rr 

t=i 


,  P^(l)i?(2) 

1  2L(5)/I<j(2) 


Using  the  definition  of  F(x,y,L )  in  (4.12),  we  find  that  the  first  integral  in 
(4.17)  is 


F 


X  U 5)  A(l)\ 

i(l)‘t(l)’t(l)/* 


Since  F(-x,2/,L)«F(x,$/t  L) ,  we  obtain  for  K(X,Y)  from  (4.17) 


(4.18) 


P/1  &(5)  >4 ( 1  )\ P°° 


II  S1  2L(5)72?(T) 


L  i-1 


Equation  (4.18)  forms  the  basis  for  .£(13)  in  the  machine  program. 


If  iV(l)»l,  expanding  the  expression  in  (4.18),  we  obtain 


K(X,Y)-F 


/  X  L(5)  /(l)\Wf2),  n/M/W2)\Fp/S?(l)/?(2)T 

(jnT'tnT'iarj^j  (-l)  (,  * 


Evaluating  the  Integral  using  Lemma  6  of  App.  D,  we  obtain 


(4.19) 


X(£,y)-F 


f  X__  U5 )_  A(  1)\^2) 
^t(I)  »*(!)»*(! 


/W(2)\ 

~p<7RU)R(2)]k 

\  k  ) 

2L(5)/2’(?(2)_ 

v<?(2)  »<?(2)/* 
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Equation  (4.19)  for  K(X,X)  forms  the  basis  for  X(3)  in  the  machine  program.  The 
expressions  X(3)  and  AX13)  are  equivalent  here  for  K  when  different  computational 
methods  are  used. 

4.2.  UNIFORM  DISTRIBUTION  OVER  A  RECTANGLE 

We  consider  a  ripple  of  dispensers  whose  subweapon  distribution  is  uniform 
in  range  and  deflection  over  a  rectangle.  Each  dispenser  is  assumed  subject  to 
ballistic  errors  according  to  a  gaussian  distribution  with  standard  deviations 
s (3)  and  s(4).  The  centers  of  impact  of  ( a.,b . )  refer  to  the  center  point  of 
the  dispenser  pattern.  Each  dispenser  is  assumed  to  release  the  subweapons  so 
that  the  distribution  of  the  subweapon  is  uniform  over  a  rectangle  [21(3) ,2£(4) ]. 

If  the  pattern  center  is  at  (0,0),  the  kill  probability  D^(u,v)  to  a  target  element 
at  (u.u)  due  to  a  specific  weapon  is 

M  3)  £(4)  ,, 

(4.20, 

where  D(u-x,v-y)  is  the  damage  function. 

If  N{ 2)  weapons  are  released  and  if  the  center  point  of  the  pattern  is  at 
(0,0),  the  kill  probability  from  a  dispenser,  D^(.u,v),  is 

(4.21)  DD'(u,v)-l-[l-pDt(w,y))W(2), 

where  D^(m,v)  is  given  by  (4.20),  and  p  is  the  subweapon  reliability  factor. 

The  dispenser  is  subject  to  a  gaussian  ballistic  error,  so  that  D^(.u,v) ,  the 
dispenser  kill  probability,  is  given  by 

(4.22)  Dd(u,V)~  J  J_mDD{ (M~*  ’  V'y )0(eW): B  (3)8(4) 


For  a  fragment-sensitive  target,  substituting  the  damage  function  (2.1)  in 


(4.20),  we  obtain 


Db(u,v) 


M3)  L(  4) 

f  f 

J-M3)J-£.(4) 


dxdy 

4L(3)L(4) 


n/?(l)i?(2),/  U  L( 3)\.  /  V  £(4)\ 
Yd)Y(2r\Y(l),Y(l)r\Y(2),Y(2);* 


where 


Y2(l) 


,/?2d) 

2O0  * 


Y2(2) 


/?2(2) 

"  2Z?0  ' 


The  expression  for  £>^(w,y)  from  (4.22)  is  thus 


(4.W 


r,,vOiV(2)) 


dxdy 

8(3)a(4)- 


a.  Approximation  1:  Ballistic  Error,  No  Edge  Effects.  The  computational 
methods  based  on  (2.27)  will  not  work  if  the  expression  for  D^(u,v)  contains  an 
integral.  Thus,  expression  (4.23)  must  be  simplified.  One  method  is  to  replace 
the  density  functions  in  (4,23): 


VyO)^!?/ 


Y(l) 


and 


v  LW\ 

\y(2)  *  v  ( 2 ) ) 
Y<2) 


by  the  uniform  density  functions 


2L(5) 


and 


a(VtL(6)) 

223(6) 
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where  a(x,f,)  is  as  defined  in  (2.9)  and  variances  are  equal.  Thus,  &( 5)  and  l,(b) 
are  given  by 


(4.24) 


L-  ( 3 )  -I..2  ( 3 )+3 y2  ( 1 ) =L2  ( 3  )4  , 

£2(6)=L2(4)+3y2(2)=/.?(4)4^2)  . 

2Dq 


hasentlally,  we  are  assuming  the  ratios  /?(1)/L(3)  and  /?(2)/£(4)  are  siral  1 .  We 
enlarge  the  pattern  |.2£(3)  ,27.(4)  ]  to  be  I 2L( 5) , 2L(6) )  and  assume  no  edge  effects 
in  the  enlarged  pattern.  Thus,  we  have 


(4.25) 


DD(u, v)= 


V"L  4L(5)£i(6)  j 


_ 1_ _ 

4£(5)L(6) 


If  s(3)  and  8(4)>0  and  we  use  the  definition  of  f(x,y)  in  (3.16),  the  expression 
for  Dq(u,V)  becomes 


(4,26) 


^(u.y)* 


it  ft  \r(u  &(5)W  i>  £(6)> 

\  l  4i(5)£(6)J  y\e(3)t8(3))J\eW,eW)’ 


When  we  use  (4.26)  for  a  ripple  of  N(  1)  dispensers  with  Cl's  at  (a. ,(>.),  the 

v  v 

kill  probability  D^(.a,v)  due  to  the  whole  ripple  is 


N(  1) 

(4.27)  Cp(w,y)»l-  jj  1 l-rP^fu-a^ ,v-b^)] 

i*  1 


'7iff(l)ff(2)l'V(2)l  JU'ai  L(5)\JV~bi  L(6)\ 
4£(5)£(6)J  j’' \s(3)  ’6(3)/' \s (4)  ’8(4)^’ 


where  r  is  the  dispenser  reliability  factor. 
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Substituting  the  expression  for  D  (u,v)  as  given  by  (4.27)  in  (2.27), 


obtain 


(A. 28)  K(  X 


,»-f  f  L  Ji-r,-easnM2)f  <2)l/t<1  t.sA 

j|  I  L  AZr(5)L(6)J  y\  a(3)  *8(3)  J 


m. 

J\  8(4)  ’s(4) 


e(A)  ’« (4)/  l  ^  (^TaT)^ (y cic<^ * 


Equation  (4.28)  forms  the  basis  for  *(103)  in  the  machine  program.  For  only 
one  dispenser,  we  can  integrate  the  expression  in  (4.28),  using  Lemma  6  of  App.  D: 


(4.29)  A(*,Y)»or>|l- j\-i 


4£(5)£(6) 


?(JL.  M5)  4(1  )\J  y  L( 6)  >1  (2 )\ 
\£(3) *t(3)  **(3)^(4) 


where  F(x,y,L )  is  defined  in  (4.12)  and  *(3>  and  t( 4)  are  defined 


(4.30) 


*2(3)=*2(l)+82(3), 

t2(4)-*2(2)+82(4). 


The  expression  for  the  coverage  K(X,Y)  in  (4.29)  is  for  the  case  of  a  single 
dispenser  rectangular  pattern  where  the  ratio  of  the  subweapon  MAE  (mean  effec¬ 
tive  area)  to  the  pattern  area  f4i(3)L(4)J  is  small,  it  can  also  be  considered 
as  an  approximation  to  the  case  discussed  in  Sec.  4.1.3,  "Approximate  Method: 

Fragment-sensitive  Target."  For  the  range  dimension  L(5),  we  use  the  value  ob¬ 
tained  in  (4.16),  i.e.. 


I2  (5)-L2  (i)+3^2  (3)«L2  (l)+3<72  (0)+^lU 

22  o 


and  equivalently 


(4.31) 


£2(3)«L2(l)+3?2(0). 


We  make  a  similar  approximation  in  deflection,  again  equating  variances,  and 


obtain 
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i2(  6)-3<72(2)-382(2)+^p-, 


and  equivalently 


(4.32) 


L2(4)-3a2(2). 


Equation  (3.29)  with  L( 3)  and  L( 4)  given  by  (4.31)  and  (4.32)  form  the  basis  for 
X(10)  in  the  machine  program, 

b.  Approximation  2;  Edge  Effects.  No  Ballistic  Error.  A  second  approxima¬ 
tion  for  D(u,V )  in  (4.23)  may  be  made  if  a(3)  and  a(4)  are  small  with  respect  to 
L(l)  and  L( 2).  We  then  assume  a(3)*8(4)»0,  so  that  Dp(u,V)  is  given  simply  by 


(4.33) 


:2> 


For  a  ripple  of  N(  1)  dispensers,  the  kill  probability  D^[u,v)  due  to  the  whole 
ripple  is 


N(  1) 

(4.34)  Dd(u,v)~ t-  jj 
i- l 


.  f,  I".  vmi)h(.2).(u'ai  LO )\JV~bi  M4)\F(2)l 

!"r \l~ L  7('i'h (TTATnT’TuTATm’^T/J  j  . 


Substituting  (4.34)  In  (2.27),  we  obtain 


r  „,r  n  f  f  V  Nil)  .  j.  T,  pnR(\)R(2)h(Xtil)+X~ai  L( 3)\ 

(4.35)  J-  f[  y(i')V<'2)  j -Y(ir 


x  h( 


rytiD+Y-b,  £(4)\l*(2)g 


Y(2)  ,y(2)/J 

Equation  (4.35)  forms  the  basis  for  Af(104)  ir.  the  machine  program.  The  expressions 
for  K(X,Y)  in  (4.28)  and  (4,35)  are  arswers  to  the  same  problem  but  different  ap¬ 
proximations  are  used;  i.e.,  Af(103)  and  7f(104)  are  equivalent. 


4,2,2,  Impact-sensitive  Target 

For  an  impact-sensitive  target,  by  substituting  the  damage  function  (2.8;  in 
(4.20),  we  obtain  for  D^{u,v) 
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(4.36) 

We  will  assume  that  B(3)/B(1)>1  and  £(4)/B(2)  >1 .  Define  6 (x,B,Z,)  as  the  integral 
of  two  a  functions,  i.e., 


(4.37)  0(s,B,£)  =  J  a(.x-y,B)a(.y,L)dy. 

—  CO 

Thus,  we  have 

/I  \x\<L-B, 

(4.38)  8(x,B,Z;WB+B-x  Is-Bz\x\zb>rB  > 

(o  | x | >L+B . 

The  function  6(x,B,I«)/4B£  is  a  density  function  with  variance  B2+B2.  When  we  use 
(4.37)  and  (4.38)f  D^(.u,v)  in  (4.36)  becomes 

p„8(utB(l),L(3))6(u,B(2),L(4)) 

(4.39)  Db(u,V)*  4B(3)£(4)  * 


For  a  single  dispenser,  D^{u,v)  is  given  by 


(4.40)  DD{utV)m  1-J  J  ^l-p/8(n-x,B(l)  ,B(3) ^ 


N(  2) 


x  (3))  g{fa)lV0w' 
where  p'-pp^  as  in  (3.19). 

a.  Approximation  1:  Ballistic  Error,  No  Edge  F.ffects.  As  for  the  frag¬ 
mentation  case  in  Sec.  4.2.1,  we  must  make  an  approximation  for  D ^  in  (4.40). 

As  one  approximation,  we  will  replace  the  density  function  8(u,fl(l) ,B(3))/4B(1)B(3) 
in  (4.40)  by  the  density  function  or(w,L(5) )/2£(5) ,  again  requiring  equal  variances. 
Thus,  we  let 


-36- 


e(M.fl(l).L(3)).a(«t£(5)) 
4B(1)L(3)  21(5)  * 


where 


(4.42) 


L2(5)-L2(3)+B2(l). 


Similarly,  we  let 


B(i>.B(2).£(4)).tt(i>t£.(6)) 
4B(2)L(4)  2L(6)  * 


where 


(4.44)  £2(6)«£2(4)+B2(2). 

Then,  in  place  of  (4.40),  we  have 


00  oo  p 

(4.45)  Dd(u,v)«  1-J  J  l-p'a(u-x,L{5))a(v-yM(>))j$~$i 


-\N(2) 


'  9{e  (3))  9(^Ui)eO)B%) 


4l-  <2)lf/JL„  M5)\  /  V 

I  L  /•  (5)£(6)  J 


£(6)\ 

8(4)/' 


We  note  that  this  expression  for  D^(u,v)  is  the  same  as  for  a  fragmentation 


weapon  in  (4.26)  except  that  4B(1)B(2),  the  area  of  the  target,  replaces -the 
MAE-nB())B(2) ,  and  the  definitions  of  L(5)  and  £,(6)  in  (4.42)  and  (4.44)  are 
slightly  different  from  those  in  (4.24).  Thus,  the  expression  for  K(X,Y)  is 
similar  to  (4.28): 


(4.46) 


K{x,y) 


p'B(l)B(2) 

£(5)£(6) 


JxtCD+X-aj  l(5)\A 

J\  o(3)  ’e(  3)/;\ 


fyt(2)+X-b. 


8(4) 


ML 

’8(4) 
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Equation  (4.46)  is  the  basis  for  £(105)  In  the  machine  program.  For  a  single 
dispenser,  we  obtain 

(4  47)  m  V)-Jl-  [j  p/5(l)5(2_)7(2)ln/  £,  UV  AO)\  JJ__  UU  MV.) 

14.47)  *(£,£>  jl  L(5)i(6)  J  (F\tO),tO)’tO))F\tW,t(.U),tw)* 


where  t(3)  and  t( 4)  are  as  in  (4.30)  and  L(3)  and  L( 4)  as  in  (4.42)  and  (4.44). 
b.  Approximation  2:  Edge  Effects,  No  Ballistic  Error.  A  second  approxima¬ 


tion  is  to  assume  e(3)«a(4)“0.  Then  from  (4.40)  D^(u,v)  is 


(4.48) 


For  a  ripple  of  JV(1)  dispensers,  we  obtain,  as  in  (4.35), 


r  r  1  *(1) 

(4.4„  i- 

im\ 


L-r  jl-|\~p' 


6(a:i(l)+£-a.,fl(l),L(3))6(j/t(2)+y-i.lB(2),L(4)) 


4L(3)L(4) 


'  W(2)|' 


Equation  (4.49)  forms  the  basis  for  £(106)  in  the  machine  program.  Thus,  expres¬ 
sions  (4.46)  and  (4.49),  i.e.,  £(105)  and  £(106),  are  solutions  to  the  same  problem 
using  different  approximations. 

4.3.  UNIFORM  DISTRIBUTION  OVER  AN  ELLIPSE 

We  consider  the  same  problem  as  in  the  previous  section  except  that  the  sub¬ 
weapons  are  assumed  to  be  uniformly  distributed  over  an  ellipse  with  semiroajor  axis 
T(l)  and  semiminor  axis  T{ 2).  The  kill  probability  D^{u,v)  is 

(4.50)  ^(M”x»v~J/)y(i)y^2)1T* 

C\ 
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where  C\  is  the  area  inside  the  ellipse  [x/Td)  12+[y/2’(2)  ]2»1,  The  kill  prob¬ 
ability  Dd(u,v)  for  a  single  dispenser  against  a  target  element  at  (u,v)  is  thus 
given  by  (4.22)  and  (4,21)  with  D^{u,x>)  from  (4.50). 

4.3.1,  Fragment-sensitive  Target 

For  a  fragment-sensitive  target,  substituting  the  damage  function  (2.1)  in 
(4.50),  we  obtain 


(4.51) 


where 


D^(u,v) 


R(l)R(2) 

T(1)T(2) 


P~M\  dxdy 


y(2)/ y(I)y(2)  ' 


Y2(l) 


_/?2(l) 
2Da  ■' 


Y2(2) 


fl2(  2) 

2 D0  4 


In  Sec.  3.3,  the  .->.1861  ellipse  function  P(.A,B;x,y)  was  defined  in  (3.26)  by 

(*•52)  P(4,B;x*)-JJ  gtt)g{T0dm, 

C2 


where  C2:  [  (C-x)//l]2+L  (Vl/)/S]2^l .  We  can  thus  express  D^{u,v)  in  (4.51)  as 


(4.53) 


D^(u,v) 


/?(!)/?( 2)  mi)  T(Xt  _u _ v_  \ 

rornTr  \y(1),y(2)’y(1)*y(2)7* 


with  P(.A,B;z,y)  as  in  (4.52).  Thus,  Dp(.u,v) ,  which  is  similar  to  (4.23),  is 


(4.54) 


p/?(i)ff(2)rmi)  T(2)  u-x  v-u  \1W(2)) 
J’d)2’(2r  Vy(1)*y(2)’y(1)*y(2);J  ( 


\  dxdy 
)/e  (3)s  (4) ' 
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a.  Approximation  1:  Ballistic  Error,  No  Edge  Kffects.  As  before,  we  muse 
simplify  (A. 54).  If  TXO/yO)  and  T(2)/y(2)  are  large,  we  make  the  following  ap¬ 
proximation  in  (4.53): 


(4.55) 


D^(u,v)> 


(  fld)/?(2) 

\  r(i)?(2) 


r u  v 

■Ml 

[nn 

[_y(2)  j 

otherwise . 


Again,  this  is  the  approximation  often  used  when  we  ignore  edge  effects:  The 
kill  probability  is  the  ratio  of  the  MAF.  to  the  pattern  area  inside  the  pattern 
and  zero  outside.  Under  this  approximation,  we  obtain  from  (4.22) 


dxdu 

(3)e(4)* 


where 


[x-w  I2 

Lm)] 

|/<2>J 

The  integral  in  (4.56)  is  also  an  offset  ellipse  function.  Thus,  we  have 


(4.57) 


pn(u  vJi-r<  em>m1y(2)lp(m  m i-jj_  .jl_\ 
UDKU'V>  y  (_  .'(l)f(2)  J  j  \fl  (3)  *e  (4)' 8  (3)  *8(4)/’ 


where  P(A,B;X,y)  is  defined  in  (4.52). 

For  a  ripple  of  0(1)  dispensers  with  Cl's  at  (a .,£.),  the  kill  probability 

t*  t 

£>p(u,o),  which  is  similar  to  (4.27),  is 


(4.58)  0p(u,t>)-l 


TU- TA2) _u-“i  vsh) 

II  r:  I*  y(l)T(2)J  (  \«(3) *«(4) *a(4) *8(4)/  • 

£■1  L  '  ’  J 


Substituting  the  expression  for  0p(«, u)  given  in  (4.58)  in  (2.27),  we  obtain 
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xpfill) 

\«(3) 


(i)  n 2) 


(3) ’s(4) ' 


xt(l)+£-a. 
8(3) 


equation  (4.59)  is  the  basis  for  .£(107)  in  the  machine  program. 

If  the  pattern  is  circular,  7(  l)*7’(2)*=2’,  and  the  dispenser  ballistic  errors 
are  equal,  s(3)*fl(4)>*fl,  the  expression  for  £(£,y)  in  (4.59)  becomes 


</\xt(l)+X-a. )2+lyt(2)+Y-b .  ]2\ 


where  P(/?,r)  is  the  offset  circle  function  discussed  in  Sec.  3.3  and  defined  in 
(3.28)  as 


(4.61) 


P(/?,r)«  JJ  g{x-r)g(y)dxdy. 

x2+y2sR2 


As  discussed  in  Ref.  1,  the  circular  coverage  function  P(P,r)  may  be  approxi¬ 
mated  for  R  large: 


(4.62) 


P(*,r)*J 

r- 


g(u)du*\-G(r-'ffi-\) . 

//?2- 1 


This  approximation  is  very  close  for  P£5  but  should  not  be  used  for  P< 3.  Using 
this  approximation  in  (4.60),  we  obtain  a  second  approximation: 


(A. 63) 


m,Y)> 


»  .«  (  iv(i)r 
f  f  ;i-  Tfl  i-r> 

1  i~\  L 

(  L  T2  J  ) 

H 


Axt{  1  )+X-a .  l2+lyt(2)+Y-b.  I2-  /t2-b~2 ' 


xK*4of)Hy4m>^- 

b.  Approximation  2:  Edge  Effects,  No  Ballistic  Error.  A  second  approxima¬ 
tion  in  (A. 54)  is  to  set  8(3)**e(4)“0.  Then  D^(u,v)  from  (A.5A)  is 


(A.6A) 


D  (tt  uV<1  T  Pj.(l)/?(2)r/M  _JL-\f(2> 

VU,UJ  1  L1  -2’(1)T(2)^\y(1  ) * y(2)  ,y(I),y(2)/J 


where  y( 1)  and  y(2)  are  as  in  (A. 51).  For  a  ripple  of  N( 1)  dispensers,  Dp(u,v) 
is  thus 


Dp(u,v) 


pR(l)R(2)Jni)  T(2)u~ai  tf'blVly(2)l 

T(l)7’(2)  \Y(1)  *y(2)  *Y<1)  *y(2) /J  j 


Substituting  in  (2.27),  we  obtain 


(4.65) 


K(X,Y) 


pff(l)/?(2) 

T(1)T(2> 


xp 


in  1)  T(2)  »t<2)M-trt-]»<2>( 

\V(1)'Y<2)!  Yd)'  ’■  y(2)  1\  ( 


Equation  (4.65)  is  the  basis  for  ^(108)  in  the  machine  program.  Equations  (4.59) 
and  (4.65)  are  two  approximations  to  the  same  thing  by  different  methods. 

If  the  pattern  is  circular,  7’(1)*7’(2)»!T  and  if  /?( 1  )=/?i2)-f?.„  so  that  y(1)*y(2)= 


Kq.  (4,65)  for  K(X,Y)  becomes 
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(4.66) 


Axt(l)+X-a.  ]^fCyt(2)+y-b .  ]2 

y 


where  P(fl,r)  is  the  offset  circle  function  as  in  (4.61).  Equation  (4.66)  corre¬ 
sponds  to  Eq.  (4.60). 

Finally,  if  we  use  the  approximation  for  P(/?,r)  in  (4.62),  we  obtain  the 
second  approximation,  corresponding  to  (4.63)  for  K(X,Y): 


(4.67) 


K(X,Y) 


A  ^+x-ai i2+\-yt(2H*-bi ]2-/r2-Y2  ^(2)j^ 


4.3.2.  Impact-sensitive  Target 

For  an  impact-sensitive  target,  substituting  (2.8)  in  (4.50),  we  obtain 

(4.68)  £>i ,(».*>)« (2)  1  I  a(w”®>®U))ar(y-j/,S(2))dixij/, 

where  Cu :( x/T( 1 ) J2+[ y/T(2) \?< 1 .  If  T(l)/B{l)  and  f(2)/S(2)  are  large,  we  can 
make  an  approximation  similar  to  that  used  in  (4.55),  i.e.,  let 


(4.69) 


'  4fi(l)P(2) 

\rjk  ttT(  1)7’(2) 


Pp(M, w)= 


.  0 


otherwise. 
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Thus,  If  (4.22)  is  used,  D^{utV)  becomes 


d  tu  _v\ 

VdKU,V)  ^  v?am2)  J  |%(3)»8(4)'a(3)*e(4);' 

This  is  exactly  the  same  as  D^(u,v)  in  (4.57)  for  the  fragment-sensitive  case, 
except  that  we  have  replaced  ft(l)R(2) /T(l)?(2)  by  the  ratio  4B(l)B(2)/7t!T(l)7’(2) 
and  p  by  p'.  Thus,  we  obtain  an  expression  that  corresponds  to  (4.59): 

“  (  iV(Df  {  r 

«.70) 


Ad)  mi? 

AflO)^^)*- 


zt(l)+X-a.  yt(2)+Y-b. 


'  «  V  '  '  '  4\ 

S(3)  *  TUT)  7 


Equation  (4.70)  forms  the  basis  for  Af(109)  in  the  machine  program.  For  the  cir¬ 
cular  case,  Eq.  (4.60)  applies  with  the  above  replacement,  i.e., 


(4.7., 

f  im l  L  ' 


Kxt  ( 1  )+X-a.  ]  2+f  yt  ( 2  )+Y-b .  J  2 


Finally,  for  T/s  sufficiently  large,  we  obtain  an  expression  similar  to  (4.63); 


(4.72) 


(A, xt ( 1 )  +X-a.  J 2+r wt ( 2 )  +Y-b .  ]2_/r2.82  'i 
- "  8 - - - / 
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We  note  that  the  approximation  ignores  the  edge  effects.  If  the  edge  effects 
are  important,  i.e.,  r(l)/B(l)  and  T(2)/B(2)  small,  we  may  transform  the  elliptical 
pattern  into  an  equal  area  rectangular  pattern  and  use  the  results  of  Sec.  4.2.2. 

In  this  case,  we  would  let 


(4.73) 


L(  3) 


U  4) 


/ng’(2) 

2  * 


and  use  Eq.  (4.49). 


4.4.  UNIFORM  DISTRIBUTION  OVER  AN  ELLIPTIC  ANNULUS 

We  consider  the  same  problem  as  that  of  Sec.  4.3,  "Uniform  Distribution  over 
an  Ellipse,"  except  that  the  subweapons  are  assumed  to  be  uniformly  distributed 
over  an  elliptic  annulus,  the  area  included  between  the  inner  ellipse  with  semi¬ 
axes  VO)  and  V( 2)  and  the  outer  ellipse  of  seraiaxes  TO)  and  T( 2).  The  kill  prob¬ 
ability  D^(u, v)  is  now 


(4.74) 


n  (u  ul*  ..  r  m.K.  a.'^nra - 

VU,U'  J.L  nCr(l)2,(2)-fc7(l)l/<2)  3* 

^5 


D(u-x,v-y)dxdy 


where  Cj  is  the  area  between  the  ellipses  [x/r(l) J2+[j//jT(2)]2»1  and  [x/V**(l)l+ 

]«l.  The  kill  probability  Dp(u,V)  for  a  single  dispenser  against  a  target 
element  at  («, t>)  is  thus  given  by  (4.22)  and  (4.21)  with  D^{u,v)  now  given  by 

(4.74). 


4.4,1,  Fragment-sensitive  Target 

For  a  fragment-sensitive  target,  substituting  the  damage  function  (2.1)  in 

(4.74),  we  obtain  for  P^(.u,v)  an  expression  similar  to  that  of  Eq.  (4.53): 


Db(u,v) 


Rj  l)/?(2) 

=7’(1)7’(2)-V(1)V(2) 


[i 


T (I)  T{ 2)  u  v 
Y(1),y(2);y(1)>y(2) 


(/(I)  V(2)  u  V  \] 
Y(1),Y(2)‘Y(1)’Y(2);J’ 
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where  P(A,B;u,v)  is  the  offset  ellipse  function  in  (3.26)  and  y(1)  and  y(2)  are 
as  in  (4.51),  i.e., 


Y2(l) 


Y2(2) 


X2(  2) 

2  D0  * 


Thus,  as  with  Eq.  (4.54),  D^{u,v)  i3  now 


(4.75) 


Dd(u,v ) 


pR(l)B(2) 

T(l)T(2)-V(l)f/(2) 


2'(1)  T{ 2)  u  V  \ 
,Y(1)’Y(2);Y(1),Y(2)/ 


t/(D  V(2)  u  V  \ 
1Y(1),Y(2)’y(1),Y(2)/ 


dxdy 

e(3)a(4) * 


a.  Approximation  1:  Ballistic  Error,  No  Edge  Effects.  If  2’(1)/y(1)» 
J’(2)/y(2),  V(1)/y(1) ,  J/(2)/y(2)  are  large,  we  make  the  same  approximation,  as 
in  (4.55),  that 


(4.76) 


Db(u,v ) 


ff(l)/?(2) 

y(l)J'(2)-t/(l)t/(2) 


if  (u,t>)  is  inside  the  elliptic  annulus  C 5  and  otherwise  zero.  Under  this  approx¬ 
imation,  Dq(u,v)  becomes 


(4.77) 


Dp(.u,v) 


t  p/?d)/?(2) _ 

2’(l)2’(2)-l/(l)t/(2) 


dsdu 

8(3  )8d)' 


The  Integral  in  (4.77)  is  the  difference  of  two  offset  ellipse  functions;  thus 


(4.78) 


0p(u,v) 


.u,. 


p/?(  1  )/?(2) 


l//(2) 


2’(1)J’(2)-V(l)f/(2)J 


(Ki 


TU2  ^(2)  u  U  \  p/t/d)  >/(2) .  » _ y_ 

l8(3),8(4)'8(3)’8(4)/l'/As(3)»«(4)*«(3)’0(4) 


1  ii  this  equal  ion,  ;  t/1  ;l< , is  t  lit*  offset  ellipse  fund  ion  given  in  Kq.  (3.26). 

it  we  proceed  as  in  (.i.V>),  lor  a  ripple  of  A*(l)  di-.per.sers  with  (M's  at 
Uj.  ,;  .'  tnd  a  dispenser  reliability  factor  >•,  the  fractional  coverage  K(X , >’)  is 


a) 


iilUAA’  \ . .. 

."( 1  )'”td’l  •»■(!)!>■(’) 


’.■cm 


.:•(  (  II-*.'-.}  . 

_  l 

" '  7<  n 


y_j. 

■ ' '  7r*T 


\  ./;»•  ti)  h’l-’) 
)  '  \?(  !> ’s(.f) 


(D+X-11. 
_  J. 
"Vt  3)" 


yt(2)+y-b 
a  (!>) 


I.quation  ('>.79)  is  the  basis  for  ,Y(110)  in  the  machine  program.  It  allows  a 
ballistic  error  hut  considers  no  edge  effects.  Thus,  to  he  valid,  the  ratios 
.”( 1  )/£(  l)  and  ?(2)/R(2)  should  he  reasonahlv  large,  i.e.,  •'3. 

If  the  pattern  is  a  circular  annulus,  T(  1  )-'/’(2)  =  r  and  W(  1  )=I*,(2)=iv',  and  the 
ballistic  errors  are  equal,  s  (!)“.>?( '>)=a ,  the  expression  K(X,Y),  which  is  similar 
to  (•'*.(>()),  becomes 

>V(  I ) 

•If 

£=i 


(4.8o)  A(A',n=; 


.  'ixft  ( 1 )  +X-Q;  1  7+]yl(2)  ry-.!. . \ 


-Y 


4xf  ( 1  )+X-a . '  ?-4 !  ut  ( 2  )+Y-l  .  l 


f(** 


A(\J 

t(l) 


)h(*'  irioW^ 


where  Yd-,”)  is  the  oifsel  >  ircle  f 
•\s  a  second  approx i ma »  ion  for 


tine  lion  as  in  (4  61). 

T/s  and  W/n  large,  using  tite  approximation  for 


r)  in  ('*.62),  we  obtain  an  expression  that  is  similar  to  (4.63): 
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Equation  (4.81)  is  the  basis  for  *(100)  in  the  machine  program.  Again,  no  edge 
effects  are  considered  and  min  ( T/R(l) ,T/R(2)  I  should  be  greater  than  5.  Further, 
in  order  that  approximation  (4.62)  be  valid,  27 j  and  W/e  should  be  greater  than  3. 

b.  Approximation  2;  Edge  Effects,  No  Ballistic  Error.  The  second  approxi¬ 
mation  is  to  set  8(3)«a(4)-0.  Then,  proceeding  as  in  Sec.  4.3.1a  for  Eqs.  (4.64) 


and  (4.65),  we  obtain  for  *(*,*)  an  equation  that  corresponds  to  (4.65): 


(4.82) 


Sflmi  m 

[P\Y(l)*r(2)‘  Y(l)  ’  y(2)  ) 


tt(l)+X-a.  yt{2)+Y-bi 

yTT)  •  7(2)  ' 


■)]i  ]  }V  '7$)^ 


» 


where  y(1)  and  y(2)  are  defined  in  Eq.  (4.51).  Equation  (4.82)  forms  the  basis 
for  *(111)  in  the  machine  program.  Equation  (4.82)  is  used  when  the  edge  effects 
are  important,  while  Eq.  (4.79)  is  used  when  the  ballistic  errors  are  important. 
If  the  pattern  is  a  circular  annuluB  2'(1)«2’(2)«2’  and  V(  1 )  —1/(2' we  obtain 

s 

'  an  expression  that  corresponds  to  (4.66): 
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CO  t 

(4.83)  K(X,y)=J  J 


_  on  _  a» 


(V(l)/  ( 

1-  rf  1-r  1- 

Li=i  ' 


L-EfiL 
L  T?-W2 


■14 


f,p  Axt{[)±X-ai  ]2+iyt(2)+Y-bi 


A,  Axt(\)+X-a.  \2+\ut(2)+Y-bj  l?\P 

m  2)h 

hL  40)\, 

A(2)\ 

\Y*  Y  )\\ 

y 

Finally,  a  second  approximation  for  T/y  and  W/y  large  that  is  similar  to 
(4.67)  is 


00  cc 

(4.84)  K(Jr,Y)«J  f 


Nil)  / 


L  i-1 


I'T 


!_-ESL 

t2~w2 


■H- 


/( *C ( 1  )+#-a.  j *+Cy i(2)  +¥-b .  J 2- Aw2- y2 


-G 


( Axt  ( 1  h-X-a^+lyt  (2)+y-bi  32- M-y‘ 


l 


N(  2) 


Equation  (4.84)  forms  the  basis  for  #(101)  in  the  machine  program.  In  this  case, 
T/R  and  W/R  should  be  larger  than  5, 

4.4,2.  tmpact-sensitlv e_Ta  rjje  t 

Fo.  <n  impact-sensitive  target,  D^(m,p)  is  the  same  as  in  (4.68)  for  the  el¬ 
lipse,  except  that  the  elliptic  area  C4  is  replaced  by  the  elliptic  annulus  C ■,  of 
(4.74),  if  T(l)/Z?(l),  T(2) /B(2) ,  W(l)/B(l),  and  (7(2 )/0(2)  are  sufficiently  large, 
we  make  the  same  approximation  as  in  (4.69)  over  C 5,  i.e., 
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Db(u,v )» 


p-4B(l)B(2) 

iO'om^Pvo) v{2)  ]  if  (K,y)  is  inside  c5. 


otherwise. 


U.\der  this  approximation,  K(.X,X) ,  which  is  similar  to  (4.70),  is 


(4.85)  K(X,T) 


.f  r  _ p'-u-b^b-& _ p 

J  J  I  (|  )  ir[2*(l)2’(2)-l7(l)W(2))j 

"  "  \  i-1  '  L  J 

_/t(  1 )  r(2)  y*(2)+Y-2>A 

[V(3)'fl(4l-  s(3)  *•  s(4)  / 


TU)  •  i(4)  / 


.feiimi. 

\a(3) *a(4) * 


xt(l )+X-a.  yt(2)+l’-b> 

V  u 


8(3)  ’  8(4) 


We  note  that  (4.85)  is  identical  with  (4.79)  for  the  fragment-sensitive  target 
except  that  the  ratio  7?(l)/?(2)/[!P(l)r(2)-V(l)V(2)  ]  is  replaced  by  4fl(l)B(2)/ 
ir[2’(l)2’(2)-W(l)V(2)].  Equation  (4.85)  is  the  basis  for  Y(112)  in  the  machine 


program. 


For  the  circular  case,  2!(1)«T(2)-Z’  and  V(1)*V(2)-V,  and  the  ballistic  errors 
*(3)*«(4 )"B,  the  expression  K(,X,X),  which  corresponds  to  (4.81),  becomes 


(4.86)  K(X,X ) 


V  f  1  "  \  * 

•LLrff 


P  LrL  rI.pi^ai^.2,?1jy(2) 

il'j  (  L  n(T2-V2)  J 


(  ( 1  )+X-Si ]2+Lyt (2)+*-^ ] 2- /W2~e2  ^ 


/[«t(l)+^a£)2+CJ/t(2)+T-fc.32-^ 


KX,^T )M*4(!)-p* 


Again,  we  note  that  (4.81)  and  (4.86)  are  identical  except  for  the  factors 


and  mm>, 

k(,T2-V2)  T2-V 2 
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Equation  (A. 86)  forms  the  basis  for  X(10Z)  in  the  machine  program.  Since  no  edge 
effects  are  considered,  27B(1)  and  T/B( 2)  should  be  greater  chan  5;  further,  in 
order  that  approximation  (A. 62)  be  valid,  Tie  and  f//fl  should  be  greater  than  3. 

If  a  case  occurs  in  which  edge  effects  are  important,  i.e.,  T/B(l)  or  T/B{2), 
some  other  method  should  be  used,  e.g.»  X(109). 
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5.  THE  COMPUTER  PROGRAM 

The  simplified  weapons  evaluation  model  described  here  is  sufficiently  broyd 
to  cover  almost  all  problems  that  arise  in  the  field  of  nonnuclear  weapons  evalua¬ 
tion.  The  FORTRAN  program  is  designed  to  compute  the  solutions  to  most  of  the 
effectiveness  problems  discussed  in  Secs.  3  and  4.  After  describing  the  inputs 
and  computational  methods  that  are  generally  used  for  all  problem  types,  we  con¬ 
sider  each  problem  type  that  involves  different  inputs  or  different  methods. 

This  section  has  been  planned  to  permit  the  user  to  run  the  program  without 
referring  to  previous  sections.  In  the  event  that  further  investigation  is  desired, 
however,  the  appropriate  references  are  given  for  each  case  discussed. 

5.1.  GENERAL  CONSIDERATIONS 

The  following  list  provides  the  inputs  required  for  the  program,  with  ref¬ 
erences  to  the  sections  in  which  these  inputs,  their  uses  and  restrictions are 
described.  Odd  integers  in  the  inputs  denote  a  parameter  in  the  x  (range)  direc¬ 
tion;  even  integers  denote  a  parameter  in  the  y  (deflection)  direction. 


FORTRAN 

Inputs  Description  Reference 


A3,A4  Target  area  dimensions  5.1.1 

B3,B4  Impact-type  target  dimensions  5.1.5 

B1,B2  Fragment- type  target  parameters  5.1.5 

Dl,D2  Fragment-type  target  parameters  5.1.5 

CLl,Q0  Dispenser  range  ballistic  parameters  5.3 

GCL3,GCLU  Rectangular  pattern  half  dimensions  5.4 

ETl,ET2  Elliptic  pattern  semiaxes  5.4 

Wl,W2  Elliptic  annulus  inner  semiaxes  5.4 

Ul,U2  Number  of  integration  steps  5.1.4 

5Bl,5i/2  Target  offsets  5.1-1 

51,52  Ballistic  error  standard  deviations  5.2 

53,54  Ballistic  error  standard  deviations 

for  dispenser  5.4 

Tl, T2  Aiming  error  standard  deviations  5.1.2 

A(J) ,B(J)  Individual  weapon  aiming  point  5.1.3 

HXAJ,NXBJ  Flags  for  aiming  point  pattern  5.1.3 

D,DF  Spacing  between  successive  weapons  5.4 

N,NB  Number  of  weapons,  number  of 

subweapons  5.4 

5P  Probability  of  kill  if  hit, 

weapon  reliability  5.1.5 

R  Dispenser  reliability  5.4 
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5.1.1.  Target  Area:  Inputs  A3,  A4,  SUlt  SU2 

The  target  area  in  all  cases  considered  is  a  rectangle  with  dimensions  /13  in 
the  x  (range)  direction  and  44  in  the  y  (deflection)  direction.  The  inputs  43 
and/or  44  may  be  zero.  The  center  of  the  target  area  may  be  offset  from  the  aim 
point  at  SUl,  SV 2,  whore  SUl  is  the  range  offset  and  SU2  is  the  deflection  offset. 

5.1.2.  Aiming  Error:  Inputs  7T ,  T2 

The  aiming  error  distribution  is  assumed  gaussian  (see  Sec.  3.2.1)  with  stan¬ 
dard  deviations  2T  and  T2  in  x  and  y,  respectively.  In  terms  of  REP  and  DEP, 

AREP=. 674421,  ADEP-. 674422.  If  21-22-2’,  then  CEP-1. 17742. 

5.1.3.  Aiming  Point  Array:  Inputs  4(«7),  B(J) ,  D,  DF,  N 

The  aiming  pattern  is  the  array  of  the  K  desired  center  of  impacts  [4(«7)  ,B(«T)], 
,7-1,2,. . .  ,N.  (See  Sec.  3.1.)  The  points  are  so  defined  that  ZA  («7)-0,EB(,7)«0.  The 
pattern  may  be  a  direct  input  or  may  be  calculated.  If  the  flag  MTAt7»0,  a  direct 
input  for  the  A(«/)'s  is  called  for  and  D  must  be  1;  if  the  flag  NXBJ- 0,  the  B(J)’s 
are  called  for  and  DF  must  be  1.  If  NXAJml,  the  4(J)'s  are  calculated,  using  the 
uniform  spacing  D;  if  NXBJ~lt  the  B(«7)’ s  are  calculated,  using  the  spacing  DF . 

5.1.4.  Integration  Routine:  Inputs  V 1,  U2 

The  basic  integration  problem  (see  Sec.  2.3)  is  the  evaluation  of  the  double 
integral 

(5.D  i-J  J  )h(v * 

where  K(xty)  is  a  nonnegative  integrable  function  never  larger  than  1;  4(1)  and 
4(2)  are  the  target  half  dimensions;  and  h(x,A)  is  the  function  defined  in  Eq. 

(2.19)  as 

i  / 

(5.2)  /i(*,4)«^j  g(a*z)dz. 


where  g(x)"  exp(-x?/2)//2:r. 
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We  note  the  limits 


(5.3) 


lim 
T-*  o 


h(x/T,A/T) 

T 


a(x,A)  $2A 

u  S 


if  \x\</ . 
ir  irj>/5, 


lim  h(xsA)~g(x) . 
A-*  0 


In  order  to  une  a  numerical  integration 
grain  in  (5.1)  through  an  appropriate  choice 
E( 2) ,  obtaining  the  expression 


routine,  we  meat  truncate  .•>«  inte- 
oC  the  truncation  limits  £(1)  and 


(5.4) 


.  £(DJ7f2) 
~E(.\y-E(2) 


where  s,  is  the  truncation  error.  To  obtain  £(1)  and  E( 2)  we  consider  in  turn 
two  possible  choices  for  truncation  limits.  The  first  set,  £*(3)  and  l'{ 4),  are 
based  on  the  fact  that  K(x,y)< 1.  For  truncation  limits  £(3)  and  E( 4),  the  corre¬ 
sponding  truncation  error  e2  is 


(5.5) 


where  H(r.,A)  is  defined  as 


x 

(5*6)  ff(x,A)-(’  h(y,A)dy. 

and  is  evaluated  by  Eq.  (2.20).  We  determine  £( 3)  and  J?(4)  satisfying  the 
inequalities 
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Thua  che  truncation  error  e2  is 


(5.7) 


r.  <4(10)-8. 


The  second  set  £2( 5)  and  E( 6)  is  based  on  the  fact  that 

CD 

f  h(.x,A)dx*l . 

^  —  CD 

Thus  for  truncation  limits  27(5)  and  27(6) ,  the  truncation  error  e3 
<5-8) 

We  determine  £2(5)  and  E( 6)  for  eacn  function  K(x,y)  sr.tisfying  the  inequalities 

|x|S(5)X<X^)il0"4> 

(5.9) 

max  *»  /  _  v  i  «**4 

|i/|>£'(6)X(x^)i10  » 

so  that  the  truncation  error  e3  is 

e3s4(l0)"4. 


The  truncation  limits  2(1)  and  £2(2)  are  chosen  as 


(5.10) 


2(l)-rain(£2(3),£2(5)), 
£2(2)*min(2(4)  ,2(6)) . 


The  truncation  error  ci  is  thus 
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1 


if 


When  we  use  the  expression  (5.6)  for  //(/,' i),  wo  obtain  for  1 


(5.15) 


21/(1)  21/(2) 


-[»(»  A*  («-eTi’5#i')]  • 


We  define  L’ne  function  F(x%y,L)  as 


(5.16) 


F(x,y,[j)=H (x+^  ,L)-I1  ( x-y  ,L)  . 


Then,  I  from  Kq.  (5.15)  is 


21/(1)  ?.!/( 2) 


(5-17)  I-  /f(x(j),J/(k))F(x(J),^,^)F(y(/c),.€m,|I|I)i 


In  terms  oi  he  l-'OKTRAN  symbols,  the  integral  J  is 


2U\  ?.U'> 

(5.18)  1=  lKb{x(,j)  ,y(k))Ft  (x(j) 


If  K(z,y )  is  symmetric  in  either  x  or  i/,  the  sum  in  (5.17)  need  only  be  taken  to 
1/(1)  or  1/(2).  If  symmetric  in  both  x  and  t/,  we  obtain 


1/(1)  l'(2) 


(5.19)  I- 


%?,  » F(xU) •fT)--yni)r(y(J:) 


Kquations  (5.17)  and  i5.19)  are  the  basis  for  the  integration  routine,  subroutine 
/150.  for  one  dimension,  the  expressions  are  single  sums  and  are  the  basis  for 


the  single  integration  routines. 
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5.1.5.  Damage  Functions:  Inputs  SP,  B 3,  B4,  SI,  S2,  D 1,  02 

Two  types  of  damage  functions  are  used  (see  Sec.  2.1).  We  define  an  impact - 
sensitive  target  as  one  for  which  there  is  a  definite  geometric  figure  (in  our 
case  a  rectangle  of  dimensions  S3,  B4),  which  must  be  impacted  by  the  weapon  or 
subweapon.  A  target  for  which  a  vulnerable  area  is  given  thus  belongs  in  this 
category.  The  input  parameter  Sb  is  the  probability  of  damage  if  hit.  When  we 
define  SI"  B3/2,B2'»B4/2,  0  is  given  by 

(5  20)  \sp  if-fllSxSBl,-B2%SS2, 

(O  otherwise. 

We  note  that  S3  and  B4  are  the  dimensions  in  the  ground  plane. 

We  define  a  fragment-sensitive  target  as  one  for  which  a  significant  portion 
of  the  damage  effect  is  due  to  the  fragmentation  of  the  weapon.  Usually,  the 
damage  function  is  determined  as  a  damage  matrix  using  a  "lethal  area"  program. 

This  damage  function  is  approximated  by  an  exponential  function  (see  Sec.  2.1)  of 
the  form 

(5.21)  D(x,y)-SP-M'D2  exp{-— +^^-Y 

\Sl2  S22  / 

The  parameters  SI,  02,  SI,  and  R2  are  Inputs  to  be  obtained  from  empirical  data. 
First,  we  require  that  tiS1',B2=MAF,  where  the  MAE  has  been  determined  by  some  other 
method.  The  ratio  of  SI  to  S2,  or  the  ellipticity  of  the  damage  function,  may 
be  calculated  or  estimated.  The  initial  values  0 1  3nd  02  will  usually  be  set  equal 
to  one.  In  this  case  the  input  parameter  SP  is  the  reliability  factor  for  the  in¬ 
dividual  weapon. 

5.2.  RIPPLE  OK  BOMBS:  INPUTS  SI,  S2 ,  N 

We  consider  a  ripple  of  N  weapons  delivered  with  a  common  aiming  error  ac¬ 
cording  to  an  aiming  point  pattern.  Each  weapon  is  nubject  individually  to  a 
ballistic  error  according  to  a  gaussian  distribution  with  standard  deviations 


51  and  52.  The  Inputs  51  and  S?.  may  be  zero.  Again,  in  terms  of  REP  and  DEP, 
BREP-. 674451,  BDEP-. 674452. 

Figure  1  presents  a  flow  diagram  of  the  computer  program  for  the  ripple  of 
bombs.  The  fractional  coverages  2f(l),  Af(ll),  a'(5),  and  Af (15)  are  answers  to  the 
problems  considered  in  Sec.  3.  The  program  is  described  as  follows: 

•  The  fractional  coverage  for  a  ripple  of  bombs  against  a  fragment- sensitive 
target  is  A- ( 1 1 ) .  Pertinent  inputs  are  PI,  D2,  51,  52,  51,  and  52,  which 
are  obtained  from  empirical  data.  The  program  uses  Eq.  (3.7). 


Fig.  1 — Ripple  of  Bombs  Delivered  witn  a  Common  Aiming  Error 
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•  The  fractional  coverage  against  an  iirpaot-oeneitive  target  is  *(15>.  This 
problem  uses  B 3  and  S4  in  the  computations  rather  than  Rl ,  R 2,  D 1,  and  D 2. 
The  program  uses  Eq.  (3,20). 

•  The  fractional  coverages  *(i)  and  *(11)  are  answers  to  the  same  problem 
when  different  but  equivalent  expressions  are  used;  *(5)  and  *(15)  are 
also  equivalent.  However,  *(1)  and  X(5)  are  restricted  to  cases  in  which 
the  number  of  weapons,  N,  is  less  than  or  equal  to  12.  The  computation 
for  N  very  small  is  faster  and  more  accurate  for  X(l)  and  *(5)  than  it  is 
for  *(11)  and  *(15). 

•  The  flag  NO  is  used  to  skip  parts  of  the  program  to  compute  *(1)  or  X(5). 

The  code  is  restricted  to  the  computation  of  X(l),  X(2) ,  *(3) ,  X (4),  or 
X(5)  when  N0»0;  when  N0=1,  *(1)  through  *(5)  are  omitted;  there  are  no  re¬ 
strictions  when  N0«2.  The  subroutine  XR0  computes  *(I)  through  *(3). 
Equation  (3.13)  is  used  to  compute  7(1)  and  Eq.  (3.21)  to  compute  X(5). 

a  The  basic  flow  of  the  computer  model  is  the  same  for  *(11)  and  *(15).  In 
subroutine  PR,  the  P(x,j )  function  has  been  redefined  for  *(15).  The  flag 
K0UT,  which  is  automatically  set  and  is  not  a  data  input,  determines  which 
P(x,j)  function  to  use. 

5.3.  RIPPLE  OF  FIXED  DISPENSERS:  INPUTS  NB,  S2,  TABLE.  OR  CL\ ,  Q0 

We  consider  a  ripple  of  N  fixed  dispensers  actuated  with  a  common  aiming  error 
according  to  an  aiming  point  pattern.  From  each  dispenser,  NB  subweapons  are  re¬ 
leased.  In  range,  the  ballistic  dispersion  for  the  subweapons  within  each  dispenser 
is  given  by  an  empirical  table.  In  deflection,  each  subweapon  is  assumed  subject 
to  a  gaussian  ballistic  error  with  standard  deviation  S2.  The  range  pattern  is 
fitted  by  an  approximating  function,  the  stick  distribution  H(x/Q0 ,CL\ /Q0) ,  dis¬ 
cussed  in  Sec.  2.2.2c  and  App.  B,  through  the  use  of  the  two  parameters  CL\  snd  <?0. 

Figure  2  presents  a  flow  diagram  of  the  computer  program  for  a  ripple  of 
fixed  dispensers  The  fractional  coverages  *(2),  *(12),  *(3),  Y(13),  *(4),  and 
*(14)  are  answers  co  the  problems  in  Sec.  4.1. 
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The  program  is  described  as  follows: 

•  The  fractional  coverage  for  a  ripple  of  fixed  dispensers  agai.ist  a  fragment- 
aeneitive  target  is  *(12).  The  inputs  are  R 1,  R 2,  01,  and  02;  the  appro¬ 
priate  reference  is  Eq.  (A. 9). 

•  The  approximation  to  *(12)  is  *(13).  The  computation  involves  only  a 
single  integration.  It  is  thus  faster  than  the  computation  for  *(12),  but 
a  degree  of  accuracy  is  lost.  The  appropriate  reference  is  Eq.  (4.18). 

e  The  fractional  coverage  for  an  irrpaat-censitive  target  is  *(14).  The 
Inputs  03  and  04  are  necessary.  The  appropriate  reference  is  Eq.  (4.13). 

•  The  fractional  coverages  *(12),  *(13),  and  *(14)  use  Oil  and  QO,  which 
are  dispenser  range  ballistic  parameters.  If  the  ballistic  data  are  'n 
tabular  form,  i.e.,  CtfPT-1,  and  the  table  is  read  an  an  input  array,  3et 
001-1  and  Q0-0.  The  actual  values  will  then  be  computed. 

•  The  fractional  coverages  *(2),  *(3) ,  and  *(4)  are  the  same  as  *(12),  *(13), 
and  *(14);  however,  they  are  restricted  to  .V-l  and  OOS30. 

5.4.  RIPPLE  OF  PATTERNS :  INPUTS  S3,  54,  00,  R 

We  consider  a  ripple  of  N  patterns,  each  containing  00  subweapons,  delivered 
with  a  comuon  aiming  error  according  to  an  aiming  point  array  for  the  patterns. 

The  A(J),  B(J)  of  Sec.  5.1.3  are  now  the  pattern  aiming  points.  Each  pattern  is 
subject  to  a  gausslan  ballistic  error  with  standard  deviations  53  and  54  and  to  a 
reliability  factor  0.  The  NB  subweepons  are  assumed  to  be  uniformly  distributed 
within  the  pattern.  In  the  computations  we  had  to  make  one  of  two  assumptions. 

The' first  was  that  the  damage  level  was  uniform  throughout  the  pattern.  This  as¬ 
sumption  Ignores  the  edge  effects  and  is  valid  whenever  the  pattern  area  is  much 
larger  (say,  16/1)  than  the  target  area  4«S1»S2  or  the  MAE-it W?l*/?2.  Cases  for  which 
this  assumption  is  made  will  be  designated  as  "no  edge  effects."  For  some  cases, 

the  edge  effects  may  be  considered  if  the  ballistic  error  standard  deviations  53  and 

♦ 

#  54  are  ignored.  These  cases  will  be  designated  as  "edge  effects." 
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5.4.1.  Rectangular  Patterns:  GCL3t  GCU 

Consider  a  ripple  of  N  rectangular  patterns  of  dimensions  2GCL3t  2GCLb,  de¬ 
livered  with  a  common  aiming  error  according  to  an  aiming  point  array. 

Basically,  the  flow  diagram  for  a  ripple  of  patterns  is  the  same  as  that  for 
a  ripple  of  fixed  dispensers.  (See  Sec.  5.3.)  The  coverage  *(103)  for  a  fragment- 
sensitive  target,  ignoring  edge  effects,  is  computed  according  to  .'.4.  (A. 28).  The 
inputs,  which  are  derived  from  empirical  data,  are  GCL'i  and  GCL4,  which  are  the 
rectangular  pattern  dimensions;  51,  R2,  01,  D2,  R;  and  53,54*0.  The  coverage  ( 103 ) 
should  be  restricted  to  the  cases  in  which  GCL3/RI  and  GCLU/R2  are  greater  than  5. 

If  these  ratios  are  not  greater  than  5,  the  coverage  *(104)  for  a  fragment-sensitive 
target,  including  edge  effects  (i.e.,  53=54*0),  should  be  used,  The  coverage  *(104) 
is  computed  according  to  Eq.  (4,35). 

The  coverage  *(105)  for  an  impact-sensitive  target,  ignoring  edge  effects, 
requires  inputs  53,  54,  and  53,54*0.  The  ratios  2GCL3/B3  or  2G5L4/B4  must  be 
greater  than  5.  If  they  are  less  than  5,  the  coverage  *(106)  for  an  impact- 
sensitive  target,  including  edge  effects,  should  be  used.  The  coverage  *(105) 
is  computed  according  to  Eq.  (4.46)  and  *(106)  according  to  Eq.  (4.49). 

5.4.2.  Ripple  of  Elliptic  Patterns 

The  flow  diagram  for  the  following  three  cases  is  the  same  as  that  for  the 
corresponding  rectangular  cases.  The  coverage  *(107)  for  a  ripple  of  N  elliptic 
patterns  against  an  area  of  fragment-sensitive  targets,  ignoring  edge  effects,  is 
computed  according  to  Eq.  (4.59).  The  coverage  *(107)  should  not.  be  used  if  the 
ratios  £71/51  or  ET2/R2  are  less  than  5.  The  coverage  *(108)  for  a  ripple  of  N 
elliptic  patterns  against  an  area  of  fragment-sensitive  targets,  including  edge 
effects,  should  be  used  instead.  The  coverage  *(108)  is  computed  according  to 
Eq.  (4.65). 

The  coverage  *(109)  against  an  area  of  impact-sensitive  targets,  ignoring  edge 
effects,  is  computed  according  to  Eq.  (4.70).  The  coverage  *(109)  should  not  be 
used  if  the  ratios  £71/51  or  ET2/B2  are  less  than  5;  53  and  54  should  not  be  zero. 
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5.4.3.  Elliptic  Annulus  Patterns:  ET 1,  ET 2,  VI,  V2 ,  fl 

We  consider  a  ripple  of  V  elliptic  annulus  patterns  that  are  delivered  with 
a  common  aiming  error  according  to  an  aiming  point  array.  The  axes  of  the  outer 
ellipses  are  ET i  and  ET2,  while  those  for  the  inner  ellipse  are  VI  and  V2.  The 
computation  for  the  elliptic  patterns  involves  the  elliptic  coverage  function, 
which  is  time  consuming.  In  each  of  the  following  cases,  a  circular  approximation 
is  given,  a t;  approach  that  is  fast  and  in  many  cases  gives  a  very  close  approxi¬ 
mation.  Further,  if  we  are  dealing  with  circular  patterns,  we  may  use  these 
approximations  in  place  of  the  elliptic  patterns  of  the  previous  section  by  setting 
the  inner  radius  to  zero. 

The  coverage  *(110)  for  a  ripple  of  N  elliptic  patterns  against  an  area  of 
fragment-sensitive  targets,  ignoring  edge  effects,  is  computed  according  to  Eq. 

(4.79)  if  the  flag  AItlO-1.  If  the  ratios  ETl/Rl.  ET2/P2,  Vl/ffl,  or  W2/F2  are  less 
than  5,  it  is  preferable  to  use  *(111),  which  takes  edge  effects  into  account. 

For  the  special  case  in  which  the  ellipses  are  circles,  /?l*/?2  and  S3*S4,  the 
special  approximation  may  be  used  if  the  ratios  ETl/S'i  and  V1/S3  are  larger  than 
3.  In  this  case,  the  coverage  *(100)  for  a  ripple  of  N  circular  annulus  patterns, 
ignoring  edge  effects,  is  computed  according  to  Eq.  (4.81)  if  the  flag  IV10CM. 

Again,  we  may  use  *(100)  at  times  in  the  elliptic  case  by  converting  the  ellipses 
to  circles  of  the  same  area  and  converting  the  ballistic  errors  to  the  circular 
case  as  above.'  Likewise,  *(100)  may  be  used  in  place  of  *(107)  under  the  same 
restrictions  by  setting  the  inner  radius  to  zero.  The  coverage  *(111)  for  a  ripple 
of  N  elliptic  patterns  against  an  area  of  fragment-sensitive  targets,  considering 
edge  effects,  is  computed  according  to  Eq.  (4.82),  if  the  flag  7V1 11*1. 

For  the  special  case  in  which  the  ellipses  are  circles  and  7?l-7?2,  the  coverage 
*(101)  for  r  ripple  of  N  circular  annulus  patterns  against  an  area  of  fragment- 
sensitive  targets  is  computed  from  Eq.  (4.84)  if  the  flag  77101*1.  We  can  use  *(101) 
under  these  conditions  in  place  of  *(108)  by  setting  the  inner  radius  equal  to  zero. 

The  coverage  *(112)  for  a  ripple  of  N  elliptic  annulus  patterns  against  an 
area  of  impact-sensitive  targets,  ignoring  edge  effects,  is  computed  according  to 
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Eq.  (4.85)  if  the  flag  $112=1.  In  general,  Af ( 1 i 2)  should  not  be  used  if  the  ratios 
ETl/Bl,  ET2/B2  are  less  than  5.  If  the  ratios  Vl/Bl  or  W2/B2  are  less  than  5,  it 
is  better  to  set  f/1  and  W 2  to  zero  and  use  X(109). 

For  the  special  case  in  which  the  ellipses  are  circles,  i.e.,  ETlmET2  and 
Wl=W2,  and  the  ballistic  error  is  circular,  i.e.,  !‘?3=54,  a  special  approximation 
may  be  used  if  the  ratios  ETI/S2  and  W2/S3  are  larger  than  3.  In  this  case,  the 
coverage  X(102)  for  a  ripple  of  $  circular  annulus  patterns  against  an  area  of 
impact-sensitive  targets,  ignoring  edge  effects,  is  computed  according  to  Eq.  (4.87) 
if  the  flag  $102=1.  The  same  restrictions  apply  as  for  X(112)  for  the  ratios 
ETl/Bl,  ET2/B2,  Vl/Bl,  and  W2/B2.  Since  the  computation  time  for  X(102)  is  much 
smaller  than  for  K(112),  it  is  sometimes  worthwhile  to  use  A((102)  even  in  the  el¬ 
liptic  case  by  converting  the  ellipses  to  circles  of  the  same  area,  and  by  con¬ 
verting  the  ballistic  error  to  a  circular  case,  letting  S=S3=S4  be  the  equivalent 
ballistic  standard  deviation  in  both  directions.  Likewise,  we  can  use  K(1C21 
in  place  of  X(109)  under  the  same  restrictions  by  setting  the  inner  radius  to 
zero . 


5.5.  SUMMARY 

The  following  tabular  listings  summarize  the  problems  that  are  coded  and  pro¬ 
vide  a  reference  for  further  information. 


RIPPLE  OF  BOMBS,  SEC.  5.2 

Output  Description 

X(U)  Fragment-sensitive  target 

X(l)  Fragment-sensitive  target;  limited  to  the 

case  of  $2 12 

X(15)  Impact- sensitive  target 

X(5)  Impact-sensitive  target;  limited  to  the 

case  of  $212 


Reference 

Eq.  (3.7) 

Eq.  (3.13) 
Eq.  (3.20) 

Eq.  (3.21) 


NOTE:  The  computations  for  X(l)  involve  no  numerical  integrations, 
merely  a  finite  sum.  It  is  very  fast  and  more  accurate  than  *(11)  if 
$  is  small. 


-65- 


Output 

*(12) 

*(2> 

*(13) 

*(3) 

*(14) 

*(4) 


Output 

*(103) 

*(104) 

*(105) 

*(106) 


Output 

*(107) 

*(108) 

*(109) 


Output 

*(110) 

*(111) 

*(112) 


Output 

*(100) 

*(101) 

*002) 


RIPPLE  OF  FIXED  DISPENSERS.  SEC.  5.3 
Description 

Fragment-sensitive  target 
Frogmen l-sens i t ive  target;  limited  to  the 
case  of  17*1  and  NB£ 30 
Approximation  for  *(12);  faster  in  compu¬ 
tation  time,  but  net  as  accurate 
Approximation  for  *(12);  limited  to  the 
case  of  17*1  and  1VSS30 
Impact-sensitive  target 
Impact-sensitive  target;  limited  to  the 
case  of  Nm 1  and  NB£ 30 


RIPPLE  OF  RECTANGULAR  PATTERNS.  SEC.  5.4.1 
Description 

Fragment-sensitive  target,  no  edge  effects 
Fragment-sensitive  target,  with  edge  effects 
Impact-sensitive  target,  no  edge  effects 
Impact-sen3itive  target,  with  edge  effects 


RIPPLE  OF  ELLIPTIC  PATTERNS,  SEC.  5.4.2 
Description 

Fragment-sensitive  target,  no  edge  effects 
Fragment-sensitive  target,  with  edge  effects 
Impact-sensitive  target,  no  edge  effects 


RIPPLE  OF  ELLIPTIC  ANNULUS  PATTERNS,  SEC.  5. 
Description 

Fragment-sensitive  target,  no  edge  effects 
Fragment-sensitive  target,  with  edge  effects 
Impact-sensitive  target,  no  edge  effects 


RIPPLE  OF  CIRCULAR  ANNULUS  PATTERNS,  SEC.  5 
Description 

Fragment-sensitive  target,  no  edge  effects 
Fragment-sensitive  target,  with  edge  effects 
Impact-sensitive  target,  no  edge  effects 


Reference 

Eq.  (4.9) 

Eq.  (4.11) 

Eq.  (4.18) 

Lq.  (4.19) 
Eq.  (4.13) 

Eq.  (4.14) 


Reference 

Eq.  (4.28) 
Eq.  (4.35) 
Eq.  (4.46) 
Eq.  (4.49) 


Reference 

Eq.  (4,59) 
Eq.  (4.65) 
Eq.  (4.70) 


Reference 

Eq.  (4.79) 
Eq.  (4.821 
Eq.  (4.85) 


4.3 


Reference 

Eq.  (4.81) 
Eq.  (4.84) 
Eq.  (4.86) 
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5.6,  QDHP  PROGRAM  INPUTS  . 

The  entire  data  deck  for  the  QDHP  program  is  read  as  an  array.  Space  has 

been  provided  for  up  to  i200  entries.  A  brief  description  of  the  entries  follows. 

Data  (l)«P/?0«number  of  problems  to  be  processed. 

Data  (2)»0iV»number  of  weapons. 

Data  (3)»CWB«number  of  bomblets  per  dispenser. 

Data  (4)«05K>number  of  intervals  in  the  ballistic  pattern  density  distribution 
table  plus  1. 

Data  (5)«0£4.f*>l  if  the  aiming  pattern  offsets  in  range  are  uniform; 

-0  if  the  aiming  pattern  offsets  are  not  uniform  and  must  be  read  in. 

Data  (6)~CXBJ"1  if  the  aiming  pattern  offsets  in  deflection  are  uniform; 

■0  otherwise.  Must  be  read  in. 

Data  (7)*CNPTml  if  there  is  a  ballistic  distribution  table  to  be  read  in: 

-0  otherwise. 

Data  (8)«CW,rM  if  the  distribution  of  the  density  values  in  the  ballistic  toble 
is  cumulative; 

>0  otherwise. 

Data  (9) -PIMP* the  number  of  weapons  per  DPI  (desired  point  of  impact).  This  al¬ 
lows  more  than  one  weapon  to  be  assigned  per  DPI;  i.e.,  CN/  (FTMP"AJ<!) 
must  be  an  integer  i  1.  FIMP*AJJ  must  equal  CN  if  the  weapons  are 
being  dropped  one  at  a  time. 

Data  (lO)»j4AW0*=a  special  flag  to  skip  parts  of  the  program.  Set  *  0  to  compute 
any  £  between  £(1)  and  £(5).  Set  «  1  to  compute  any  other  X  and 
set  »  2  to  compute  a  combination  of  £(1)  through  £(5)  and  any  other 
X  values. 

Data  (ll)-Ml-l.  Compute  £(1).  Otherwise,  set  to  0.  This  definition  of  1  th3t 
means  "compute"  and  0  that  means  "do  not  compute"  is  followed  for 
all  flags  on  the  X  values.  The  AN  number  corresponds  to  the  X  num¬ 
ber  to  be  computed.  £(1)  is  used  for  fragment-sensitive  targets. 
Restriction  CN& 12.  For  larger  CN,  use  £(11). 

Data  (12)»A/Y2“1 .  Compute  £(2).  Case  of  a  long,  narrow,  fixed,  dispenser-type 
pafern.  Restrictions  (CN*l,  CNBZ 30).  Otherwise,  use  £(12). 

Data  (13)-M3*1.  mpute  X(3).  An  expansion  case  of  £(13).  (CW«  1,  CNSs 30). 

Data  (14)'»/UV4»1.  Cou._  *e  £(4).  Used  for  a  dispenser  pattern  against  a  hard  tar¬ 
get.  (CN<-  CNBZ30).  Otherwise,  use  £(14). 

Data  (15)»A/l/5=l.  Compute  a,''.  Impact-sensitive  target.  (CNs 4).  Otherwise, 
use  £(15). 

Data  (16)»/J510»l .  Compute  £(10,  Hand  method.  Rarely  used. 

Data  (17)*/.Afll*l.  Compute  £(11).  -d  for  fragment-sensitive  target  with  6W£l3. 

Data  (18)»AW12“1 ,  Compute  £(12).  Us*.  eor  a  fixed  dispenser  weapon  against  a 

fragment-sensitive  target,  »ds  R\,  R2,  Cl,  1,  00,  51,  and  52  inputs. 

Data  (19)»M13*1.  Compute  £(13).  An  approx,  ‘ion  to  £(12).  One  assumes  a  uni¬ 
form  distribution  of  the  bomblets  range  over  the  distance  2*CL2 
and  a  gaussian  distribution  in  deflei. 

Data  (20)»A514=*1.  Compute  £(14).  Similar  to  £(12)  for  .  ’xed  dispenser  pattern 
against  a  hard  target. 

Data  (21) 5=  1 .  Compute  £(15).  Used  for  impact-sensitive  tai,,  Needs  53 
and  54  inputs,  the  target  dimensions.  (6W>4). 

Data  (22)’*A100»1.  Compute  £(100).  Case  of  a  ripple  of  dispenser?  with  ..  •<*!>- 

nut  or  circular  pattern.  Used  for  fragment-sensirive  target 
no  edge  effects. 

Data  (23) 10 1— l .  Compute  £(101).  Fragment-sensitive  target  with  edge  effects. 

Data  (24)«A102«1 .  Compute  £(102).  Vulnerable  area,  no  edge  effects 

Data  (25)»A103'>1 .  Compute  £(103).  Series  of  rec  jangles,  fragment-.-,  isitive  target, 
but  no  edge  effects. 
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Data  (26)-A104-l.  Compute  *(104).  S3  and  S4  must  equal  0.  Series  of  rectangles, 
fragment- sensitive  target  with  edge  effects. 

Data  (27)-A105-l.  Compute  *(105).  Series  of  rectangles,  vulnerable  area,  but  no 
edge  effects. 

Data  (28)-j4106«1.  Compute  *(106).  Series  of  rectangles,  vulnerable  area  with  edge 
effects.  S3  and  S4  set  to  0. 

Data  (29)-A107-l.  Compute  *(107).  Must  have  values  for  S3  and  S4.  Series  of 
ellipses,  fragment-sensitive  target,  no  edge  effects. 

Data  (30W108-1.  Compute  *(108).  S3  and  S4  equal  0.  Series  of  ellipses, 
fragment-sensitive  target  with  edge  effects. 

Data  (31)-A109-1.  Compute  *(109).  S3  and  S 4  must  be  >  0.  Series  of  ellipses, 
vulnerable  area,  no  edge  effects. 

Data  (32)-A110-l.  Compute  *(110).  S3  and  S4  must  be  >  0.  Elliptical  annulus, 
fragment-sensitive  target,  no  edge  effects 
Dcta  (33)-i41 I 1-1 .  Compute  *(111).  S3  and  S4  must  be  0.  Elliptical  annulus, 
fragment-sensitive  target  with  edge  effects. 

Data  (34)-All2-l.  Compute  *(112).  S3  and  S4  must  be  >  0.  Elliptical  annulus, 
vulnerable  area,  no  edge  effects. 

Data  (35)-0-desired  spacing  in  range  in  feet  between  bomba. 

Data  (36)-0F-deeired  spacing  in  deflection  in  feet  between  bombs. 

Data  (37)-SP-the  product  of  the  reliability  factor  and  the  probability  of  a  kill 
if  hit,  providing  it  is  known.  Otherwise,  1. 

Data  (38)-i4«7J-number  of  wing  stations  i.  1. 

Data  (39)-A3-tsrget  area  dimension  (length).  If  you  are  considering  a  point  target, 
A3-A4-0. 

Data  (40) »X4“ target  area  dimension  (width) . 

Data  (41)-fl3-t«rget  dimension  for  a  hard  target.  Length  z  1. 

Data  (42)-fl4-target  dimension  for  a  nard  target.  Width  z  1. 

Data  (43)-Sl-bomb  ballistic  standard  error  in  range.  (A  fixed  dispenser  is  con¬ 
sidered  a  bomb.) 

Data  (44)-S2-bomb  ballistic  standard  error  in  deflection. 

Data  (45)»S3"dispenser  ballistic  standard  error  in  range. 

Data  (46)"S4-diepenser  ballistic  standard  error  in  deflection. 

Data  (47)-ri«aim  standard  error  in  range.  ox-RF.P1.4828. 

Data  (48)-r2-aim  standard  error  in  deflection. 

Data  (49)»Si/l-coordinate  of  the  aim  point  in  the  x  direction  for  an  offset  target. 
Data  (50)-S£/2-aame  as  Data  (49)  but  in  the  y  direction. 

Data  (51)-£/l-the  number  of  steps  used  in  integration  over  a  quarter  of  the  total 

space.  In  general,  01-8.  If  CL1  is  too  small,  an  increase  in  U 1  will 
solve  the  problem. 

Data  (5?.)-i/2-t/l. 

Data  (53)-CDl-one-hslf  the  length  of  the  dispenser  ballistic  table. 

Data  (54)»40"standard  error  of  the  tabular  values  in  range  of  the  dispenser  bal¬ 
listic  table. 

Data  (55)-ffl-effectlveneaa  radii  in  range  in  the  case  of  an  MAE  type  of  effective¬ 
ness  index,  i.e.,  MAE"'  *fll*/?2.  Tflil. 

Data  (56) "^"effectiveness  radii  in  deflection.  R2zi, 

Data  (57)-01-initial  value  constane  when  a  Carlton-type  damage  function  is  used. 

In  general,  01-1. 

Date  (58) -02-01. 

Data  (59)«#-reliability  factor  for  a  dispenser,  if  known.  Otherwise,  1. 

Data  (60)«W»inner  semiaxis  of  a  dispenser  elliptical  annulus  in  the  x  dlrrctlon. 
Data  (61)-i72-inner  semiaxis  of  a  dispenser  elliptical  annulus  in  the  y  direction. 
Data  (62)-£71"outer  semiaxis  of  a  llspenser  elliptical  pattern  in  the  x  direction. 
Date  (63)-£3T*outer  semiaxis  of  a  dispenser  elliptical  pattern  in  the  y  direction. 
Data  (64)-GC0i*half  the  dispenser  rectangular  pattern  dimension  in  the  x  direction. 
Data  (65)-G6’M-h*lf  the  dispenser  rectangular  pattern  dimension  in  the  y  direction. 
,Dsta  (66)-005-numbsr  of  integration  eMps  f or  offset  ellipse  function.  Uee  40. 


Data  (67)»r,-r  .-arcing  point  for  integration  for  offset  ellipse  function.  Use  0. 
Data  (68)-P.<vyw*l  if  dispenser  deliveries  against  rectangular  target  areas  are  to  be 
computed.  ,Y(100)  through  X(112).  Otherwise,  0. 

Data  (101)=/1»[/1(J)  ,J*l,iV]  array.  Table  of  aim  points  if  spacing  (D)  is  not  uni¬ 
form  in  range,  CXAJm0. 

Data  (301)-B«*[B(tO  ,J»1,/V]  array.  Table  of  aim  points  if  spacing  ( DF )  is  not  uni¬ 
form  in  dei lection.  CXBJ^O. 

Data  (501)=SKT=Sl/J(J) ,J~1,NVI.  Intervals  on  the  ballistic  table.  First  value 
must  be  0.  NPT* 1. 

Data  (601)  =BPJ»BPJ (v)  ,t7*l,JVW.  Density  values  on  the  ballistic  table.  CNCP* 0. 

Data  {701)”T1"T1 (J) ,J“1,NVI .  Cumulative  density  values  on  the  ballistic  table. 
CVCP* 1. 


5.6,1,  Input  Method  Used  in  QDHP  Program 

The  first  card  is  the  number  of  the  first  problem  to  be  solved  and  must 
be  an  integer.  It  is  read  in  on  15  format.  If  the  user  is  cut  off  because  of 
interval  timer  overflow  before  he  has  processed  all  hie  cases,  the  process  allows 
minimum  reabsemblage  of  the  data  deck  without  interrupting  the  run  sequence.  The 
entire  data  array  is  read  in  or.  format  (5I1,I7,5P12.11) . 

The  first  five  entries  are  flags  to  read  or  skip  a  data  field  on  that  card. 

A  blank  or  zero  means  rjad  and  store,  and  a  one  means  skip  the  field.  This  system 
of  reading  data  allows  the  user  to  read  in  only  the  inputs  thac  he  wishes  changed 
for  subsequent  runs.  In  the  second  field  of  each  data  card  is  the  index  of  the 
first  data  item  on  that  card.  The  data  are  read  five  items  per  card.  For  ease 
of  keypunching,  left  adjust  all  data  items.  A  minus  sign  in  column  8  of  the  last 
data  card  of  the  data  deck  is  essential  to  mark  the  end  of  the  set.  To  run  ad¬ 
ditional  cases,  add  the  data  cards  with  the  chr.nged  parameters  and  the  appropriate 
flags  in  columns  1  through  5.  The  last  changed  card  will  have  a  minus  sign  in 
column  8. 

The  listing  on  the  opposite  page  presents  typical  data  entries  for  three 
cases;  the  computer  program  for  the  simplified  weapons  evaluation  model  follows: 
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5.7.  COMPUTER  PROGRAM  FOR  SIMPLIFIED  WEAPONS  EVALUATION  MODEL  (For  System  360/65) 


C  UPDATED  MASTER  DECK  QDHP  PROGRAM  RM-5677  AS  LOADED  ON  DISC  AT  RAND 
DIMENSION  A(200),B(200),A22(200),B22(200),D*TA( 1200) 

DIMENSION  SVI(100),8PI(l00),Tm00) 

COMMON  DATA 

COMMON  SE4,SE3,U3,U4,Z,SQ1,SQ2«Q1«Q2,CL2,PI,SP1 , SQ3,C ,CL3, ATI , AT2 
COMMON  8S1 • 8S2 ,SE12, SE22 ,S5,C45, SW3, SW4, ST1, SU3,SU4, AL 1, AL2 , SHI 
COMMON  FCL3,FCL4,E SI ,ES2 «ER1 »ER2  »ETB1 » ETB2, KR1 » WR2,WB1 » WB2, WS1 »WS2 
COMMON  A22,B22,A1, A2,B1,B2 

COMMON  SQR2,SPI2,PI2,AJKX,BJMX«TU1,TU2,T12»T22 
COMMON  R32,RD1 »RD2»SPR»SW12»SW22»SRR 
COMMON  SN100,SN10l,SNI02,SNil0ySN112,Xll2 

COMMON  X100,Xl01,X102,Xi03,X104»X105,X106,X107,X108,Xl09,XUO,Xlil 
COMMON  N100,N101,N102,Nil0,N112, INUM 
COMMON  KK , KOUT , KOO , KO , K 1 1 ,  LU3 , LU4 , N • NB , LU5 

EQUIVALENCE  ( DATA( 1 ) , PRO) , (DATA( 2) ,CN) , (DATA! 3) ,CNB) , (DATAI4) ,CNVI 

1  ) ,  (DATA(i)  ,CXAJ) ,  (DATA(6)  .CXfjJ  i,  (DAT A ( 7 ) ,CNPT ) ,  ( DATA( 8 ) ,CNCP ) , 

2  (DATA(9),FIMP),(0ATA(10) , ANNO) , (DATA! 1 1 ) , AN1 ) , (DAT A(12 ) » AN2 ) , 

3  ( CAT A( 13), AN 3), (DATA! 14) , ANA) « ( DATA( 15) ,  AN5) ,(DATA(16),ANI0), 

4  i DATA ( 17), AN 11), (DATA( 18 ) , AN12  )  r (DATA( 19), AN 13), ( DA TA( 20), AN 14), 

5  ( OATA (21 ) , AN1 5) , ( DATA (22 ) , A100 ) , (DATA 1 23) , AIOI ) , (DATA(24> , AI02 ) , 

6  ( DATA( 25) « A103) ,(DATA(26),A104),(DATA(27),A105), (DATA (28) , A106) , 

7  ( DATA (29) , A107) , ( OATA < 30) ,A108) , (DATA( 31) ,A109) , (DATA ( 32) ,AI 10) 
EQUIVALENCE  ( DATA (33) , All 1 ) , (DATA(34 )» A1 12 ) , (DATA! 3$ ) ,D) , (DATA! 36) 

1  ,  DF ) ,  ( DAT  A(  37 ) » SP ) ,  ( OAT  A  ( 38)  ,AJJ) ,  (DATA09)  ,A3 ) ,  (DATA(40)  ,A4)  , 

2  ( DATA (41 ) ,83 ) ,( OATA( 42) ,84) , (DATA! 43),SI),(DATA(44),S2) »  OATA  (45) 

3  , S3 ) , (DAT A (46 ) ,S4 ) , (OATA (47 ) , Tl ) , (DATA (48 ) »T2) » (OATA( 49 ) ,SU1 ) , 

4  ( DATA ( 50) , SU2 )»(DATA(51)»U1)»( DAT A ( 32 )cU2),(DATA(53),CLl), 

5  ( DATA( 54) »Q0 ) » ( DATA( 55),Rl),(DATA(56),R2) , (DATA! 57) , D 1 ) , (DATA (58) 

6  , 02 ) , ( DATA ( 59 ) ,R ) , ( DATA (60 ) , W1 ) , (DATA ( 61 ) ,H2) , I  DATA ( 62 ) , ET 1 ) , 

7  ( OATA ( 63 ) . E  T2 ) , ( DATA ( 64 ) ,GCL3 ) , ( DATA ( 65 ) , GCL4 ) 

EQUIVALENCE  (DATA(66)«UU5),(0ATA(6?)«VV5) ,( DATA (68) ,ANUM) 
EQUIVALENCE  ( DATA ( 101 ) «A)« ( OATA ( 301 ) ,8) , ( DATA ( 501 ) *SVI ) , (DATA (601) 

1  ,BPI)  ,  (DATA  ( 701 )  <,TI ) 

2  FORMAT  (8F10.4) 

3  FORMAT  ( 1H0*6X» 1H09X, 2HDF, 8X»  2HSP»9X» IHN,9X»  2HNB , 7X*  3HNVI ,7X,3HNN4 
1/3F10»2,4I10) 

4  FORMAT  ( 1H03X, 13HALL  All)  *  0.) 

5  FORMAT  ( 1H03X, 13HALL  6(J)  *  0.) 

6  FORMAT (6E20. 8) 

7  FORMAT (2A4, IX, I 1,615) 

8  FORMAT  (1H1, IX, 2110, 615) 

9  FORMAT  ( 1H05X, 3HCL1, 7X, 2HQ0, 8X«  2HU3, 8X«  2HU4«8X, 3HSE3o?X, 3HSE4,9X, 1 
1HE,9X,1HC/1H  2X,8F10.5) 

10  FORMAT  (lH0,4X,2HA3,eX,2HA4,8X,2H83,8X,2HB4,8X,2HSl,8X,2HS2,8X,?HS 
13,8X,2HS4/8F10.4) 

11  FORMAT  ( 1H0,4X,2HT1 ,8X,2HT2, 7X, 3HSU1,7X, 3HSU2, 8X, 2HUI, 8X,2HU2, 8X,2 
1HLI»8X,2HQ0/8F10.4) 

12  FORMAT  ( IH03X,4HA(J ) ) 

13  FORMAT  ( 1H03X,4HB( J ) ) 

14  FORMAT  (1615) 
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15  FORMAT  (  IH05X,  4HX<  1  )  ,6X,4HX( 2) ,6X, 4HX(  3)  »6X*4HX(4)  ,6X,4HX!5 ) //4X*5F 
110.6) 

16  FORMAT  ( 1H05X* OF  10.6* 3E20. 8 ) 

17  FORMAT ( 1H05X*  5HX( 10) , 5X, 5HX! 1 1 ) » 5X, 5HX( 1 2 ) ,5X,5HX ( 13 ) , 5X, 5HX !14  )  *  5 
IX,5HXC15),5X,5HX{30)//4X.7F10.6//1H  5X.6HXI  100)  ,4X*6HX!l0l)  *4X,6HX 
2 ( . 02 ) »4X*6HX ( 103 ) *  4X  * 6HX (104)*4Xt6HX(105)*4X*  6HX I  106 )//4X*7F 10.6// 
31H  SX,6HX(  107) ,4X,6HX! 108) ,4X*6HX! 109) ,4X,6HX! 110) ,4X,6HX< 111 ),4X* 
46HX ( i 12  )//4X» 7F10.6) 

19  FORMAT  ( IH0.4X *3HfcT l * 7X, 3HET2 * 6X , 4HGCL 3, 6X,4HGCL4, 7X, 3HUU5, 7X, 3HVV 
15/8F10.4) 

22  FORMAT! 1H02X,2HN0,3X, 2HN l, 3X*2HN2*3X*2HN3,3X,2HN4*3X*2HN5,5H  N10, 

15H  N1 1 f  5H  N12.5H  N13,5H  N14.5H  N15/16I5) 

24  FORMAT  ( IHO, 4X ,?HR1 *8X ,2HR2 * 8X , 2H01 , 8X , 2HD2, 9X, 1HR, 8X» 2HW1, RX, 2HM2 
1/8F10.4) 

25  FORMAT  IlH02X,3H  N0«5H  N100.5H  N101.5H  N102*5H  N1C3.5H  N104,5F  N10 
15»5H  N106 f 5H  N107.5H  N108.5H  N109.5H  N110.5H  NUU5H  N112/1615) 

2006  FORMAT < 1H011IH  EQUATIONS  ARE  INCONSISTENT.  Ctl  MUST  SATISFY  90TH 
ICLl  <  OR  *  TO  C*SQRT ( 3. )  AND  CL l  >  OR  «  TO  1./2./AJJ.  QO-O./fciH  IF 
2THfc  TWO  ABOVE  VALUES  FOR  CL1  ARE  CLOSE*  CHOOSE  ONE  AND  60.) 

2007  FORMAT ( 1H0.61H  CHECK  INPUTS  AJJ  AND  FIMP.  CN/ < AJJ*FIMP )  MUST  BE  AN 
1  INTEGER.) 

732  FORMAT  I1H0.76H  THERE  IS  NO  SOLUTION.  THE  BOMBLET  IMPACT  POINTS 

1  ARE  BUNCHED  MORE  THAN  IN/lH  *66H  A  GAUSSIAN  DISTRIBUTION.  ME  HAVE 

2  ASSUMED  A  GAUSSIAN  DI STRI8UTI0N,/IH  *87h  I.E.  ME  HAVE  SET  CLl-O. 

3  AND  QO-C,  THE  STANDARD  DEVIATION  OF  THE  TABULAR  DISTRIBUTION.) 

74 2  FORMAT  (1H0»70H  THERE  IS  NO  SOLUTION.  PROBABLY  THE  TABULAR  D1STR 
1 1  BUT  ION  IS  8IM0DAL./1H  ,77H  ME  HAVE  USED  A  UNIFORM  DISTRIBUTION 
2 FOR  THE  RANGE  DISPENSER  BOMBLET  IMPACT/IH  ,47H  PATTERN*  I.E.  WE  S 
3ET  00-0.  AND  CL1-C*SQRT ( 3) . ) 

32  FORMAT! 12) 

34  FORMAT  ! 1H0*6X.6HSVI( J) ) 

35  FORMAT  ( 1H0,6X  *5HT  I  (<! ) ) 

36  FORMAT  ( IH0,6X,6HBPI ( J) ) 

39  FORMAT (IH03X*4HK0  -!5«F10.6) 

00  1003  I  »  1*1200 

DATA! I )  *  0. 

1003  CONTINUE 

PI  -  3.14159265 
PI2  «  PI  *  2. 

•*PI  «  SORT(PI) 

SP 12  *  SP1  ♦  .5 
SQR3  «  SORT ! 3.  ) 

SQR2  *  SORT !2. ) 

READ  (5*14)  NUM 
1111  CALL  OECRD  I0AVA) 

40  NPRO  «  PRO  *  .000001 

N  «  CN  ♦  .000001 

NB  *  CNB  ♦  .000001 

NVI  -  CNVI*  .000001 
NXAJ  -  CXAJ*  .000001 
NX8J  *  CXBJ*  .000001 
NPT  -  CNPT*  .000001 
NCPT  «  CNCP*  .00000) 

NIMP  «  FIMP*  .000001 
NN4  n  AJJ  *  .000001 

I NUM  -  ANUM  ♦  .000001 
NO  -  ANNO*  .000001 
N1  «  AN1  *  .000001 

N2  >  AN2  *  .000001 
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N3 

X 

AN3  + 

.000001 

N4 

9 

AN4  + 

.000001 

N5 

X 

r:~ 

.00000! 

N10 

x 

»  •  „  0+ 

.000001 

Nil 

X 

anii+ 

.000001 

N12 

* 

AN12+ 

.000001 

N13 

X 

AN13+ 

.000001 

N14 

X 

AN  14+ 

.000001 

N15 

X 

AN15  + 

.000001 

IF ( INUH.fcQ.O )  GO  TO  55 
NlOO  *  AIOO+  .COOOOl 
NlOl  *  A101  +  .000001 
N102  *  Al’’2  +  .000901 
N103  =  A103  +  .000001 
N104  a  A104+  .000001 
N105  *  A105+  .000001 
N106  «  A106+  .000001 
N107  «  A107+  .000001 
N108  «  A108+  .000001 
N109  *  A109+  .000001 
Nl  10  *  Al 10+  .000001 
Nlll  *  All  1+  .000001 
Nl 12  <  At  12+  .000001 

55  LU5  «  UU5  ♦  .000001 

WRITE  16,3)  NUM,NPRO«NXAJ,NXBJ,NPT »NCPT»NIHP» INUP. 

IF  1 iNUN.GT .0)  GO  TO  57 

WRITE  (6,72)  N0,N1,N2,N3,N4,N5,N10,NII,N12,N13,N14,N15 
GO  TO  56 

57  WRI?E(6,25)  NO.N100,NICI,N102,N103,N104,N105,N106,N107,N1Q8,N109, 
INI 10,NI 11 ,NII2 

56  WRITE16,3)  D,DF, SP,N,N8,NV! ,NN4 
ANNA  a  AJJ 

ANN3  *  CN/ ( ANN4*F IMP) 

IZ  *  1FIX1ANN3) 

AIP  *  IZ 

FP3  »  ANN3  -  AIP 
IF1FP3.EQ.0.)  GO  TO  58 
WRITE  (6,2007) 

GO  TO  392 

58  WRITE (6, 10)  A3,A4,B3,B4,S1,S2,S3,S4 
WRITE  (6, 1 1 ) T1 ,T2, SU1 , SU2,U1 ,U2,CLl , QO 
WRITE  (6,24)  R1,R2,01,02,R,W1,W2 
WRITE  (6,19)  ETl «ET2,GCL3«GCL4,UU5,VV5 
IF  (riPT.EQ.  0)  GO  TO  93 

WRITE! 6, 34 ) 

WRITE  (6,2)  (SVI ( Jl »J«1,NVI ) 

IF (NCPT.EQ.O )  GO  TO  94 
WRITE<6,35) 

WRITE(6,2)(TI( J),J«I,NVI) 

SPIUI  ■  e. 

00  92  R  «  2,NVI 

BP! (K)  «  TI(R)  -  TI(R-l) 

92  CONTINUE 
94  WRITE(6,36) 

WRITE  (6,2)  ( BPI ( J ) , J* 1,NVI ) 

TUI)  «  0. 

00  27  K  >  2.NVI 
TUR)  «  TI(K-i)+BP!(K) 

27  CONTINUE 
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93  SN100  «  0. 

SN101  *  0. 

SN102  *  0. 

SN11G  *  0. 

SNI12  *  0. 

IF (T l.EQ.O. )  Tl«  .00000001 
IF (T2.EQ.0. 1  T2*  .00000001 
A1  *  A3  *  .5 

A2  *  AA  *  .5 

81  *  83  *  .5 

82  *  BA  *  .5 

CNl  »  CN  -  1. 

CN12*  <CN*1.)*.5 
S12  «  Sl*Sl 

$22  *  $2*S2 
TUI  »  T l*Ul 
TU2  *  T2*U2 
ST  1  «  SUl/Tl 
T12  *  Tl*T 1 
T22  *  T2*T2 
ATI  «  Al/Tl 
AT?  «  A2/T2 
ST2  «  SU2/T2 
BT1  »  Bl/Tl 
R12  »  P.l*Rl 
R22  *  R2*R2 
BS1  *  0. 

BS2  «  0. 

IF ( Sl.NE.O. !  BSl  =  Bl/Sl 
IFIS2.NE.0.)  BS2  »  B2/S2 
10112  012  »  01*01 
022  »  02*02 
NKl  «  N  -  l 
ROl  «  R12/(2.*01) 

R02  *  R22/(2.*02) 

All  *  SQRT(ROl) 

AL2  *  SORT (R02  5 
XI  *  0. 

X2  *  0. 

X3  *  0. 

XA  *  0, 

X5  *  0. 

X10  *  0. 

Xli  *  0. 

X12  *  0. 

X13  «  0. 

X1A  »  0. 

X15  «  0. 

X30  «  0. 

X100  *  0. 

X101  =  0. 

X102  *  0. 

XI 03  *  0. 

XlOA  *  0. 

X105  *  0. 

X 106  «  0. 

X10?  »  0. 

X108  >  0. 

X 1 09  *  0. 
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xuc  -  0. 

Kill  «  0. 

XI 12  >  0. 

E  *  0. 

C  *  0. 

K0  *  0 

XPC  «  0 

NRST  «  0 

£0  300  J  >  1,N 

IP  «  1  ♦  tJ~lt/tN!NP*NNA) 

CIP  «  IP 

IFtNXAj.NE.lt  GO  TO  1798 
AIJ)  *  C*tCI«-ANN3*.S-.5t 
SD3  «  0 

1708  A22  ( J )  «  MJimj) 

IFtNXbJ.NE.lt  GO  TO  1799 
AJ  «  J 

JF  «  l AJ- 1. ) /ANNA 
FPC  »  JF 

Bl J)«I (2.*AJ-l.)-ANNA*(2.«FPC«l . ) )*DF*.5 
SOA  •  OF 

1799  622(  J )  *  e<Jl*BU) 

300  CONTINUE 
Ir4NXAJ.EQ.lt  GO  7C  1795 
IFlN.NEoltGO  TO  3C1 

SC3  «  0. 

6C  TC  1795 

301  A JNX  «  All) 

A  JNN  «  All) 

OC  302  J  «  2fN 

4F(A( Jt.GT.AJNX)  GC  TO  303 

IFCAIJ t.LT.AJHN)  A JNN  «  AtJ) 

CO  TO  302 
303  A JNX  «  AtJ) 

302  CONTINUE 

SC3  ■  ( AJNX-AJNN )/CNl 

1795  IFtNXej.EG.lt  GC  TC  199A 
IF (N.NEi 1)  GC  TC  1796 
SOA  *  0. 

GC  TC  199A 

1796  B JNX  «  Bit) 

BJNN  *  Btl) 

CO  30A  J  «  2»N 
IFtBtJt.GT.6JNX)  GC  TC  3P5 
IF t B t J J.lT.BJMN)  BJNN  «  BtJ) 

GC  TO  30A 
305  B JNX  «  BtJ) 

3CA  CONTINUE 

SCA  *  tBJNX-BJNN)/CNl 
199A  IF  t SC3 .EG .0  .  )  GC  TC  1996 
WRIYE(6,12) 

HR ITEt 6t  2 ) t At J),J«1,N) 

1797  IFtSCA.EQ.O.)  GO  TO  1997 
WR!TEt6,13) 

WRITE! 6»  2)tBt J)fJ«l,N) 

GC  TO  1995 
1996  WRITE  16 » A  ) 

GO  TO  1797 
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1997  WRITE  16,5) 

1995  IF1CL1.NE.I.)  GO  TO  312 

700  SZ  *  0. 

IF (NPT .EQ.O)  GO  TO  312 
711  SUSK  *  0. 

SML*  *  0. 

KPO  ■  1 

00  716  LL  *  2»NV  I 

V2  «  ( SVI  (LL)  ♦  SVIILL-in*.5 

SUSM  «  SUSM  ♦  V2*BPI ( LL) 

v3  «  svmu*sviau+sviau*svi(LL-i)*svicLL-i)*sviaL-i) 

SMLM  »  SMLM  ♦  V3/3.*BP  ULL ) 

716  CONTINUE 

C  «  SORT (  SNLH— SUSM* SUSM) 

00  701  K  >  1.NVI 
IF(TI(K).GE.*25)  GO  TO  702 

701  CONTINUE 

702  J  *  K 

SKI  *  SVIIJ-U  ♦  <SVUJ)-SVl(J-m*f  l.25-THJ-l))/ITHJ)-TI(J-I))) 
00  703  K  »  1 , NVl 
IF(TI(K).GE.. 75)  GO  TO  704 

703  CONTINUE 

704  I*K 

sx?  »  sviu-m(svim-sviu-m*<<.75-T!(i-n)/n!<i)-Ti<i-n> > 
SX3  *  (SX2-SX1>*.5 
SX3C  »  SX3/C 

IF1SX3C.GE..6743)  GO  TO  T028 
WR1TE(6,732) 

GO  TO  7085 

7028  1 F ( 3X3C.LE. ( SQR3*.5+.0001 ) )  GO  TO  7030 
WRITE(6,742) 

GO  TO  708 

7030  IF ( ABS ( SX3C-SQR3*.5) .LE. .0001 )  GO  TO  708 
IF ( ABS ( SX3C- • 6744 1 .LE. .0001 )  GO  TO  7085 
£1  «  Q. 

Z2  *  r. 

SZ  *  C*.5 

7040  CSZ  «  SORT ( 3.*IC*C-SZ*SZ 1 ) 

CALL  FF(0.,SX3,CSZfSZ,AF) 

CK  *  AF-.5 

IFCABS1CK1.LE.. 00011  Gn  TO  707 
IF1CK.LT. 0.1  Zl«SZ 
IF (CK.GE.O. >Z2=SZ 
SZ  *  IZl+Z2)*.5 

IF( ABS (Z1-Z2 ) .LE . .01 )  GO  TO  707 
GO  TO  7040 

707  QO  «  SZ 

CL1  *  SORT  1 3.*(C*C-SZ*SZ » ) 

GO  TO  709 

708  QO  *  0. 

CL1  *  SQR3*C 
GO  TO  709 
7085  CLl  »  0. 

QO  -  C 

709  KPO  «  1 

312  QQ2  *  QO  *  QO 

SSR  *  SP  ♦  SPI  *  R1 
SPR  *  SP*PI*R14R2 
SQl  *  SQRT1RDI  ♦  S12) 
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SQ3  *  SORT (CO?  ♦  ROll 
Q1  -  SQR7(SQl*S01  4  T12) 

SQ2  *  SORT (R22/(Z.*D2 )♦  S22) 

Q2  «  SQRT(SQ2*SQ2  ♦  T22) 

1800  SSU3  >  0. 

SSU4  «  0. 

A  JMX  «  ADS  (Sill  -  AID) 

B  JHX  *  ABSISU2  -  8(1)1 
IF(N.EQ.l)  GO  TO  3031 
00  3030  J  *  2yN 
ASUl  «  ABStSUl  -  A( J) ) 

BSU2  ■  ABS( SU2  -  BUI) 

IF ( ASU1 .GT . AJHX)  AJMX  -  ASU1 
I F ( BSU2.GT  cBJHX)  8JHX  »  BSU2 

3030  CONTINUE 

3031  DO  306  J  >  1 « N 
NJi  «  N-J+l 

SSU3  «  SSU3  4  ABS( A(  NJ«|4A(JI) 

S5U4  »  SSU4  4  ABS(B( NJI )4B( J) ) 

306  CONTINUE 

Sl.13  *  SU1  4  SSU3 
SU4  «  SU2  4  SSU4 
KK  «  1 

IF (SU3.EQ.0. )  GO  TO  3165 
U3  «  2.*U1 
GO  TO  3166 

3165  U3  «  U1 

3166  IF ( SU4.EQ.0. )  GO  TO  3167 
U4  «  2.*U2 

GO  TO  317 

3167  U4  »  U2 

317  LU3  «  U3 
LU4  *  U4 

IFUl.EQ.. 00000001)  GO  TO  3192 
SE1  *  (4.4AT1I/U1 
SE12  *  SE1  *  .5 
AK  *  U3 
GO  TO  313 
314  AK  *  AK  -  1* 

313  AX  -  FXJ(AK) 

CALL  H( AX* ATI «HXLI 
AA  *  1.-  HXL 

IF ( AA.GT..0001 )  GO  TO  3175 
GO  TO  314 

3175  SE3  *( SE12  4  AXl/Ui 

318  SE2  *  (4.4AT2J/U2 
SE22  »  SE?  *  .5 
AK  «  U4 

GO  TO  3161 
3182  AK  «  AK  -  1. 

3161  AY  «  FYJ(AK) 

CALL  H( AVy  AT2»HXL I 
AA  «  l.-HXL 

I F ( AA.GT, .0001 5  GO  TO  3185 
GO  TO  3182 

3185  SE4  *( SE22  4  AYI/U2 
3192  IF(NO.EQ.O)  GO  TO  340 
KOUT  »  0 

3002  IF( INUM.GT .0  )  CALL  Z100(N103yNl04,Nl05yN106yNi07yNi08yNl09,Nlll I 
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IF ( Nl 1 . EQ. 0 )  GO  TO  )24 
KO  =  11 

1F(T1. Eg. .OOOOOOOl)  GO  FO  613 
GSH  =  (  AJHX+4.4SQ1 >/TUl 
GS(2  =  (BJHX  <•  4.*S02)/TU2 
S612  =  AMINKSH3,G$E1)*.5 
Sfc22  =  AMINl(SE4,GSE2)*.5 
611  CALL  A50 
6  10  XI  l  *  l 

WRITE  (6,39)  KH.Xlt 
324  IF (N12 .EQ.O)  GO  TO  326 
K0  =  12 

IF(TI. Eg.. OOOOOOOl )  GO  TO  723 
GStl  =  (AJMX  ♦  CL1  ♦  4,*SQ3)/TUl 
GSE?  =  (BJHX  ♦  4.*S02)/TU2 
SE  12  -  AHINH  St3,GSEl  )*.5 
SE22  *  AMINl(S£4»GSE2)*«5 
72.1  C  *  SSR/SQ3 
CALL  A50 
71**  X12  »  l 

WRITE  (6,39)  KO, X12 
125  IFIN13.EU.0)  GO  TO  332 
820  CL 3  *  SURT (CL  1**2  ♦  3.*SU3**2) 

!F(T1. CO. .OOOOOOOl)  GO  TO  825 
GSE2  *  (BJHX  ♦  4.*SQ2)/TU2 
SC22  =  AHIN1(SE4,GS£2)*.5 
825  C  *  (SSR  *  .5J/CL3 
KOO  *  l 
KO  *  13 
CALL  A41 
840  XI)  «  Z 

WRITE  (6,39)  KO, XI 3 
332  KOUT  *  l 

IF (N15.EQ.0)  GO  70  3^4 
KO  »  15 

IFCTl.EO.. OOOOOOOl)  GO  TO  3033 
GSfcl  *  (AJHX  ♦  B1  ♦  4.*SI)/TU1 
GSE2  *  ( BJHX  ♦  B2  ♦  4.*S2)/TU2 
SE 12  *  AHlNl(St3,GSEl  )*.5 
SE22  *  AH)N1(SE4,GSE2)*.5 
3033  CALL  A60 
X15  «  Z 

WRITE  (6,39)  KO, X15 

334  IF (N14.EQ.0)  GO  TO  335 
KO  *  14 

IFIT1.EQ.. OOOOOOOl)  GO  TO  3341 
GSfcl  «  (AJHX  ♦  CL1  ♦  81  ♦  4. *00) /TUI 
GSE2  *  (BJHX  ♦  82  ♦  4.*S2)/TU2 
SE 12  *  AN1NI ( St 3,GSE l )*. 5 
SE22  *  AHIN1 ( S£4,GSE2 )*.5 
3341  CALL  A50 
X14  »  Z 

WRITE  (6,39)  KO, X14 

335  IF (N10.EQ.0)  GO  TO  336 

1515  CL3  “  SQRT(CL1**2  ♦  3.*SQ30*2) 

CL4  »  SQ2*$QR3 

CALL  FF(SU1,CL3,AI,T1,AF> 

ZO  *  AF 

CALL  FF(SU2,CL4,A2,T2.AF) 
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Z  =  ZO  *  AF 

SKC  *  SR*Pl/4.*Rl/CL3*R2/CL< 

153  XIO  *  Z  *  K  *  (l.-(l.-SRC)**NB> 

CALI  FF(SUl,Bl»CL3»Tl»AF) 

ZO  *  AF 

CALL  FF(SU1,82,CL4,T2,AF) 

Z  *  ZO  *  AF 
IFUi.NE.O.)  GO  TO  1560 
X30  «  1.  -  ( l.-SP*Z)**NB 
GU  TO  336 
1560  X30  *  0. 

336  IF(NQ.EQ.l)  GO  TO  392 
340  IF ( SD4#Nt • 0. >  GO  TO  392 

IF(T1.EQ.. 000000011  GO  TO  3400 
GSEl  *  ( AJMX  ♦  BU4.*S1)/TU1 
SE 12  *  AMIN1(SE3,GSE1I*.5 
3400  CALL  XR0(Xl»X2»X3fX4»X5»5D3,SD4) 

WRITE  (6,15)  XI, X2»X3»X4*X5 

392  WRITE  (6,171  X 10, XI 1 , X12 ,X! 3,X14,X15 ,X30,X100,X 101,X102, X103, X104, 
1X105,  X 106,  yi07,X10L,X109,r.klO,Xlll,Xll2 
IF (KPO.EQ.O )  GO  TO  1112 
1113  CLl  »  1. 

QO  *  0. 

11)2  IFJNUM.EQ.NPRO)  GO  TO  2222 
CALL  DECRO  (OATA) 

NUM  «  NUM  ♦  1 

IFISNIOO.NE.O  ,)  N100  *  SNIOO 
IFISN101.NE.0.1  NIG1  *  SN101 
IFISV102.NE.0.)  N102  *  SN102 
IF  ( SN1 10.NC.0.  )  N1 10  *  SNUO 
IF  ( 5N1 12*NE*0*  )  Nt  12  «  SNl.i2 
GO  TO  40 
2222  CALL  EXIT 
ENO 


SUBROUTINE  2 100(N103«N10*,N105,N106,N107,N108«N109,N1 1 1 ) 

DIMENSION  A(200>, Bl 200 ), A22I200), B22( 200), DATA! 1200) 

COMMON  DATA 

COMMON  SE4,SE3,U3,U4,Z,S0l,SQ2*Qlf02,CL2,PI,SPI , SQ3»C,CL3, ATI , AT2 
COMMON  BS1 ,BS2, SE12.SE 22, S 5, C45,SW3, SWA, ST l , SU3,S04«AL1«AL2 , SW1 
COMMON  FCL3,FCL4,ES1,ES2,ER1,ER2,ETB1,FTB2,MR1,WR2,HB1,WB2,WS1,WS2 
COMMON  A22,822,A1,A2,B1*B2 

COMMON  $032, SP I2,P 1 2, AJMX.BJMX, TU1.TU2, T12.T22 
COMMON  R32  »R01 ,RD2 ,SPR,SW12,SW22,SRR 
COMMON  SN100,SN101,SN102,SN110,SN112,X112 

COMMON  X100,X101,X102,X103,X104,X105,X106,X107,X106,X109,X110,X111 
COMMON  N100«N10l,N102,N110,N) 12, INUM 
COMMON  KK , KOUT ,K00,K0,K11, LU  3 , L'J4 , N , NB , LU5 

EQUIVALENCE  ( DATA!  33)  , AU1 ) ,  I OATAI  34 )  ,A1 12)  ,1  DATA! 35)  ,0) ,  (DATA! 36) 

1  , OF ) , (DAT  A(37),SP),( DATA( 38) ,AJ J ), I  DATA! 39) ,A3) , (OATA (405  ,A4) , 

2  (0ATA(4l) ,63) , (DATA(42>  »B4), (DAT A (43), SI ), (0ATA(44),S2),(0ATA(45) 

3  , S3),(0ATA(46),S4),(0ATA(47),Tl),(0ATA(48),T2), (DATA (49 ) , SU1 ) , 

4  ( OATA (50 ) , SU2 ), (OATA( 51), Ui), (DATA1 52) ,U2) , 1 DATA (53) ,CL1 ) , 

5  (0ATA(54) ,Q0) , ( OAT A( 55 ) , R1 ) ♦ ( DATAI56) ,R2 ) , ( 0ATA( 57) ,01 ) , (OATAI 58) 

6  , D2),(0ATA(59),R),(0ATA(60),W1), (DATA (61 ) ,W2 ) , (OAT A (62 ) ,ET1 ) , 

7  (DATA(63),6T2),(DATA(64 ) ,GCL 3 ) ,(0ATA(65),GCL4) 

EQUIVALENCE  ( DATA( 101 ) , A ) , (OATAI 301) »B ) 

2  FORMAT  (8F10.4) 


6  FORMAT (6E2Q.8  ) 

15  f DRMAT  ( IH0,3X,35HGCL3  AND  GCL4  MUST  BF  >  0.  FOR  KO  *14) 

)0  FORMAT  (  lH0jX,3HSWl»8X»2HS5»fiX,?HR3/8F10.4) 

30  f  ORMAN  IH01X.4HK0  *IS,F10.6> 

*i}  FORMAT  ( 1M02X,  3F.HS 3  OR  S4  INRtiTS  ARE  NOT  VALID  FOR  KO  *14) 

44  FURMAN  IH02X.4HK0  =  1 5 ,  2X ,  65HS  3  AND  S4  MUST  BE  0.  FOR  THIS  CASE  TO 
1BF  VALID.  USf  KO- t  CASF.  ) 

pool  FORMAN  IH02X,  B2H  N100» N 101 »N 1 10, M 1 1 2  HAVE  BEEN  SET  TO  ZERO  BECAUSE 
1?M.  MINIET1/R1  .ET2/R2I  WAS  LF  4.1 

2 OOP  FURMAT  ( 1H02X, Z7H  N100 ,N1 0 1 , Nl 02 ,NU 0  HAVE  RFtN  SFT  TO  ZFRU  BECAUS 
U:  THE  M1NIWI/RI.W2/R2)  LE  5.) 

20 0)  FORMAT  I1H02X,S7H  N100,N102  HAVE  BFFN  SFT  TO  ZERO  BECAUSE  SW2/S5 
IIS  5  E  3. ) 

2004  FUKMAI ( IH02X, 105H  N102.Nl 12  HAVE  BEEN  SET  10  ZERO  BECAUSE  EITHER 
l THE  MINIET 1/B1.ET2/B2)  OR  THE  M1NTW1/B1 , W2/B2 )  WAS  LE  5.1 
>005  FORMAT ( 1H02X.78H  RATIO  OF  SW1/S5  LT  3.  BUT  GT.  0.  GIVES  TOO  GREAT 
IAN  ERROR  FDR  X 100 , XlO 1 , X  102. ) 

RJ  =  SORT i R l*R2 I 
R 32  M3  ♦  S3 
S5  =  SOR  T( S  )*S4) 

IF  ( S5.E0.O. )  S4  *  .1 

S52  *  S5  *  S5 

SW2  *  SORT ( E T 1 *E T2 ) 

SW1  =  SORT  1 W1*W2 ) 

SWr.  *  SW2/S5 
WS 14  *  SWl/SO 
If (WS14.EU.0.I  GO  TO  1001 
IF (WS14.GE . 3.  )  GO  TO  1001 
IF (N100.EQ. 1 )  GO  TO  100? 

IF (N101.EU. 1 )  GO  TO  100? 

I F (N10? .NE • 1 )  GO  TO  1001 
1002  SN100  =  NlOO 
SN101  *  Nt  01 
SN 102  *  N 102 
NlOO  *  0 
N101  *  0 
NIO?  =  0 
WRITE  (6,2005) 

1001  SW12*  SW1*SWI 
SW22*  SU2*SW2 
WRITE  (6, 30 ) SW 1 , S5, R 3 
SSR  «  SP*SP I*Rl 
SRR  *  SP4R 1*R2 
S4BB  *  SP44.*B1*B2 
SWSN  *  SW22-SW12 
3002  ER1  «  ET1/R1 
tR2  *  ET2/R2 
ETHl*  ET1/B1 
ETB2  *  ET2/B2 
WRl  >  Wl/Rl 
WR2  »  W2/R2 
WB1  *  Wl/Bl 
WB2  «  W2/B2 
ER  >  AMIN  1 ( ER 1»£R2 ) 

WR  >  AM1NKWR1.WR2) 

EB  *  AMIN1 (ETB1 ,E  TB2 ) 

WB  *  AMIN1 ( W6 l t WB2) 

IF (S3.EQ.0. )  GO  TO  3000 
IFIS4.EQ.0.)  GO  TO  3000 
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WS 1  =  W1/S3 
WS2  «  W2/S6 
ESI  «  ETI/S3 
ES2  *  FT2/S6 
GO  TO  3001 

3000  ESI  *  0. 

ES2  *  0. 

MSI  *  0. 

W$2  «  0. 

3001  IFCER.GT.5.1  GO  TO  3086 
IHN100,EQ.l)  GO  TO  3006 
IFCNlOl.EQ.il  GO  TO  3006 
IF(NUO.NE.l)  GO  TO  3086 

3006  SN100  =  N100 
SN101  *  N101 
SNllC  >  N1 10 
N100  =  0 
N101  =  0 
N110  *  0 
WRITE  (6*2001 ) 

3086  IFCWR.GT.5.1  GO  TO  3087 
IFCWR.EQ.O.)  GO  TO  3087 
IFCSWl.Eti.O.)  GO  TO  3087 
IF(NlOO.EQ.l)  GO  TO  3007 
IFtNlOl.EQ.l)  GO  TO  3007 
IFCN102.E0.il  GO  TO  3007 
IFCNllO.NE.il  GO  TO  3087 

3007  SN100  =  N100 

SN101  »  Nf.Ol 

SN102  -  N.'i02 

SN110  *  NllO 

N100  =  0 
N101  *  0 
N102  »  0 
NllO  *  0 
WRITE  C6»200 21 

i'087  1FCSW5.GT.3. 1  GO  TO  3088 
IFCN100.E0.il  GO  TO  3008 
IFCN102.NE.il  GO  TO  3088 

3008  SN100  -c  N100 

SN102  *  N102 

N100  *  0 
N102  *  0 
WRITE  (6*2003) 

3088  IFCN102.E0.11G0  TO  3090 
IF (N1 12.NE • l 1  GO  TO  3192 

3G?J  IF (EB.LE.5. 1  GO  TO  3089 
IF ( WB.GT.5. 1  GO  TO  3192 
IFCSWl.EQ.O.r  GO  TO  3192 

3089  3N102  *  N10? 

i'N 112  «  N112 

N102  *  0 

Nl  12  «  0 
WHITE  (6*2006) 

3192  IF (N100.E0.0 1  GO  TO  3265 
KO  *  100 

C65  *  (l.-Cl.-SP*R32/SWSW)**NB)  *  R 

5210  IFCSW1.EC.0.)  GO  TO  6531 
SW3  *  SORT (-1 . ♦SW12/S52) 
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GO  TO  4532 

4531  SW3  =  0. 

4532  SW4  =  SURT(-l.*SW22/S52) 

!F(T1.FQ.. 00000001!  GO  TO  4533 
GSEl  «  UJMX  ♦  ET1  ♦  4,*$3)/TUl 
GSF2  =  (BJMX  ♦  ET2  ♦  4,*S4)/TU2 
SE12  »  AMIN1(SE3,GSEI )*,5 

SE22  *  AMINl(SE4,GSE?)*,5 

4533  IFJSWl.EQ.O.)  KO  *  1001 
IF(NIOO.EQ.O)  GO  TO  5211 
CALL  A50 

X100  *  l 

WRITE  (6,39)  KO, X 100 
3265  I F (N102 . tQ.O )  GO  TO  3270 
KO  *  102 

C45  *  R*  ( 1 •- ( 1 .-S4HB/ ( PI*SWSW) ) **NB ) 
IF(NIOO.EO.O)  GO  TO  5210 
5211  IF (SW1.EU.0. )  KO  *  1021 
CALL  A50 
X102  =  / 

WRITE  (6,39)  KO, X102 
3270  IF(NIOI.EO.O)  GO  TO  333 

IF(Tl.tt). ,00000001)  GU  TO  4534 
GSEl  =  (AJMX  *  ETl  ♦  4,*ALl)/TUl 
GSE2  =  (BJMX  ♦  ET2  ♦  4,*AL2)/TU? 

St  12  =  AMINKSC3,GSE1)»,5 
SE22  =  AMIN1(SF4,GSE2)*.5 

4534  KO  =  101 
IF(S3,GI.O.)  GO  TO  4620 
IHS4.Gr.O.)  GO  TO  4620 
C45  =  ( SF*R  32 ) /SWSW 

AL 12  =  AL1*AL? 

I F  ( SW  l  .  F(,  ,0  «  )  GO  TO  4615 
SW3  *  S0RT1SW12/AL12-1.) 

GO  TO  4616 

4615  KO  *  1011 
SW  3  =  0. 

4616  SW4  *  SURT(SW??/AL12-1.) 

CALL  A60 

XIOl  =  7. 

WRITE  (6,39)  K 0 , X 1 0 1 
GO  TO  333 

4620  WRITE  (6,53)  K( 

33  3  FCL  3  =  SORT ( GCL3*GoL  3  ♦  3. ♦KOI) 

FCL4  =  SQRT (GCl4*GCL4  f  3.*R02) 

I F  ( N 10  3.  t'Q.O  )  GO  TO  337 
IF(Tl,fcO.. 00000001)  GO  TO  4536 
GSEl  *  (AJMX  ♦  FCL  3  ♦  4,*S3)/IU1 
GSE2  =  (BJMX  ♦  FCL4  ♦  4.*S4)/TU2 
SE 12  =  AM  INI ( SE  3«GSE 1 )*»5 
SE  2?  =  AMIN1  (St4,GSt'2)*,5 

4535  KO  =•  103 

C45  =  R  +  ( l»-( l.-SPR/l 4.*FCL  3*FCL4))**NB) 
CALL  A60 
X 103  =  l 

WRITE  (6, 39)  KO, X 103 
3  37  IF (N104.bQ«U )  GO  TO  33fl 

IF(T1. tO. .00000001)  GO  70  4536 
GSi  1  -  (AJMX  ♦  GCL3  ♦  4,*ALl)/TUi 
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GSE2  *•  ( BJMX  ♦  GCL4  ♦  4.*AL2)/TU2 
SE 12  *  AMIN1(SC3*GSE1)*«5 
SE22  *  AM1NKSE4,GSE2)*«5 

4536  KO  «  104 

IF ( S3.GT oO. )  GO  TO  3371 
I F ( S4.GT . 0. )  GO  TO  3371 
IF(GCL3.EQ.0.)  GO  TO  3372 
C45  «  SPR/(4.*GCL3*GCL4) 

CALL  A50 
X104  «  l 

WRITE  (6*39)  KO*  X 104 
GO  TO  338 

3371  WRITE  (6*54)  KO 
GO  TO  338 

3372  WRITE(6,15)  KO 

338  FCL3  »  SORT ( GCL3*GCL3  ♦  61*81) 

FCL4  *  SORT ( GCL4*GCL4  ♦  B2*B2) 

I F ( N105 . EO.O )  GO  TO  339 

IF (T l*EQ« « G0000001 )  GO  TO  4537 
GSE1  «  (AJKX  ♦  FCL3  ♦  4.*S3)/TUl 
GSE2  *  (BJMX  ♦  FCL4  ♦  4.*S4)/TU2 
SE12  *  AM1N1(SE3»GSE1)*«5 
SE22  «  AM1N1(SE4»GSE2)*«5 

4537  KO  «  105 

C45  «  R*(l.-(l.-(SP  *81*62 )/ (FCL3*FCL4) ) **NB ) 
CALL  A50 
X105  «  Z 

WRITE  (6*39)  K0.X105 

339  IF(N106.EQ.0)  GO  TO  341 
IF(T1.EQ. .00000001)  GO  TO  4538 
GSE1  «  (AJMX  ♦  GCL3  *  BD/TU1 
GSE2  «  (BJMX  ♦  GCL4  ♦  B21/TU2 
SE 12  *  AMIN1(SE3*GSE1)*.5 
SE22  *  AM  INK  SE4*  GSE2  )*«5 

4538  KO  *  106 

IF ( S3.GT.0. )  GO  TO  3391 
IF1S4.GT.0.)  GO  TO  3391 
IF(GCL3„EQ.0.)  GO  TO  3392 
CALL  A50 
X 106  *  Z 

WRITE  (6,39)  KO* X 106 
GO  TO  341 

3391  WRITE  (6,54)  KO 
GO  TO  341 

3392  WRITE(6»15)  KO 

341  IF1N107.EQ.0)  GO  TO  342 

IF1TI.EQ.. 00000001)  GO  TO  4539 
GSfcl  «  (AJMX  ♦  ET1  ♦  4.*S3)/TU1 
GSE2  *  (BJMX  ♦  ET2  *  4.*S4»/TU2 
SE  12  *  AM!NKSC3,GSE1)*.5 
SE22  «  AM!Nl(SE4,GSE2)*.5 

4539  KO  »  107 

IF ( S3«tQ«0. )  GO  TO  3421 
IF(S4.E0.0.)  GO  TO  3421 
C45  *  R*(l.-(l.-SRR/SW22 )**NB) 

CALL  A50 
X107  *  l 

WRITE  (6,39)  KO* X107 
GO  TO  342 
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3421  WRITE  (6* *>31  K0 

342  IF(N108.fcQ.0)  Gil  TO  343 
1F(T1.EQ..OOOOCOCI)  GO  TO  4540 
GStl  *  ( A JMX  ♦  FT  1  ♦  4.4AL1I/TUI 
GSE2  «  ( BJMX  ♦  ET2  ♦  4„*AL2)/TU2 
SE 12  *  AMIN1(SE3,GSE1 )*.5 

SE22  »  AMIN1(SE4,GSE2)*.5 

4540  KO  -  108 
IFIS3.GT.0.)  GO  TO  3451 
IF ( S4.GT.0. )  GO  TO  3451 
C45  =  SRR/SW22 

CALL  A50 
XI 08  *  Z 

WRITE  (6,39)  KO, X 108 
GO  TO  343 

3451  WRITE  (6,53)KO 

343  IF(N109.EQ.0)  00  TO  3271 
1FIT1.E0.. 00000001)  GO  TO  4541 
GSEl  *  ( AJMX  ♦  ET1  ♦  4.*S3)/TUl 
GSt-2  -  (BJMX  ♦  ET2  ♦  4,*S4)/TU2 
SC  12  =  «MINKSL3,GSE1>*.5 

SE22  =  AM?Nl(SE4,G$E2)*.5 

4541  KO  =  109 

IF( ii.fO.O. )  GO  TO  3431 

IF(S4.Ew.O. )  GO  TO  3431 

C45  =  R*(1.-(1.-S4BB/(SW22*PI))**NR) 

CALL  A50 

X 109  =  l 

WRITE  (6, 39)  KO* X 109 
GO  TO  3271 

3431  WRITE  (6,53)  KO 

3271  IF(NllO.EQ.O)  GO  TO  3273 
IF(T1. EO.. 00000001 )  GO  TO  4542 
GSL1  =  (AJMX  ♦  FT  1  ♦  4.*S3)/TU1 
GSE2  =  (BJMX  ♦  ET2  ♦  4.*S4)/TU2 
SE 12  =  AMINl(S£3,GSfcl )*.5 

SF22  =  AMIN1(SC4,GSE2>*.5 

4542  KO  =  110 

IF (S J.EQ.O. )  GO  TO  3272 
IF ( S4.LQ. 0. )  GO  TO  3272 
C45  =  R*( !•-( I ,-SRR /SKSW ) **NB ) 

CALI  A50 
XUO  =  Z 

WRITE  (6,39)  KU, X  1 10 
GO  TO  3273 

3272  WRITE  (6,53)  KO 

3273  IF(Nlll,fcQ.O)  GO  TO  3275 

IF ( T l.EQ.. 000000/1)  GO  TO  4543 
GSEl  -  (AJMX  «•  ET1  ♦  4.*ALl)/TUl 
GSt?  *  (BJMX  ♦  ET2  ♦  4.*AL2)/TU2 
SE I ?  =  AMINl(SE3,GSFl)*.5 
SE22  =  AMINi(Sl4,GSE2>*.5 
4643  KO  =  111 

IFIS3.GT.0.)  GO  TO  3274 
IF ( S4  »GT .0. )  GO  TO  3274 
C45  *  SRR/SWSW 
CALL  A50 
Xlli  =  7 

WRITE  (6,  )9 )  KO, X 1 1 1 
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GO  TO  3275 

3274  WRITE  (6,53)  KO 

3275  IFINU2.EQ.O)  GO  TO  332 
IFCT1.EQ.. OOOOOOOl!  GO  TO  4544  ' 

GSE1  *  I AJMX  ♦  ET1  ♦  4.*S3)/TU1 
GSE2  «  ( BJMX  ♦  ET2  ♦  4.*S4)/TU2 
SE 12  «  AMiNl(SE3»GSEl)*»5 

SE22  *  AMINl! S£4» GSE2 )♦• 5 
4544  KO  »  112 

IF(S3.E0.0.)  GO  TO  3276 

IFIS4.EQ.0.)  GO  TO  3276 

C45  ■  R4<l.-ll.-S4BB/{Pl*SWSW))**NB) 

CALL  A50 

XI 12  *  Z 

WRITE  (6,39)  KO, XI 12 
GO  TO  332 

3276  WRITE  (6,53)  KO 
332  RETURN 

END 


SUBROUTINE  A50 
REAL  SS(150>/150*0.0/ 

DIMENSION  0ATAU200) 

DIMENSION  A22( 200) «B22(200! 

COMMON  OAT  A 

COMMON  SE4 1 SE3»U3*U4* Z»SQ1 * SQ2, 01 ,Q2 »CL2 »P l » SPl * SC3»C,CL J* AT1.AT2 
COMMON  BSI,BS2,SE12,S£22,S5,C45,SW3,SW4,ST1,SU3,SU4,AL1,AL2,SWI 
COMMON  FCL3,FCL4, ESl«ES2t£Rl,EK2,ET8l*E?B2,WRl,WR2«WBl «W82,W$1 ,WS2 
COMMON  A22»B22fAl»A2,BltB2 

COMMON  SQR2«SPI2tPI29AJMX»BJMXtTUltTU2.T12tT22 
COMMON  R32,RD1 «RD2«  SPR, SW12,  SW22»  SRR 
COMMON  SNld0,SNl01«SN102,SN110«SNll2,XH2 

COMMON  X100tX101tX102«X103fX104fX105vX106»X107tX108«X109vX110tXUl 
COMMON  N100,NIG1,N102,N110«N112« INUM 
COMMON  KK,KOUT,KOO,KO,XlitLU3»LU4,N,N8,LU5 

EQUIVALENCE  ( OAT A(33! , M 11 ) , < 0ATAI34 ) , A1 12 ) , (0ATA( 35) ,D> , (DATA( 36) 

1  ,DF), (DATA137) >SP) , ( DATAi 36) ,AJ j) , <0ATA( 39) , A3 ) » (0ATA(40),A4), 

2  (DATA<41),B3),(DATA(42),B4),(0ATA(43),S1),(0ATA<44),S2)*(0ATA(45) 

3  *  $3 ) * ( DAT  A(46 ) ♦ S4 ) ,  ( DATA (47) »T  1 )  v (DATA(48) » T2) » (DATA! 49) » SUl ) » 

4  ( DATA! 50) « SU2 ) * I0ATAI 51 ) ,Ul ) , ( DATA (52 ) ,U2 ) , ( DATA (53 ) ,CL1 ) , 

5  (DATA(54),Q0),(DATA(55),Rll,(DATA(56)»R2) ,(0ATA(57) ,Dl) ,(0ATA(56l 

6  »D2)»(DATA(59)»R)»(0ATA(6Q)»Wl)»(DATA(61)»W2)»(0ATA(62)fETt)» 

7  (DATA (63) ,ET2 ) , ( DATA (64) ,GCL3 ) , ( DAT A (65 ) * GCL4) 

1  FORMAT  ( 1H02X,3HSE1»7X,3HSE2/2F10.6) 

IFITl.EQ. .00000001)  GO  TO  5057 
505  ZSUM  «  0. 

DO  51  J  *  1,LU3 
AJ  -  J 
SSUM  «  0. 

X  -  FXJ(AJ) 

DO  51?  K  *  l,LU4 
AK  *  K 

Y  «  FYJ(AK) 

CALL  K6(X«Y,PK6) 

CALL  FF(Y*SE22»AT2»l.»AF) 

SSUM  «  SSUM  ♦  PK6  *  AF 
512  CONTINUE 
53  SS(J)  =  SSUM 
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CALL  FFIX,SE12,ATl,l..AF) 

ZSUM  =  7 SUM  ♦  SSI J)*AF 
51  CONTINUE 

5057  IFIT1.GT.. 00000001)  GO  TO  5055 
SE1  *  Al/Ul 
SE12  *  SE1*.5 
S62  *  A2/U2 
SE22  *  SE2*. 5 
ZSUM  *  0. 

DO  5100  J  =  1 > LU3 
AJ  *  J 
SSUM  *  0. 
x  =  rxj(Aj) 

DO  5200  K  *  l.LUA 
AK  *  K 
Y  *  FYJIAK) 

CALL  K6 1 X  ♦  Y » PK6 ) 

SSUM  *  SSUM  f  PK6 
5200  CONTINUE 

SSIJ)  *  SSUM/U5 
ZSUM  *  Z SUM  ♦  SSIJ) 

5100  CONTINUE 

Z  »  ZSUH/U3 
GO  TO  507 
5055  Z  =  ZSUM 

IFISU3.EQ.0.)  Z  »  ?.*Z 
IFISU4.EQ.0.)  GO  TO  50A 
GO  TO  507 

506  l  ~  2.*  I 

507  RETURN 
END 


SUBROUTINE  PP(X,J,PXJ) 

DIMENSION  0ATAI1200) 

DIMENSION  A22I 200) (8221200) 

COMMON  DATA  * 

COMMON  SEAfSE3tU3tUA»Z«S0l,S02»0l;02tCL2tPI»SPIfSQ3»CtCL3fATI»AT2 
COMMON  SSltBS2tSE 12tS£22»$5»C45»SW3*  SWA* ST l * SU3 * 5U4, All ♦ AL2 *  SWl 
COMMON  FCL3, FCL4, ESI , ES2 ,ER 1 , ER2,ETB 1* ETB2 , WR l , WR2 »WBl , WB2 , WS1 , WS2 
COMMON  A22*822«A1»A2«B1»B2 

COMMON  S0R2,SPI2,P!2»AJMX,BJMX,TUI,TU2,T12,T22 
COMMON  R32.ROl»RD2tSPR,SMl2tSVr22,  SRR 
COMMON  SN100*SN101«SN102,SN110iSN112'Xll2 

COMMON  X100»XIOI»X102»X103*XIOA*X105*X106*X107»X108*XI09»X110*X111 
COMMON  N100»N10l»N102«Nll0*N112r INUM 
COMMON  XK»KOUTtKOO,KO,Kll,LU3tLU4»N,NBtLU5 

EQUIVALENCE  I  DAT A 143) , SI ) , I  DATA (A* ) , S2 ) , I  DAT At  55 1 , R 1 ) , (OAT A ( 56 ) t R2 
1) 

PXJ  *  0. 

IF IKOUT.cO.l )  GO  TO  1 
IFIJ.EQ.2)  GO  TO  5 
SQ7  *  7.*SQl 
IFIX.GT.SQ7)  GO  TO  1G0 
PXJ  *  Rl/SQ1*EXP{-(X/S0l)**2*.5)/SQR2 
GO  TO  100 
5  SOB  »  7.*SQ2 

IF  I X.GT.S08)  GO  TO  100 
2  PXJ  ■  R2/SQ2*EXP(~( X/SQ2 )**2*. 5 ) /3QR2 
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GO  TO  100 

1  IF(J.E0.2)  GO  TO  3 
CALL  SF(X,ftl»Sl,SMF) 
PXJ  =  SMF 
GO  TO  100 

3  CALL  SF  ( X,  B2 ,  S2,  SHF ) 
PXJ  *  SMF 
100  RETURN 
END 


SUBROUTINE  XR0(X1 ,X2, X3,X4,X5, SD3,S04) 

DIMENSION  0ATAI1200) 

DIMENSION  A< 200), BI 200 ),A22( 200) ,6221200) 

DIMENSION  1(200), I  STOP! 200 ) « LSTOPI 2001 »KSTQP ( 2001 
COMMON  DATA 

COMMON  SE4fSE3,U3,U4,Z,SQl,SQ2,Ql,Q2,CL2,Pt,SPI ,SQ3,C,CL3,ATl,AT2 
COMMON  BSI,BS2,SE12,SE22,S5,C45,SW3,SW4,ST1,SU3,SU4,AL1,AL2,SWI 
COMMON  FCL3,FCL4,LS1, ES2 ,ER1 »ER2 * ETB1 » ETB2,WRl, WR2, WB1 ,WB2, WS1 , WS2 
COMMON  A22,B22,AI,A2,B1,B2 

COMMON  SQR2,SPI2,Pl2,AJMX,BJMX,TUt,TU2, T12,T22 
COMMON  R32,RDl,RD2,SPR,SW12,SW22,SRR 
COMMON  SN100,SN101,SN102,SN110,SN112,X112 

COMMON  Xi00,X101,Xl02,X103,Xl04,X105,Xl06,Xl07,X108,X109,Xll0»Xlll 
COMMON  N100,N101,N102«N110«N112, INUM 
COMMON  KK,K0UT,K00,K0.KU,LU3,LU4,NtNB,LU5 

EQUIVALENCE  ( DATA! 1 ) »PRO ) * (DATA! ? ) *CN) * ( DATA! 3) ,CN8) ft  DATA! 4) » CNVI 

1  ) ,<DATA(5),CXAJ),(DATA(6),CX8J),(DArA(T),CNPT) , (DATAI 8) ,CNCP) , 

2  ( DATA! 9) »FIMP) » I  DATA  1 10) ,  ANNO ) * (OAT  A 1 1 1 ) * AN1 )» (DATA! 12 ) » AN2) » 

3  (DATA(13)«AN3)«(DATA(14) » AN4) , (DATA! 15) » AN5 ) 

EQUIVALENCE  ( DATA! 33) , A1 II ) , I DATAI 34) , Al 12) , I  DATA! 35) ,0) , IDATAI 36) 

1  »DF ) » (DATAI 37 ) »SP) » (DATA! 38 )  , AJ J ) « (OAT  At  39) , A3 ) , (DATAI40) , A4) , 

2  I  OAT A (41 )  ,B3 ) , I  DATA! 42) ,B4) , I DATA  1 43) »  SI ) » I  OAT  A (44) «32 ) * (DATA (45  t 

3  ,S3),(DATA(46),S4),(DATA(47),Tl),(0ATA(48),T2) , (DATAI49) ,SUl ) , 

4  ( DATA  1 50) » SU2 ) , IDATAI 51 ) ,Ul ) , (DATAI 52) «U2) * ( DATAI 53 )  ,CL 1 ) * 

5  ( DATAI 54), QO) , (DATA( 55) ,Rl > , (DATAI56) ,R2) , I  OAT A(57) ,D1 ) , (DATA(  58 ) 

6  ,02), (DAT A 1 59), R) 

EQUIVALENCE  I DATAI 101) , A ), (DATAI 301 ) ,B) 

N1  «  ANI  ♦  .000001 

N2  =  AN2  ♦  .000001 

N3  *  AN3  ♦  .000001 

N4  *  AN4  ♦  .000001 

N5  *  AN5  ♦  .000001 

AU1  *  ABS(SUl) 

AU2  =  ABSISU2) 

CN1  *  CN  -  1. 

NK1  *  N  -  1 
00  l  J  *  1 ,  N 
I(J)  «  0 
ISTOP(J)  *  o 
LSTOP(J)  »  o 
KSTOP(J)  *  0 

1  CONTINUE 
IF(Nl.NE.O)  GO  TO  2 
IF (N5.EQ.0)  GO  TO  362 

2  DO  390  KK  =  1,N 

I F (N.GT • 12 )  GO  TO  392 
ONEG  *  l. 

LOT  *  0 
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KA  =  KK-1 
FKK*  KK 
WXI  *  0. 

H5  *  0. 

I  F(M0D <KAf ?l )  1  321, 1322,  1321 
VJ21  ONtrG  *  -ONEG 
1322  KOO  «  2 

I F ( T 1 . EQ a  *00000001 )  GO  TO  1324 
GSE2  *  <AU2*B2*4.*S2)/TU2 
SE22  *  AMTNl(SE4»GSE2)*»5 
1324  CALL  A41 
1430  W7  «  l 

344  Q2  *  SORT ( SQ2**2/FKK*  T22 ) 

Ql  *  SORT (SQ1**2/FKK*T12) 

CALL  V ( 2, SU2»  VXJ  ) 

W6  *  VXJ 

350  IF(N5.Nt.O)  GO  TO  510 
LOT  *  1 

352  IF(Nl.EQ.O)  GO  TO  388 
410  IF(SD3.NE.O.)  GO  TO  425 
162  CALL  VI l'SUl'VXJ) 

WXi*  VX J*YF ( KK  ) 

0  =  WX1  *  W6 
GO  TO  385 

510  IFISD3.NE.0*)  GO  TO  520 
1820  Kll  *  3 
CALL  A40 
184  W5  *  l 

0  *  W5*W7 

GO  TO  383 
520  CONTINUE 

425  00  426  J  *  1,KK 
KSTOP(J)  *  N-KK*J 
IFIJ.NE.l)  GO  TO  105 
I(J)  *  J 

LSTOP(J)  *  KSTOPIJI 
GO  TO  106 

105  IIJ)  *  KJ-ll  >  1 
LST0PIJ1  *  tm*NKl-KK*J 

106  I  STOP! J)  *  HINOIKSTOPI J) »LSTOP( J) ) 

426  CONTINUE 
IF(LOT.EQ.l)  GO  TO  429 
IF (N5*NE«0)G0  TO  192 

429  CALL  F (  I , A8AR» ASQySBAR ,RSQ) 

UMA  *  SUi  -  ABAR 
AAM  *  ASQ  -  ABAR**2 
SQU*  2.*SQl*SQt 
CALL  VI 1»UHA, VXJ) 
l  *  VXJ  *  EXPi-FKK  *  AAN/SQll) 
175  WXi  «  WXI*  l  *  SP**KK 
GO  TO  1107 
192  ZSUM  *  0. 

4026  00  193  IJ  *  1,LU3 
A I J  *  IJ 
A I  *  FXJUIJ) 

CALL  PKXI I«AIf PO) 

APO  *  PO 

CALL  FF(AItSH2*ATl»l«,AF) 

ZSUH  «  2SUM  v  APO*AF 


193  CONTINUE 

l  -  ZSUM  *  SP**KK 
IFISUl.NE.O.  )  GO  TO  195 
Z  =  2 .«■  Z 
l‘)5  W5  =  H5  ♦  Z 
1107  IS  =  KK 
GO  TO  101 
104  IS  =  IS  ♦  1 
ISS  =  IS  -  1 
IIIS)  =  KISS)  ♦  1 
IF(IS-KK)  107,108,107 

107  LSTOPl  ISM  )  =  1(1)  ♦  NK1  -  KK  MS  M 
ISTOPI  ISM )  *  MINO(KSTOP(  ISM),LSTOP(ISM>) 

108  GO  TO  102 

101  IIIS)  =  I(IS)  ♦  l 

I F ( (S  -  KK)  109, 102,109 

109  LSTOPl IS+1 )  =  1(1)  ♦  NKl  -  KK  ♦  IS  ♦  1 
ISTOPUSM)  *  MINO(KSTOP(ISM),LSTOP(ISM)) 

102  IFU(IS).GT.ISrnPUS))  GO  TO  103 
IF ( IS.NE.KK )  GO  TO  104 

IF (LOT. EO. 1)  GO  TO  428 
IF (N5.NE .0)GO  TO  427 
428  CALL  F( I , ABAR, ASQ, BBAR,BSQ ) 

UMA  *  SU1  -  ABAR 
AAM  »  A SO  -  ABAR**2 
SQlls  2.*SQ1*SQ1 
CALL  V ( 1 , UMA, VX J ) 

Z  *  VXJ  *  EXP(-FKK*AAH/S011) 

WX1  «  WX1+  Z*SP4*KK 
GO  TO  101 
427  ZSUM  =  0. 

4276  DO  194  IJ  =  1,LU3 
A I J  *  IJ 
A I  *  FXJ(AIJ) 

CALL  PKX( I ,AI ,P0) 

APO  *  PO 

CALL  FF(AI,SE12,AT1,1.,AF) 

ZSUM  *  ZSUM  ♦  APO*AF 

194  CONTINUE 

Z  =  ZSUM  *  SP**KK 
IF ( SU3.NE.0. }  GO  TO  191 
Z  =  2.  *  Z 
191  05  *  W5  ♦  Z 
GO  TO  101 

103  IS  =  IS-1 
IFUS.EQ.O)  GU  TC  382 
IF(I(IS).EO.ISTOP(IS))  GO  TO  103 
GO  TO  101 

382  IF(LOT “1)383,385, 385 

383  X5  *  X5  ♦  0NfcG*W5*W7 
0  =>•  W5*W7 

LOT  «  l 
GO  TO  352 

385  XI  *  XI  ♦  0NEG*KX1*W6 
0  =  HX I *W6 

388  IF (O.LE • .000 1 )  GO  TO  362 
390  CONTINUE 

362  IF(N.GT.l)  GO  TO  39? 
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IF(NB.EQ.O)  GO  TO  392 

KOUT  *  1000 

N4SV  *  0 

N3SV  *  0 

N2SV  *  0 

00  400  KK  3  1,NB 
IF (N4.NE.0 )  GO  TO  404 
IF (N3.NE.0)  GO  TO  404 
IFIN2.FQ.0)  GO  TO  405 
4 r*  FKK  *  KK 

Q2  *  SORT (SQ2**2/FKK+T22 ) 

01  3  SQRT(SQ1**2/FKK*T12) 

KA  *  KK— 1 
ONEG  *  I. 

I F ( NOD ( K A • 2 ) )  401 #402  *401 

401  ONEG  *  -ONEG 

402  IFIB2.EQ.0.)  GO  TO  403 
IFIN4.EQ.0)  GO  TO  403 
KOO  *  2 

IFIT1.E0.. 00000001)  GO  TO  1325 
GSE2  *  ( AU2+B2+4. *S2) /TU2 
SE22  *  AMIN1ISE4,GSE2)*.5 
1325  CALL  A41 
W7  *  l 

403  CALL  V(2,SU2,VXJj 
M6  *  VXJ 

IFIN4.NE.0)  GO  TO  1320 

347  IFIN3.NE.0)  GO  TO  1215 

348  IFIN2.NE.0)  GO  TO  1120 

GO  TO  400 

1320  Kll  *  2 

IFIT1.E0.. 000000011  GO  TO  1323 
GSE1  *  (AUUCLUB1*4.*Q0)/TUI 
SE12  *  AM  I N 1 ( St 3  *  GSE 1 )  * . 5 
1323  CALL  A40 
W4  *  l 

X4  3  X4  *  QNEG*W4*W7 
0  *  W4  *  W7 

IFJ0.GT..000U  f,0  TO  347 
N4SV  3  N4 
N4  *  0 
GO  TO  347 

1215  CL 3  *  SORT (CL  1**2  ♦  3.*  SQ3 **2) 
W3 l* ( SP 12  *Rl/CL3)**KK 
CALL  FF(SUifCL3|AifTl»AF) 

M3  *  W31  *  AF  *  YFIKK) 

X3  *  X3  ♦  ONEG  *  W3*W6 
0  *  M3  *  W6 

IF(O.GT,.OOOl)  GO  TO  348 
N3SV  3  N3 
N3  3  0 
GO  TO  348 

1120  Kll  3  1 

IFIT1.EQ. .00000001)  GO  TO  1121 
GSE  1  »  <AUUCLl*4.*SU3)/TUl 
SE 12  3  AMIN1(S£3»GSE1)*.5 

1121  CALL  A40 
MX2  3  l 

X2  =  X2  ♦  ONfcG  *  WX2  *  W6 


0  *  WX2  *  M6 
IP  ( O.GT. »00G1 )  GO  TO  ‘400 
N2SV  *  N2 
N2  *  0 

400  CONTINUE 

405  IF(N4SV.NE.0i  N4  »  N4SV 
I F ( N3SV.NE .0 )  N3  *  N3SV 
IF (N2SV.NE  .0  )  N2  *  N2SV 

392  RETURN 
END 


SUBROUTINE  A41 
DIMENSION  DATA! 12001 
DIMENSION  A22(200ltB22(200) 

COMMON  DATA 

COMMON  SE4tSE3»U3tU4»Z»SQ!»SQ2»QI»0?»CL2»PI»SPI»S03»CtCL3»ATl»AT2 
COMMON  BSl»BS2.SEl2,SE22,S5,C45,SW3,SW4,STl,SU3,SU4,ALl»AL2,SWl 
COMMON  FCL3,FCL4,ESl,ES2,ERltER2,ETBl»ETB2,WRl,WR2,W81,WB2tWSl,WS2 
COMMON  A22»B22tAl»A2«Bl,B2 

COMMON  SQR2*  SP  I2» P 1 2» A JMX» BJMX»  TUI » TU2 , T 12«  T22 
COMMON  R32,RDI,R02tSPRtSW12tSW22»SRR 
COMMON  SNIOO,$N101,SN102#SN110,SN112,X112 

COMMON  X100.X101, X102fX103«Xi04,X105,Xl06,Xl07,X108,Xl09,XllO,Xlll 
COMMON  N100,N101,N102,N110tN112, INUM 
COMMON  KK  t KOUT , KOO, KO , K 1 1 , LU3 • LU4 ,N , NB , LU5 
EQUIVALENCE  ( DAT A ( 52 ) ♦ U2 ) 

ZSUM  «  0. 

6  DO  4  J  «  1 f  LU4 
AJ  «  J 

VJ  *  FYJUJ3 
CALL  KY(YJtZKY) 

ZK1  «  ZKY 

CALL  FF(YJtSE22,AT2*l.»AF) 

ZSUM  «  ZSUM  ♦  ZK1*AF 
4  CONTINUE 
Z  *  2.*  ZSUM 
IFCSU4.EQ.0. )  GO  TO  100 
Z  =  Z  *  .5 
100  RETURN 
ENO 


SUBROUTINE  VIJ»X,VXJ) 

DIMENSION  DATA! 1200 1 
DIMENSION  A22I 200) *822(200) 

COMMON  DATA 

COMMON  SE4»SE3*U3»U4»Z*SQl*SQ2»Ql*Q2*CL2»PI*SPI*S03»CtCL3*ATl*AT2 
COMMON  BSl.BS2,SEl2»SE22,S5,C45,SH3,SW4,STi,SU3fSU4,ALl,AL2,SWl 
COMMON  FCL3,FCL4,ESl,ES?*ERl,FR2fET81,ETB2,WRl,WR?,WBl,WB2,WSl,WS2 
COMMON  A22*B22,Al»A2,Bl,B2 

COMMON  SQR2,SPI2»P12»AJMX,HJMX, TUUTU2*  T  12*  T22 
COMMON  R32»RDl»R02*SPR*SW12*SW22*SRR 
COMMON  SN100,SN101,SN102,SN110,SNU2*X112 

COMMON  X100,Xl01,X102*Xi03,X104,X105»Xl06,Xi07,X108*X109,Xll0,Xlll 
COMMON  N100,N101,N10?,N110,N112, INUM 
COMMON  KK , KOUT ,K00 , KO  *K 1 1  *  LU3 , LU4 ,N , NO , LU5 
EQUIVALENCE  ( OATA< 55 ) ,R II » 10ATA ( 56) , R2) 

VXJ  =  0. 
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AKK  =  KK 

P2GK  *  SQRT ( 2. *P I /AKK  ) 

IFU.EQ.2)  GO  TO  2 

RP  *(R1/ISQ1*SQR2) >**KK  *  P20K  *  SQ1/Q1 

AXQ  *  X/Ql 

AQ1  *  Al/QL 

CALL  HS(AXQtAOlySH) 

VXJ  «  RP  *  SH 
GO  TO  1 
2  AXQ  *  X/Q2 
AQ2  «  A2/Q2 
CALL  HS( AXQyAQ2*  SH) 

VXJ  *(R2/(SQ2*SQR2) )**KK  *  P20K  *  SQ2/Q2  *  SH 
1  RETURN 

END 


FUNCTION  SMGIXX) 

DIMENSION  DATA! 1200) 

DIMENSION  A22I200) y H22120CI 
COMMON  DATA 

COMMON  SE4,SE3,U3yU4,2,SQl,SQ2,QlyQ2,CL2,PIySPI,SQ3yC,CL3,ATl,AT2 
COMMON  8SI»BS2tSE12»SE22»S5»C43fSW3»SHA»STlf SU3y  SU4y  AL 1? AL2 »  SWl 
COMMON  FCL3#FCL4»ESl»ES2»ERl»ER2»ETBI*FT82»WRltWR2iWBl»WB2»WSl»WS2 
COMMON  A22,B22,Al,A2,BiyR2 

COMMON  SCR2,SPI2yPI2yAJMXyBJMX,TUlyTU2yT12,T22 
COMMON  R32»RDl*RD2tSPRf  SW12ySW22?SRR 
COMMON  SNI00«SN10IfSN102fSN110*SN112»Xll2 

COMMON  X100,XlOl,X102,Xl03,X104yXI05,Xl06,XlC7,X108,XlQ9,XllO,Xlll 
COMMON  N100*N101.N102tNL10»NU2»INUM 
COMMON  KKyKOUT  *KOO»KO»K11 y  LU3 y LU4  y N y NR  y LU5 
SMG  =  EXP(-XX*XX*.5»/SQRT(PI2) 
l  RETURN 
END 


FUNCTION  YF(KKK) 

DIMENSION  DATA! 1200) 

DIMENSION  A22( 200 ) y B22( 200 ) 

COMMON  DATA 

COMMON  St4»SE3yU3yU4»2ySQl,SQ2,Ql,Q2,CL2yPIySPI,SQ3,C,CL3,ATlyAT2 
COMMON  BSl,BS2,SEl2,SE22yS5,C45,SW3ySW4,SU,SU3,SU4,ALlyAL2ySWl 
COMMON  FCL3,FCL4,FSiyES2,ERiyER2»ETBlyETB2,WRl,WR2,WBl,WB2yWSl,WS2 
COMMON  A22yB22,AlyA2yBlyB2 

COMMON  SQR2ySPI2yPI2»AJMX,BJMX,TUL,TU2yTl2,T22 
COMMON  R32yR01yR02ySPRySW12ySW22y  SRR 
COMMON  SNl00ySN101ySN102ySNll0«SNH2yXtl2 

COMMON  XlOOyX10lyXi02yXI03yX104»X105yXl06,Xl07yXl08yX109,XUO,Xin 
COMMON  NI00yNl0IyN102yN110yN112y INUM 
COMMON  KK»KnuTyK00yK0,Kll,LU3,LU4,N,NByLU5 
EQUIVALENCE  I0ATA(37I ,SP) 

ANO  *  N 

IF (KOUT «EQ. 1000)  ANO  =  NB 
YF  =  l. 

00  1  J  =  l,KKK 
AJ  *  J 

YF=YF  ♦  (ANO-AJ+I. )/AJ  *  SP 

1  CONTINUE 

2  RETURN 
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END 


FUNCTION  FXJIXX) 

DIMENSION  DATA! 1200) 

DIMENSION  A22I200)»B22<200) 

COMMON  DATA 

COMMON  S£4*S£3,U3»U4,Z*SQl,SQ2*Ql*Q29CL2,Pl»SPI*SQ3*C,CL3.ATl»AT2 
COMMON  BSl»BS2»Stl2,S£22tS5*C45,SW3»SW4,STI»SU3»SU4*ALl»AL2»SWl 
COMMON  FCL3  *  FC1.4*  ES1*  ES2»ERl*ER2*ETBl*ETB2*WRl» WR2* MBl *WB2  * WS1  *  WS2 
COMMON  A22,B22,Al»A2»Bl,82 

COMMON  SQR2* SPI2*PI2,AJMX*BJMX.TUl*TU2»T12»T22 
COMMON  R32»RDitRD2fSPR«SW12«SW22» SRR 
COMMON  SN100fSNl01,SNl02.SNll0*SNll2*Xll2 

COMMON  X100tXl01,XI02.X103,X104»Xl05.Xl06,Xl07»X108*X109*X110.Xlll 
COMMON  N100*N101*N102*N110*N1L2« INUM 
COMMON  KK*K0UT*K0n»K0*Kll,LU3,LU4,N,NB*LU5 
EQUIVALENCE  I  DATA  151 ) ,U1 ) 

IF (XX.NE.O. )  GO  TO  2 
FXJ  «  0. 

GO  TO  3 

2  IF ISU3.NE .0. )  GO  TO  4 
FXJ  *  (2«*XX-l. )  ♦  SE 12 
GO  TO  3 

4  FXJ  «  (2.*XX-1.-2.*U1 I  *  SE12 

3  RETURN 
END 


FUNCTION  FYJIXX) 

DIMENSION  0ATAI1200) 

DIMENSION  A22I200) t B22I200) 

COMMON  DATA 

COMMON  SE4,St3,U3.U4f Z ,SQl »SQ2,Ql»Q2,CL2t PI  * SPI » SQ3,C,CL3» ATI ♦ AT2 
COMMON  BS l * 8S2 * SE 12  * SE22  * S5 tC45 * SW3* SW4* ST1 * SU3* SU4* AL 1  * AL2 * SWl 
COMMON  FCL3.FCL4fESl,ES2fERl,ER2,ETBl,ETB2.WRl»WR2,WBl,WB2.WSl»WS2 
COMMON  A22*B22,Al,A2f Bi,B2 

COMMON  SQR2»SPI2»P12»AJMX» BJMX* TUI »TU2 *T 12.T22 
COMMON  R32«KD1»RD2*  SPK*  SW12*  SW22*  SRR 
COMMON  SNlOOtSNLOlf SN102* SN1 10* SNl 12* XI 12 

COMMON  X100*X101»X102»Xl03»X104»Xi05*X106*Xi07»  XI08tXt09*Xl 10»Xlli 
COMMON  NIOO*N101*N102*N110*N112* INUM 
COMMON  KK,K0UT,K00»K0.K11,LU3,LU4*N»NB,LU5 
EQUIVALENCE  I DATAI 52 ) *U2 ) 

IFIXX.NE.O.)  GO  TO  2 
FVJ  *  0. 

GO  TO  3 

2  IF ( SU4.NE.0. )  GO  TO  4 
FYJ  *  I2.4XX-I.}  *  SE 22 
GO  TO  3 

4  FYJ  =  I2.*XX-1.-2.*U2)  *  SE22 

3  RETURN 
END 


SUBROUTINE  A40 
DIMENSION  DATAI 1200) 
DIMENSION  A22(200)*0?2<200) 
COMMON  DATA 
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COMMON  SE4,SE3yU3yU4yZySQl,SQ2,QlyQ2yCL2,PI,SPI , SQ3,C,CL3y ATI y AT2 
COMMON  BSlyBS2ySE12ySE22,S5fC45,SW3,SW4,STl,SU3,SU4f AL1,AL2,SWI 
COMMON  FCL3y  FCL4f ESly ES2y ERl yER2yETBl»FTB2»WRltWR2»WBl y  WB2 * WS 1 1 WS2 
COMMON  A22,B22,Al,A2yBi»82 

COMMON  SQR2*SPI2,PI2y AJMX, 8JMXy TU1»TU2. T 12, T22 
COMMON  R32yRDlyRD2ySPK,SW12ySW22ySRR 
COMMON  SN100ySN101ySN102,SNU0ySNll2yXll2 

COMMON  XlOOyX101,X102yX103,X104,Xl05,Xl06yX107,Xl08yX109,X110,Xlll 
COMMON  N100,N101yNl02yN110yN112y 1NUM 
COMMON  KK  » KOUT y  KOOy KO ♦ K 11 , LU3 1 LU4 , N *  NB  ,  LU5 
EQUIVALENCE  IDATA<51),U1> 

ASUM  *  0. 

6  00  A  J  *  If LU3 
AJ  «  J 

XJ  *  FXJ(AJ) 

CALL  KX(XJ.ZKX) 

ZK  I  *  ZKX 

CALL  FF(XJ»SE12t ATI* l.»AF) 

ASUM  *  ASUM  ♦  ZK1*AF 
4  CONTINUE 

Z  *  2.*ASUM*YF(KK) 

IF ( SU3.EQ.0. )  GO  TO  100 
Z  *  Z  *  .5 
100  RETURN 
END 


SUBROUTINE  KX(ZA,ZKX) 

DIMENSION  DATAI12001 
DIMENSION  A22(200),B22(200i 
COMMON  DATA 

COMMON  SE4tSE3fU3fU4,ZfSQlfSQ2tQlfQ2«CL2fPIfSPltSQ3«CfCL3fATl.AT2 
COMMON  BS1 f BS2f SE 12t SE22 1 S5«C4S« SW3t SW4«  ST  1 , SU3 y SU4 y AL 1 , AL2 y SW l 
COMMON  FCL3«FCL4»ESlf ES2f ERlt ER2tETBlt ETBZfMRl* WR2tWBl «WB2  »WSl ,  WS2 
COMMON  A22tD22»Al»A2fBl«B 2 

COMMON  SQR2t SPI2t PI2tAJMX%BJMXt TUl*TU2f T 12,722 
COMMON  R32tR01fRD2f  SPRt  SM12*  SW22f  SRR 
COMMON  SNlOOt SNlOl *SN102fSN110>SNil2»X112 

COMMON  XIOtif X101 f Xl02f X103yX104f X105 tX 106? X107t X 108,X109,X110,Xlll 
COMMON  NlOOy NlOly N102y N110yNll2y INUM 
COMMON  KKfK0UTfK00fK0«KU»LU3tLU4tNfNBfLU5 

EQUIVALENCE  {DATA  143) « SI ) « (DATAI47) yTl ) y I  DATA (49 )  ,SUl ) » (DATA( 53) ,C 
1L1J# (DATA(54),Q0Jy(0ATA{55),Rl),lDATA{56),R2) 

ZKX  «  0. 

IF(Kll.EQ.l)  GO  TO  1 
IFIK11.EQ.2)  GO  TO  2 
IFIK11.EQ.3)  GO  TO  3 

1  XTS  «<ZA  •  Tl  ♦  SU1 J/SQ3 
CLS  «  CL1/SQ3 

CALL  HS(XTSyCLSySH) 

ZKX  *{ SH  •  SP1  4  R1/SQ3)**KK 
GO  TO  100 

2  XT  «  (ZA*mSUl) 

CALL  FFIXTyBlyCLlyQOy AF.) 

ZKX  «  AF  **  KK 
GO  TO  100 

3  XT  «<ZA*Tl*SUn 
CALL  SFIXTyBlySlySMF) 

ZKX  »  SMF**KK 


I 
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100  RETURN 
END 


SUBROUTINE  KY(ZZ,ZKY) 

DIMENSION  DATA (1200 ) 

DIMENSION  A( 200) t B( ?00 ) , A22 ( 200) ,622(200) 

COMMON  OATA 

COMMON  SE4,SE3,U3,U4,Z,SQ1,SQ2,01,Q2,CL2,PI,SPI,SQ3,C.CL3,AT1,AT2 
COMMON  BSI,8S2,S£12,SE22, S5»C45»SW3,SW4,ST 1 ,SU3,SU4, AL 1, AL2, SHI 
COMMON  FCL3»FCL4»ES1»ES2»ER1»ER2*ETB1»ETB2,WR1»WR2,WB1»WB2»W$1»WS2 
COMMON  A22,B22,AI.,A2,B1,B2 

COMMON  SQR2, SP 12 , PI 2 , AJMX, BJMX, TU1»TU2»TI2,T22 
COMMON  R32,RDI,RD2,SPR,SW12, SW22, SRR 
COMMON  SN100,SN101,SN102,$N110,SN112,XU2 

COMMON  X100,X101«X102,X103,X104,X105,X106,X107,X1C8,X109,X110*X111 
COMMON  Nl00,Ni01,N102,NII0,NIt2, I  MUM 
COMMON  KKfKOUT ,K00,K0,Kll,LU3,LU4,N,NB,LU5 

EQUIVALENCE  ( DAT  A (44 ) , S2 ) , (DATA  1 47) , T I) , ( DATA ( 40) , T2) ; (DATA! 49) , S 
1UI),(DATA(50), SU2 ) 

EQUIVALENCE  ( DAT A ( 101 ) , A ) , ( DATA! 301) «» ) 

ZKY  =  0. 

IF(KOO.EQ.l)  GO  TO  1 
IF1K00.EQ.2)  GO  TO  2 

1  ZKK  *  1. 

00  3  J  «  1,N 
AJS  *  SUI  -  A( J) 

CALL  FF(AJS,CL3,A1,T1,AF) 

ZAF  *  AF 

YTB  *  ZZ  *  T2  -  B(J)  ♦  SU2 
CALL  PP( YTB, Z, PXJ ) 

ZKK  *  ZKK  •  (l.  -  ( 1 .-( I .-C*PXJ (•♦NB) ) 

3  CONTINUE 

ZKY  «  (1.  -  ZKK)  *  ZAF 
GO  TO  100 

2  YTS  *  ZZ*T2*SU2 

CALL  SF ( YTS,B2»S2,SMF ) 

ZKY  »  SMF**KK 
100  RETURN 
END 


SUBROUTINE  F( I ,ABAR , ASQ, BBAR, BSQS 
DIMENSION  DATA( 1200) ,! (200) 

DIMENSION  A( 200 i, B( 200 ),A22( 200), 622(200) 

COMMON  DATA 

COMMON  SE4,S£3«U3,U4«Z ,SQ1 ,SQ2,Q1 ,Q2,CL2,PI,SPI,SQ3»C,CL3,AT1,AT2 
COMMON  BS1,BS2, SE 12, SE22, S5»C45, SW3*  SW4»ST i , SU3, SU4, AL1 , AL2 « SMI 
COMMON  FCL3, FCL4, ESI, ES2, ERl,ER2,ETBl,ETB2,WRl,HR2,WBl«WB2,WSl«WS2 
COMMON  A22,B22,A1,A2,B1,62 

COMMON  SQR2, SP 1 2, PI 2, AJMX, BJMX, TU1,TU2,T12,T22 
COMMON  R32,RD1,RD2,S?R*SW12,SM22«SRR 
COMMON  SN100,SN101,SN102,SN110,SN112,X112 

COMMON  Xi00,X101,Xi02,X103,Xi04,X105,Xl06,X107,Xl08,Xi09,Xii0,Xlil 
COMMON  N100,N101,N102,N110,N112, INUM 
COMMON  KK,KOUT *KC0,K0,K11 «LU3»LU4»N*NB,LU5 
EQUIVALENCE  (DATA! 101 )  ,A) , (DATA! 301 ) ,8) 

AK  *  KK 
ABAR  *■  0. 
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ASQ  *  0. 

BBAR  *  0. 

BSO  *  0. 

DO  l  J  *  1 ,  XX 
Ml  «  HJ) 

ABAR  *  A8AR  ♦  A(M1)/AK 
ASQ  *  A$Q  *A22 ( Mi ) /AX 
BBAR  *  BBAR  ♦  B(M1)/AX 
BSO  *  BSO  ♦B22(Ml)/AK 
1  CONTINUE 
RETURN 
END 


SUBROUTINE  PXX( I  *  A I J,PO> 

DIMENSION  DATA! 1200) « 1 1  200) 

DIMENSION  A (200) *  B(200 )*A22(200)rB22(200) 

COMMON  DATA 

COMMON  SE4,SE3,U3,U4,Z,SQ1,SQ2,Q1,Q2,C12,PI ,SPl ,SQ3,C,CL3,ATl ,AT2 
COMMON  BS1*B52*SE12«SE22*S5*C45*SM3«  SW4,  ST  1*  SU3«SU4«AL  1 » AL2  *  SW1 
COMMON  FCL3,FCL4,ES1,E$2,ER1,ER2,ETBI,ETB2,WR1,WR2,WBI,WB2,WSI,WS2 
COMMON  A22,B22,A1,A2,BI,B2 

COMMON  SQR2*  SP 12*  PI2*£JMX*  BJMX*  TU1*TU2*T12«T22 
COMMON  R32,RD1,RD2,SPR,SW12,SW22,SRR 
COMMON  SN100*SN1CI*SN102«SN110*SN1!2*X112 

COMMON  X100,X101 *  X102«X103*  X104*  X105*  X106*  XIQ7*  X108»X109*XI 10*Xl 11 
COMMON  N100«N101*N102*N110*N112* INUM 
COMMON  KKt KOUT *KOO*KO*Kl 1 *LU3*LU4*N*N8*LU5 
EQUIVALENCE  ( OATA( 43) *Sl)*iDATA(47)*Tl), I  DATA 1 49) *  SU1 ) 

EQUIVALENCE  ( DATA( 101 ) , A) , ( OATAl 301 ) ,B) 

PO  *  1. 

na  i  j  *  lyKK 

Ml  «  IU> 

SFX  «  (AIJ*T1*SU1-A(MI>) 

CALL  SF ( SFXtBI y SI t SMF ) 

PO  «  PO  ♦  SMF 
1  CONTINUE 
RETURN 
END 


SUBROUTINE  K6( XX, VY»PK6) 

REAL  XUA(200 )/2Q0*0»0/» YUB ( 200) / 200*0. 0/ 

DIMENSION  OATAUtCO) 

DIMENSION  AI200) *8(200), A22(200) ,B22 (200) 

COMMON  DATA 

COMMON  SE4 * SE3  *U3*U4*  2  * SQ1 * SG2 *  Q 1  *  02  * CL2 *  P I v  SP l *SQ3«C*CL3*AT1*AT2 
COMMON  BS1 ,BS2 *SE 12»SE22  » $5 *C45  * SW3 , SW4, ST  1 , SU3  » SU4, AL lv  AL2* SWl 
COMMON  FCL 3* FCL4,ES1*ES2*ERI,ER2*ETB1»ETB2*WR! * WR2* WBl *WB2  * WS1 *  WS2 
COMMON  A22tB22*Al*A2*Bi»82 

COMMON  SQR2*SPI2*P!2,AJMX*BJMX*TUi«TU2yT 12*T22 
COMMON  R32*R01»R02*  SPR*  SW12*  SK22*  SRR 
COMMON  SN100* SNIOI* SN102*SN1 10*  SN112*X112 

COMMON  X100*X101 ,X102*X103* Xl04*Xl05*XI06*X107*Xl08,X109*Xll0f Xll L 
COMMON  N100*N101 *N102*N110*N112* INUM 
COMMON  KKtKOUT *KOO*KO*Kl 1*LU3*LU4«N«NB*IU5 

EQUIVALENCE  ( DATA (33) , A1 II ) , (DATA! 34), Al 12 ) , ( DATA! 35 ) ,0) ♦ (DATA( 36) 
1  *DF),(DATA(37)*SP)«( DATA ( 33) *AJ J) * (DATA (39) » A3 ) » (DAT A (40 ) * A4 ) , 

'  2  ( DAT A (41 ) , B3) , (DATAI42) *84) , ( DATA! 43) , SI) * (DATA (44) ,S2) * (DATA (45 ) 


3  *S3)*(DATA(46)»S4),(DATA(47)»TI)»(DATA(48),T2)»(DATA(49)»SU1), 

4  (DATA(50),SU2),(0ATA(51),Ul) ,(DATA(52) ,U2 ) , < DAT A (53 ) , CL  1 ) , 

5  (DATA(54),Q0),(DATA(55),R1) , (OATA( 56) ,R2) ,(0ATA(57) ,D1 ) , ( DATA( 58 ) 

6  ,02) , (DATA(59),R) , (DA TA ( 60 ) , W1 ) , (OATA( 6 1 ) , W2 ) , (DATA (62) ,ETl) , 

7  ( DATA ( 63 ) , ET2 ) » ( DAT A ( 64 ) , GCL3) , (OAT A (65) , GCL4) 

EUUI VALENCE  ( DATA ( 101 ) ,A ),( DATA ( 301 > ,B ) 

PK6  *  0. 

PKP  =  l. 

00  500  J  =  1 , N 

IFCT1.EQ. .00000001)  GO  TO  501 
XUA(J)  =  XX  *  T1  ♦  SU1  -  A ( J ) 

YUBCJ)  =  YY  *  T2  ♦  SU2  -  B(J> 

GO  TO  500 

501  XUA(J)  *  XX  ♦  SU1  -  A(J) 

YUB( J)  *  YY  ♦  SU2  -  BIJ) 

500  CONTINUE 

IF (KO.EQ. I 1 )G0  TO  11 
I F (KO. EG* 12 )G0  TO  12 
IF (KO.EQ. 14)G0  TO  14 
IFtKQ.cQ. 15)G0  TO  11 
IF(KO.L:Q.100)  GO  TO  4 
IF (KO.EQ* 1001 )  GO  TO  5 
IF(KO.EO.lOl)  GO  TO  101 
IF(XO.EQ.lOll)  GO  TO  1011 
IF(K0.CQ.102I  GO  TO  4 
IF (KO.EQ* 1021 )  GO  TO  5 
I? IK0.E0.103 )  GO  TO  103 
IF (KO.EQ. 104)  GO  TO  104 
IF (KO.EQ. 105 1  GO  TO  105 
IF (KO.EQ. 106 )  GO  TO  106 
IF (KO.EQ. 107)  GO  TO  107 
IF (KO.EQ. 1081  GO  TO  108 
IF(K0.EQ.109)  GO  TO  109 
IF (KO.EQ. 1 10 )  GO  TO  110 
IF (KO.EQ. Ill)  GO  TO  111 
IF (KO.EQ. 1 12 )  GO  TO  112 

11  00  99  J  «  1 ,  N 
XTA  *  XUA(J) 

CALL  PP( XTA, 1 , PXJ) 

PXl  *  PXJ 
YTA  *  YUB(J) 

CALL  i>P( YTA, 2 »PXJ) 

PY1  *  PXJ 

PKP  »  PKP  *  (l.-PXl*PYl*SP) 

99  CONTINUE 
GO  TO  100 

12  HLS  »  CL1/SQ3 
DO  98  J  «  1,N 
YTS  *  YUB(J) 

CALL  PP(YTS»2»PXJ ) 

PP1  «  PXJ 

XTS  *  XUA(J)/SQ3 

CALL  HS(XTS,HLS« SH) 

PKP  *  PKP*(l.-PPl*C*SH)**NB 
98  CONTINUE 
GO  TO  100 
14  00  97  J  «  1  »N 
XTS  *  XUA(J) 

CALL  FF(XTS?B1,CLI,Q0,AF) 
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7 JX  *  SP  *  AF 
YTS  *  YUB ( J ) 

CALL  PP( Y7S» 2»PXJ) 

PKP  *  PKP*<l.-TJX*PXJ)**NB 
97  CONTINUE 
GO  TO  100 
^  DO  96  J  *  ltN 

FJXY  *  SORT ( XUAIJ )*A2*YU3(J)**2)/S5 
X  *  FJXY  -  SH3 

Y  *  FJXY  -  SWA 

PKP  *  PKP  ♦  (1.-  CA5*(BlGG(X)-8iGG(Y))) 

96  CONTINUE 
GO  TO  100 
5  00  95  J  *  l,N 

Y  =  SORT ( XUA( J )**2+YUB( J )**2)/$5  -  SWA 
PKP  »  PKP  *  (1.-CA5M  l.-SIGG(Y) )  > 

95  CONTINUE 
GO  TO  100 

101  AL  «  SORT (ALi*AL2) 

DO  06  J  =  UN 

FJXY  =  SQRTIXUAI J)**2*YUB< J)**2)/AL 
X  =  FJXY  -  SW3 

Y  *  FJXY  -  SWA 

PKP  *  PKP*<  1  l.-CA5*IBIGG«X»-BIGG(Ym*ANB) ) 

86  CONTINUE 
GO  TO  100 

1011  *  SORT ( AL 1*AL2) 

DO  83  J  «  UN 

Y  *  SORT I  XUA(  J )**2*YU8U  )**2!/AL  -SWA 

PKP  *  PKP*(  l.-R*{  l.-(  l.-CA5*(l.-BlGGIYm*ANBJ) 

H3  CONTINUE 
GO  TO  100 
103  DO  9A  J  *  UN 
X  =  XUAIJ) 

Y  *  YUBUl 

CALL  SFIX,FCL3,S3,3MF) 

SFA  *  SKF 

CALL  SF( Y» FCLAf SA» SMF ) 

PKP  »  PKP  ♦  (l.-CA5*$FAASNF) 

9A  CONTINUE 
GO  TO  100 
10A  DO  93  J  =  l,N 
X  =  XUAIJ) 

Y  =  YUBIJI 

CALL  SF<X,GCL3,ALl»5MF) 

SFA  *  SHF 

CALL  SFIY,GCLA,AL2»SMF ) 

PKP  *  PKP*( l.-R*( l.-I U-CA5*SFA*SHF)**N8> ) 

93  CONTINUE 
GO  TO  100 
105  00  92  J  *  UN 
X  =  XUAIJ) 

Y  *  YUMJ) 

CALL  SFIX,FCL3.S3tSMF) 

SFA  *  SMF 

CALL  SF(Y,FCLA,SA,SMF > 

PKP  *  PKP  *  II.-  CA5* SF A*SMF ) 

92  CONTINUE 
GO  TO  100 
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106  00  91  J  «  UN 
XTA  *  XUA1J) 

YTB  «  YUB( J) 

CALL  FF{XTA,BUGCL3»0.eAFl 
AFl  *AF 

CALL  FF(YTB»B2fGCLAf0.»AF) 

PKP  »  PKP  *  U.-R*(U-<U-SP*AFt*AF)**NBn 
91  CONTINUE 
GO  TO  100 

107  DO  90  J  «  UN 
XTS  «  XUAU1/S3 
YTS  «  YUBUl/SA 

CALL  A113(ESUES2«X?S«YTS»P113) 

PKP  «  PKP  *  ( U-  CA5*Pl 131 
90  CONTINUE 
GO  TO  100 

108  TA1  «  ET1/AL1 
TA2  «  ET2/AL2 
00  89  J  *  UN 
XTA  »  XUAtJl/ALl 
YTB  »  YUB( JI/AL2 

CALL  AU5(TAUTA2tXTAtYTB»P113) 

PKP  «  PKP*( U-R* C 1 I .-CA5*Pl 1 3) •*NB> 1 
89  CONTINUE 
GO  TO  100 

109  00  88  J  *  1,N 
XTA  »  XUA( JJ/S3 
YTB  *  YUB( J1/S4 

CALL  A113(6SUES2,XTA,YTB,P113» 

PKP  «  PKP  ♦  (1.  -  C*5*P1 1 3 ) 

88  CONTINUE 
GO  TO  100 

110  00  87  J  *  UN 
XTA  *  XUA(J)/S3 
YTB  *  YUB1J1/SA 

CALL  AU3(ESlffcS2fXTA,YTB»Pim 
PA  *  P113 

CALL  AU3(WSUNS2vXTAtYTB*PU3) 

PKP  «  PKP*  I U-C45*IPA~P1131 ) 

87  CONTINUE 
GO  TO  100 

111  TAl  *  ETl/ALl 
TA2  *  ET2/A'  2 
HA l  «  Ml/ALl 
WA2  >  W2/AL2 
DO  85  J  «  UN 
XTA  «  X'JA1J)/AL1 
YTB  a  YUBI J1/AL2 

CALL  A113(TAUTA2'XTA,YT8«PU3> 

PA  *  PI 13 

CALL  All3(WAUWA2'XTAfYT6,PU3> 

PKP  »  PKP  *  (U-R*<U-U.-C*5*(PA-Pmn**N0>) 
85  CONTINUE 
GO  TO  100 

112  00  84  J  *  1»N 
XTA  *  XUA1JJ/S3 
YTB  *  YUBfJl/SA 

CALL  All3(ESUES2,XTAlYTBtP113) 

PA  *  PI 13 


o  o  o  u 
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CALL  Al 13 ( MS l ,  WS2 ,XTA«YTBtP113J 
PKP  *  PKP  *  U.-  C45*<PA-P113U 
84  CONTINUE 
100  PK6  *  l •—PKP 
RETURN 
END 


SUBROUTINE  All 3UZ.82 .  XX,  YY, PI  13 ) 

01  MENS  SON  DATA! 1200) 

DIMENSION  A22I2Q0) *822(200) 

COMMON  OATA 

COMMON  SE4,SE3,U3,U4,2,SQi,SQ2,Ql,Q2,CL2,Pl»SPI,SQ3.C,CL3,AU,AT2 
COMMON  BSl,8S2,SE12,SE22,S5,C45,SW3,SW4, ST1,SUS?SU4,AL1,AL2,SWI 
COMMON  FCL3,FCL4,ESI,ES2,ER1,ER2,ETB1,ETB2,WR1,WR2?WB1,WB2,WS1,WS2 
COMMON  A22«822»A1, A2,B1,B2 

COMMON  SQR2* SP 12  «P 12 « AJMX, BJMX*  TU1»TU2,  T 12t  T22 
COMMON  R32»RD1 ,RD2 , SPR , SH12 • SW22 , SRR 
COMMON  SN100*SNl01fSN102tSN110fSN112»X112 

COMMON  X100fX101,Xl02,X103,X104,X105,Xl06,XlQ7,Xl08,X109,X110,Xlil 
COMMON  N100fN101*N102»N110«Nll2t INUM 
COMMON  KK,K0UT,K00,X0,KH,LU3,LU4,N,NB,LU5 
EQUIVALENCE  I 0ATAI66 I ,UU5» , (DATA! 671 , VV5> 

E5  *  2.*PI/UU5 
SK  *  VV5 
PSUM  *  0. 

2  00  5  J  *  1 ,  LU5 
IFU.EQ.l)  GO  TO  1 
SK  *  SK  ♦  E5 
1  SSK  «  SIN(SK) 

COSK  *  COSISKJ 

2K  «  62  ♦  XX  *  COSK  ♦  AZ  *  YY  *  SSK  ♦  A2  *  BZ 

CK  * ( AZ  ♦  COSK  ♦  XX ) **2  ♦  ( BZ  *  SSK  ♦  YY)**2 

XY2J  *  21 13 ( CK ) 

PSUM  *  PSUM  ♦  2K  *  XY2J/UU5 
5  CONTINUE 
10  PI  13  *  PSUM 
RETURN 
END 


FUNCTION  BIuGI T) 

S2  =  SQRT 12.) 

TA  *  T/S2 

BIGG  *  .5  *  (1.  ♦  ERFITAH 
1  RETURN 
ENO 


SUBROUTINE  OECROIDATA) 

OECIMAL  CARD  REAO  SUBROUTINE 

MAIN  PROGRAM  MUST  CONTAIN  A  DIMENSION  STATEMENT  AND  CORRESPONDING 
EQUIVALENT  STATEMENT • 

ARRAY  DATA  MUST  BE  FIRST  ENTRY  IN  BLANK  COMMON  IN  PROGRAM. 
DIMENSION  0RBUI5),  DATA! 1 1 ,  i 1(5) 

9  REAO  (5,101  (Il(I)vI*lt5lt 1 1STL *  I ORBUI J) , J* 1 , 51 
J  *  IABSU1STI.) 

DO  30  1*1,5 

IF  IIIID.NE.OI  GO  TO  30 
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20  DATA(J)  *  DR8UU) 

30  J  *  JH 

IF  (I1STL)  50, 40,9 

40  WRITE  (6,900)  (Il( F ),  1*1,5), IlSTL,t DRBUf J ) » J*1 , 5) 

CAU  EXIT 
50  RETURN 

10  FORMAT  (5Il,I7,5FI2.m 

900  FORMAT  ( 1HI , 86H***A0DRESS  PORTION  OF  !P  CARO  *  ZERO*  RUN  TERMINAT 
1E0.  CARD  IN  ERROR  PRINTED  8EL0W.***/lH0,5Il *I?*5F12*ll) 

END 


SUBROUTINE  FF(OX,OY,OL,S,AF) 
AF  *  0. 

XPY  *  OX  ♦  OY 

XMY  *  OX  -  OY 

0L2  *  2.*  OL 

CALL  EE(XPY,OL,S,EEE) 

El  *  EEE 

CALL  EE ( XMY, 01. ;S,EEE ) 

E2  *  EEE 
AF  *  E1-E2 
RETURN 
END 


SUBROUTINE  H(YX,YL,HXL) 

HXL  *  0. 

IF(YL.EQ.O.)  GO  TO  1 
XPL  *  YX  »  YL 
XML  *  YX  ~  YL 

HXL  *  (XPL*BlGG(XPLJ-XML*8IGGtXMLKSMG(XPL»-SNG(XMLn/(2.*YL> 
GO  TO  2 

I  HXL  *  BIGG(YX) 

?  RETURN 
END 


SUBROUTINE  EE(X,CL,S,EEE) 
IF(S.fQ.O-)  GO  TO  l 
XS  *  X/S 
SL  *  CL/S 

CALL  H( XS , SL  oHXL ) 

EEE-  HXL 
GO  TO  5 

1  AL  *  -CL 
IF(X.LT.AL)  GO  TO  2 
IF(X,GE.CL)  GO  TO  3 
EEC*  ( X+CL ) / ( 2»*CL ) 

GO  TO  5 

2  ELE=  0. 

GO  TO  5 

3  EEE=  l. 

5  RETURN 

END 


SUBROUTINE  SF (OX,OY, S, SMF ) 
IF ( S.fcO.O. )  GO  TO  1 
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xpy  »(ox  ♦  uy>/s 

XMY  »IOX-OY)/S 

SHF  *  BIGG(XPY)-BISGIXHY) 

GO  TO  5 

1  IF(ABS(OX) . LE.OYJ  GO  TO  2 
SHF  »  0. 

GO  TO  *> 

2  SHF  *  1. 

5  RETURN 

END 


SUBROUTINE  HS(OX,OL,SH) 
IFIOL.EQ.O.)  GO  TO  1 
CALL  SF(OX»OLt l.tSHF) 

SH  «  SHF/(2.*OL) 

GO  TO  2 

1  SH  *  SHGIOX ) 

2  RETURN 
END 


FUNCTION  Z113ISO) 

IF ( S0<LE<<001 )  GO  TO  1 
£113  «  ( 1.-EXPI-S0/2. »)/S0 
GO  TO  2 

1  £113  »  .5  -  SO/8.  ♦  S0**2M8, 

2  RETURN 
END 
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5.8.  TEST  CASES  FOR  A  RIPPLE  OF  BOMBS  AND  FOR  A  RIPPLE  OF  FIXED  DISPENSERS 


Test  Case  1 


1  23  1  1  0  0  2  0 


NO  N 1 

N2  N  3 

N4  N5 

N10  Nil 

N12  N13 

N14  N15 

1  0 

0  0 

0  0 

0  1 

0  0 

0  0 

0 

OF 

SP 

N 

Nft 

NVI 

NN4 

40.00 

20.00 

0.95 

24 

1 

0 

3 

A3 

A4 

83 

84 

SI 

$2 

S3 

S4 

400.0000 

200.0000 

1.0000 

1.0000 

30.0000 

20.0000 

0.0 

0.0 

T  1 

T2 

SU1 

SU2 

U1 

U2 

LI 

00 

150.0000 

100.0000 

0.0 

0.0 

8.0000 

8.0000 

0.0 

0.0 

R  1 

R2 

01 

02 

R 

Ml 

M2 

15.0000 

30.0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 

ET1 

ET2 

GCL3 

GCL4 

ilU5 

VV5 

0.0 

3.0 

0.0 

0.0 

40.GC00 

0.0 

Am 

-60 .0000 

-60.0000 

-60.0000 

-60.0000 

-60.0000 

-60.0000 

-20.0000 

-20.0000 

-20.0000 

-20.0000 

-20.0CQO 

-20.0000 

20.0030 

20.3000 

20.0000 

20.0000 

20.0000 

20.0000 

60. COCO 

60.0000 

60.C000 

60.0000 

60.0000 

60.0000 

em 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

K0  *  11  0.122739 

X(  101 

xi  in 

Xi  121 

Xi  13) 

X  i  14) 

X  ( 15 ) 

X(30) 

0.0 

0.1227 

39  0.0 

0.0 

0.0 

0.0 

0.0 

XI100) 

xt  ion 

Xi 102) 

XI103) 

XI 104) 

X (105) 

XI106I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

Xi 107  ) 

Xi 1081 

Xi  109) 

XtllO) 

X(lll)  X (112) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

*(11)  is  the  expected  fraction  damage  to  an  area  target  from  a  ripple  of  iV(l)  bombs, 
each  subject  to  a  gaussian  ballistic  error  and  the  whole  ripple  subject  to  a  gaussian 
aiming  error.  The  area  is  rectangular  and  oriented  along  and  perpendicular  to  the 
line  of  flight.  The  major  damage  is  fragmentation.  *(11)  uses  HI  and  R2  entries  and 
53-B4-1. 
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Test  Case  2 


2 

23 

1  1 

0  0 

2  0 

NO  N 1 

N2  N3 

NA  N5 

N10  Nit 

N12  N13 

NIA  N15 

1  0 

0  0 

0  0 

0  0 

0  0 

0  1 

0 

Of 

Sf> 

N 

NB 

NV1 

NNA 

AO. 00 

20.00 

0.95 

2A 

1 

0 

3 

A3 

AA 

83 

BA 

SI 

S2 

S3 

SA 

A00.0000 

200.0000 

20.0000 

20.0000 

30.0000 

20.0000 

0.0 

0.0 

T1 

T2 

SUi 

SU2 

U1 

02 

LI 

00 

150.0000 

100.0000 

0.0 

0.0 

8. 0000 

8.0000 

0.0 

0*0 

R  1 

R2 

01 

02 

R 

Ml 

M2 

1.0000 

1 .0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 

EM 

ET2 

GCL3 

COLA 

UU5 

VV5 

0.0 

0.0 

0.0 

0.0 

AO. 0000 

0.0 

A(J) 

-60.0000 

-60.0000  - 

60.0000 

-60.0000 

-6O0COOO 

-60.0000 

-20.0000 

-20.0000 

-20.0000 

-20.0000  - 

20.0000 

-20.0000 

20.0000 

20.0000 

20.0000 

20.0000 

20.0000 

20.0000 

60.0000 

60.0000 

60.0000 

60.0000 

60.0000 

60.0000 

BiJ) 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.000',) 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20.0000  - 

20.0000 

-0.0000 

20.0000 

-20.0000 

-0.0000 

20,0000 

KO  «  IS  0.0AH656 

XC 101 

XI  111 

XI  121 

XI 13) 

XllA) 

X  1 15 1 

X  (30) 

0.0 

0.0 

0.0 

0.0 

0.0 

O.OA8656  0.0 

xnooi 

XI  101  » 

XI 1021 

XI 103) 

XllOA) 

X 1105) 

XI106) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X( 107t 

XI 108) 

XI 109) 

XI110) 

Xllll) 

X  €  1 1 2 ) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X(15)  Is  the  same  method  of  delivery  as  tf(ll)  but  the  damage  is  based  on  Impact  (cooVie- 
cutter  damage  functlor).  X(15)  uses  B3  and  BA  entries  and  B1-B2-1. 


U 
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Test  Case  3 


3  23 

1 

1  0 

0  1 

0 

NO  N 1 

N2  N  3 

N4  N5 

N10  Nil 

N12  N13 

N14  N15 

0  1 

0  0 

0  0 

0  0 

0  0 

0  0 

0 

DH 

SP 

N 

NB 

NVI 

NN4 

20.00 

0.0 

0.95 

4 

1 

0 

1 

A3 

A4 

B3 

B4 

SI 

S2 

S3 

S4 

400.0000 

200.0000 

1.0000 

1.0000 

30.0000 

20.0000 

0.0 

0.0 

Tl 

T2 

SU1 

SU2 

U1 

112 

LI 

QO 

150.0000 

ICO. 0000 

0.0 

O.C 

8.0000 

a. oooo 

0.0 

0.0 

R1 

R2 

01 

02 

R 

Wl 

W2 

15.0000 

30.0000 

l.nooo 

1.0000 

0.0 

c.o 

O.C 

EU 

ET2 

GCL3 

CCL4 

UU5 

VV5 

0.0 

0.0 

o.c 

0.0 

40.0000 

0.0 

A(J  ) 

-30.0000  -10.0000  iO.OOCO  30.0000 

ALL  BIJ)  *  0. 


mil 

XI  2 1 

XI  31 

XI4) 

XI5) 

0.031374 

O.C 

0.0 

0.0 

0.0 

XC  101 

XI  111 

XI  12) 

XI 1 31 

X 1 141 

X  1 15 ) 

XI 301 

0.0 

0.0 

o.c 

0.0 

0.0 

0.0 

0.0 

X( 1001 

XI 1011 

XI 102) 

XI 103) 

XI 104) 

XI 105) 

XI 106) 

O.C 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X< 1071 

XI 1081 

XI  109) 

XIUO) 

Xf*ll» 

x 11121 

0.0 

0.0 

0.0 

0.0 

0.0 

C.O 

.K(l)  is  the  same  as  X(ll).  The  method  of  computation  is  different  and  is  faster 
than  X(ll)  provided  This  restriction  has  been  built  in  to  the  program. 
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Test  Case  4 


4 

23 

1  1 

0  0 

1  0 

NO  N1 

N2  N3 

N4  N5  N10  Nil  N12  N13  N14  N1S 

0  0 

0  0 

0  1 

0  0 

0  0 

0  0 

0 

OF 

SP 

N 

NB 

NVI 

NN4 

20.00 

0.0 

0.95 

4 

1 

0 

1 

A3 

A4 

83 

84 

SI 

S2 

S3 

400.0000 

200.0000 

20.0300 

20.0000 

30.0000 

20.0000 

0.0 

T 1 

T2 

SU1 

SU2 

U1 

U2 

LI 

ISO. 0000 

100.0000 

0.0 

0.0 

8.0000 

8.0000 

0.0 

R 1 

R2 

01 

02 

R 

Ml 

M2 

1 .0000 

1.0000 

1.0000 

l.OCOO 

0.0 

0.0 

0.0 

ET1 

ET2 

GCL3 

SCL4 

UU5 

VV5 

0.0 

0.0 

0.0 

0.0 

40.0000 

0.0 

AC  J) 

-30.0000 

-10.0000 

10.0000 

30.0000 

ALL  Si J  > 

■  0. 

XIII 

X(  2) 

X(  3) 

X(4) 

X(5) 

0.0 

0.0 

0.0 

0.0 

0.009736 

X(10) 

X<  Hi 

X(  12) 

X(  13) 

XI14) 

XflS) 

X  1 30  ) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X( toot 

X(101) 

X(  102) 

X( 103) 

XI 104) 

XI 105) 

X ( 106) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

XI 107) 

X( 108) 

X(  109) 

XI110) 

xmi) 

XI112I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1(5)  is  the  same  as  .Y (15) .  It  should  be  restricted  to  //(1)£4. 
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Test  Case  5 


5  23 

1 

1  1 

0  1 

0 

NO  N 1 

N2  N3 

N4  N5 

N10  Nil 

N12  N1 3 

N14  N15 

1  0 

0  0 

0  0 

0  0 

1  0 

0  0 

0 

OF 

SP 

N 

NB 

NVI 

NN4 

0.0 

30.00 

0.95 

2 

100 

15 

1 

A3 

A4 

83 

84 

SI 

S2 

S3 

S  4 

800.0000 

200.0000 

1.0000 

1.0000 

1.0000 

20.0000 

0.0 

0.0 

T 1 

•  T? 

SU1 

SU2 

111 

U2 

LI 

QO 

150.0000 

100.0000 

0.0 

0.0 

8.0000 

8.0000 

1.0000 

0.0 

R 1 

R2 

01 

02 

R 

HI 

W2 

7.5000 

15.000C 

1.0000 

1.0000 

0.0 

0.0 

0.0 

ET1 

ET2 

GCL3 

GCL4 

UU5 

VV5 

0.0 

0.0 

0.0 

0.0 

40.0000 

0.0 

svim 

0.0  100. 0000  200.0000  300.0000  400.0000  500.0000  600.0000  700.0000 

800.0000  900.0000  1000.0000  1100.0000  1200.0000  1300.0000  1400.0000 

BPUJ) 


0.0 

0.1000 

0.0200 

0.1000 

0.0400 

0.1000 

0.0600 

0. 0800 

0.0800 

0.0600 

0.1000 

0.0400 

0.1000 

0.0200 

ALL  At J  )  « 

0. 

Bt  J  1 

0.0 

0.0 

KO  -  12 

0.141609 

X(  10 » 

X(  11) 

Xt  12) 

XI13) 

X  f  14) 

X  1 15 ) 

Xf30) 

0.0 

0.0 

0.141609  0.0 

0.0 

0.0 

0.0 

X<  100) 

X( 101) 

Xf  102) 

Xf 103) 

X1104) 

XI 105) 

X 1 106) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X(  107) 

XI 108) 

X41 09) 

XtUO) 

Xflll) 

Xf 112) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X(12)  gives  the  expected  fraction  damage  to  an  area  target  from  a  ripple  of  fixed 
dispensers,  where  the  ballistic  dispersion  of  the  bomblets  in  range  is  given  by  a 
table.  The  deflection  dispersion  is  gaussian,  as  is  the  aiming  error.  The  major 
damage  effect  is  fragmentation.  Rh,  R2,  CL  1,  $0,  .92,  and  B 3*/i4*l  are  required  inputs. 
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Test  Case  6 


6  23  1  1  1  0  1  0 


NO  N 1 

N2  N3 

N4  N5 

N10  Nil 

N12  N13 

N14  N15 

1  0 

0  0 

0  0 

0  0 

0  1 

0  0 

0 

DF 

SP 

N 

N8 

NVI 

NN4 

0.0 

30.00 

0.95 

2 

100 

15 

1 

43 

A4 

83 

84 

SI 

S2 

S3 

S  4 

800.0000 

200.0000 

1.0000 

1.0000 

1.0000 

20 .0000 

0.0 

0.0 

T  1 

T2 

SU1 

SU2 

U1 

U2 

LI 

00 

150.0000 

100.0000 

0.0 

0.0 

8.0000 

8. 0000 

1.0000 

0.0 

Ri 

R2 

01 

02 

R 

Ml 

W2 

7.5000 

15.0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 

ET1 

ET2 

GCL3 

GCL4 

uus 

VV5 

0.0 

O.C 

0.0 

0.0 

40.0000 

0.0 

SVKJ) 

• 

0.0 

100.0000 

200.0000 

300.0C00 

400.0000 

500.0000 

600  .0000 

700.0000 

800.0000 

900.0000 

1000.0000 

1100.0000 

1200.0000 

1300.0003 

1400.0000 

8PUJ) 

0.0 

0.0200 

0.0400 

0.0600 

0.0800 

G  « 1000 

C.1000 

0.1000 

0.1000 

0.1000 

0.1000 

0.0800 

0.0600 

0.0400 

0.0200 

ALL  A( J I  -  0. 


BUI 

C.O  0.0 


KO  «  13  0.139480 


X(10) 

XC  11) 

XI 12) 

X(  13) 

X  ( 14) 

X  ( 15 ) 

X  ( 30 ) 

0.0 

0.0 

0.0 

0.139480 

0.0 

0.0 

0.0 

Xt  100) 

X< 101) 

X(  102) 

X( 103 ) 

X(104) 

X  ( 105 ) 

X ( 106 ) 

0.0 

0.0 

0.0 

0.0 

C.O 

0.0 

0.0 

X',  107) 

X( 108) 

X( 109) 

X(110) 

Xlllll 

X ( 1121 

C.O 

0.0 

0.0 

0.0 

0.0 

0.0 

AT (13)  is  an  approximation  to  A(12) 
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Teat  Case  7 


7 

23 

1  1 

1  0 

*<0  Ml 

L  0 

N2  N3 

0  0 

N4  N5 

0  0 

N10  Nil  1 

n  o 

0 

0.0 

OF 

30.00 

SF 

0.95 

N 

2 

A3 

000.0000 

A4 

20C.0000 

83 

15.0000 

84 

15.0000 

T 1 

150.0000 

T2 

100.0000 

SU1 

0.0 

S02 

0.0 

R1 

1.0000 

R2 

1.0000 

01 

l.COOO 

02 

1.0000 

ET1 

0.0 

ET2 

0.0 

GCL3 

0.0 

GCL4 

0.0 

SVIIJ) 

0.0  l'JO.OOOO  200.0000  300.0000  ‘ 

SCO. 0000  900.0000  1000.0000  U00.G00C  l; 

BPIUI 

0.0  0.0200 

0.1000  0.1000 

0.0400 

0.1000 

0.0600 

0.0800 

ALL  Aid) 

»  0. 

BIJI 

0.0 

0.0 

KO  ■  19  0.099030 

X(  101 

XI  111 

XI  121 

XI 131 

0.0 

0.0 

0.0 

0.0 

XI 1001 

X( 101) 

XI 1021 

XI 1031 

0.0 

0.0 

0.0 

0.0 

X( 1071 

X( 1061 

XI 1091 

XI1I0I 

0.0 

0.0 

0.0 

0.0 

0  1  0 


N12  N13 

N14  N15 

0  0 

1  0 

N8 

NVI 

NN4 

100 

15 

1 

SI 

S2 

S3 

S4 

1.0000 

20.0000 

0.0 

0.0 

U1 

1)2 

LI 

QO 

0.0000 

e.ooco 

1.0000 

0.0 

R 

Ml 

M2 

0.0 

c.o 

0.0 

uus 

VV5 

40.0000 

0.0 

400.0000 

500.0000 

600.0000 

700.0000 

1200.0000 

1300.0000 

1400.0000 

0.0800 

0.1000 

0.1000 

0.1000 

0.0600 

0.0400 

0.0200 

XI14) 

X  115  ) 

XI 30 1 

0.099638 

0.0 

0.0 

XI104I 

X 1 1051 

XI 106) 

0.0 

0.0 

0.0 

Xlllll 

X 1112} 

0.0 

0.0 

#(14)  la  for  the  same  type  of  delivery  os  #()2)  but  the  damage  effect  is  based  on  im¬ 
pact  (cookie-cutter  damage  function).  Required  inputs  are  53,  54,  and  51»/?2«1. 
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Test  Case  8 


8 

23 

1  1 

1  0 

1  0 

NO  N 1 

N2  N  3 

N4  N5 

N10  Nil 

N12  Nl  5 

N14  .Ml  5 

0  0 

1  0 

0  0 

0  0 

0  0 

0  0 

0 

OF 

SP 

N 

NB 

NVI 

NN4 

C.O 

0.0 

0.95 

1 

30 

15 

1 

A3 

A4 

B3 

B4 

SI 

S2 

S3 

54 

800.0000 

200.0000 

l.OCOO 

1.0000 

l.OOOC 

20.0000 

0.0 

0.0 

n 

T2 

SO! 

SU2 

U1 

112 

LI 

00 

150.0000 

100.0000 

0.0 

0.0 

a. coco 

8 .0000 

1.0000 

0.0 

Rl 

R2 

01 

02 

R 

W1 

W  2 

7.5000 

15.0000 

1.0000 

l.OCOO 

0.0 

C.O 

0.0 

EU 

ET2 

GCL  3 

GCL4 

Ul)5 

VV5 

0.0 

0.0 

0.0 

0.0 

40.0000 

C.C 

svim 

0.0 

100.0000 

200.0000 

300.0000- 

400.0000 

500.0000 

600.C000 

700.0000 

SCO. 0000 

900.0000  lOOO.COOO  U00.C000  1200.0000  1300.0000  1400.0000 

BP  I!  Jl 

0.0 

0.0200 

0. 0400 

0.0600 

0.0800 

C.1000 

0.1000 

0.1000 

0.1000 

O.ICOO 

0.1000 

0.0800 

0.0600 

0.0400 

0.0200 

ALL  A{ J  ) 

•  0. 

ALL  BUI 

»  0. 

XIII 

XI21 

XI  3) 

XU) 

X 15 ) 

0.0 

0.028578  0.0 

0.0 

0.0 

X(10) 

XI  111 

XI  121 

XI 131 

XI14) 

X  1 15 ) 

XI 30) 

0.0 

0.0 

0.0 

C.O 

0.0 

0.0 

0.0 

XI 100) 
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X(2)  Is  the  same  as  X(12)  except  for  the  method  of  computation.  It  applies  only  to 
one  dispenser  and  is  valid  for  a  small  number  of  bomblets  (NB< 30). 
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Test  Case  9 


<5 

'  23 

1  1 

1  0 

1  0 

NO  N 1 

N2  N  3 

N4  N5 

N10  Nil 

N12  N13 

N14  N15 

0  0 

0  1 

0  0 

C  0 

0  0 

C  0 

0 

OF 

SP 

N 

NB 

NV1 

NN4 

0.0 

0.0 

0.95 

1 

30 

15 

1 

A3 

A4 

B  3 

B4 

SI 

S2 

S3 

S4 

800.0000 

200.0000 

1.0000 

1.0000 

1.0000 

20.0000 

O.C 

0.0 

T  1 

r2 

SU1 

SU2 

01 

U2 

LI 

QO 

190.0000 

100.0000 

0.0 

O.C 

8.0000 

6.0000 

1.0000 

0.0 

R 1 

R2 

01 

02 

R 

Ml 

W2 

7.5000 

15.0000 

1.0000 

1.0000 

0.0 

0.0 

0.0 

ET1 

ET2 

GCL  3 

GCL4 

UU5 

VV5 

C.O 

0.0 

0.0 

0.0 

40.0000 

0.0 

SVKJi 

0.0 

100.0000 

200.G000 

300.0000 

400.0000 

500.0000 

600.0000 

700.0000 

SCO. 0000 

900.0000  IGOO.OOOC  1100. 0000  1200.0000  1300.0000  14C0.0000 

bp  in » 

0.0 

0.0200 

0.0900 

0.0600 

0.0800 

0.1000 

0.1000 

0.1000 

0.1000 

0.1000 

0.1000 

0.0800 

0.C600 

0.0400 

0.0200 

ALL  A(J  1 

*  0. 

ALL  BUI 

»  0. 

X(  11 

X*  21 

X(  3) 

X(4) 

X(5> 

0.0 

0.0 

0.0279C4  0.0 

0.0 

X(  lot 

X(  111 

X<  121 

X(  13) 

X  <  14 1 

X  (15) 

X(  30  1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 

X( 1001 

x.(  ion 

x< 1021 

Xf 103) 

XI 104) 

X  ( 1C  5 ) 

X ( 106) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

X< 107) 

X(  10B) 

XI  109) 

XfllO) 

X(lll) 

X (112) 

0,0 

0.0 

0.0 

0.0 

0.0 

0.0 

X(3 )  <s  the  same  as  *(13)  except  for  the  method  of  computation.  It  applies  only  to 
one  dispenser  and  is  valid  for  a  small  number  of  bomblets  (NB< 30). 
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Test  Case  10 
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tf(4)  is  the  same  as  except  for  the  method  of  computation.  It  applies  only  to 

one  dispenser  and  is  valid  for  a  small  number  of  bombiets  (NB< 30) . 
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5.9.  TEST  CASES  FOR  A  RIPPLE  OF  DISPENSERS,  RECTANGULAR  BOMBLF.T  PATTERN 

The  series  ,Y(103)  through  A'(106)  applies  to  ripples  of  dispensers  released  from  the 
carrier,  which  provide  a  rectangular  pattern  of  bomblet  impact  points.  AT ( 1 03 )  ap¬ 
plies  to  the  case  in  which  the  major  damage  effect  is  frcm  fragmentation,  while 
A'(105)  is  for  the  cookie-cutter  damage  effect  case.  The  ballistic  dispersion  is 
taken  into  account  for  these  cases.  If  the  MAE  of  an  individual  bomblet  is  appreci¬ 
able  (say,  at  least  a  tenth)  of  the  total  dispenser  pattern  area,  the  edge  effects 
(effect  of  a  bomblet  outside  the  pattern  area)  become  important.  In  these  cases, 
A'(!0A)  is  used  for  the  fragmentation  case  and  Y(106)  for  the  cookie-cutter  case. 


In  both  the  latter,  the  dispenser  ballistic  dispersion  is  Ignored. 


i’  1  =£>"2  =5  3  =S*i  “0 . 
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Test  Case  11 
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Test  Case  12 
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Test  CasJ  13 
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Teat  Case  14 
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5.10.  TEST  CASES  FOR  A  RIPPLE  OF  DISPENSERS.  ELLIPTIC  BOMBLET  PATTERN 

(107)  through  If (109)  apply  to  ripples  of  dispensers  released  from  the  carrier, 
which  provide  an  elliptic  pattern  of  bomblet  impact  points.  jf(107)  applies  to  the 
case  in  which  the  major  damage  effect  is  fragmentation,  while  Af(109)  applies  to  the 
cookie-cutter  type  damage  effects.  The  ballistic  dispersion  is  taken  into  account. 
If  the  MAE  of  an  individual  bomblet  is  appreciable  compared  with  the  total  dis¬ 
penser  pattern  area,  the  edge  effects  (effect  of  a  bomblet  outside  the  dispenser 
pattern  area)  become  Important.  In  this  case,  *(108)  is  used  for  a  fragmentation 
effect.  For  the  cookie-cutter  effect,  one  may  use  an  equivalent  rectangular  pat¬ 
tern  (equal  area  and  ratio)  in  *(106)  above. 
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Test  Case  15 
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Test  Case  16 
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Test  Case  17 
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5.11.  TEST  CASES  FOR  A  RIPPLE  OF  DISPENSERS,  ELLIPTIC  ANNULUS  BOMBLET  PATTERN 


£(110)  through  £(112)  are  the  seme,  respectively,  as  £(107)  through  £(109)  except 
that  the  dispenser  pattern  is  an  elliptic  annulus.  £(110)  applies  to  the  fragmen¬ 
tation  case  while  £(112)  applies  to  the  cookie-cutter  case,  £(111)  is  the  f ragmen' 
tation  case  when  the  edge  effects  must  be  taken  into  account. 
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Test  Case  18 
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Test  Case  19 
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Test  Case  20 
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5.12.  TEST  CASES  FOR  A  RIPPLE  OF  DISPENSERS «  CIRCULAR  ANNULUS  BOMBLET  PATTERN 

#(100)  rhrough  ,Y(102)  are  the  same,  respectively,  as  IfQlO)  through  Y (112)  except 
that  the  bomblet  patterns  are  circular,  further,  an  approximation  is  made  that  is 
valid  only  when  the  ratio  of  the  pattern  raaiuf  to  the  ballistic  error  is  suffi¬ 
ciently  large.  X(100)  applies  to  the  fragmentation  case  while  If (102)  applies  to 
the  cookie-cutter  case.  X(lOl)  is  the  fragmentation  case  when  edge  effects  must 


be  taken  into  recount. 
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Test.  v/a^e  21 
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Teat  Case  22 
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Test  Case  23 
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Appendix  A 
TRAIN  DISTRIBUTION 

Consider  a  train  of  N  weapons  delivered  in  a  pattern.  The  ballistic  centers 
of  impact  forming  the  pattern  are  at  points  (a. ,b .) ,i**l ,2, . . . ,N.  The  origin  is 

selected  auch  that  0.  Each  weapon  is  subject  independently  to  &  ballistic 

i.  i, 

error  whose  distribution  B  ie  assumed  to  be  of  the  form 

where  BiOf)  is  the  distribution  of  ballistic  errors  in  range  (a:)  and  il2(y)  is  c^e 
distribution  in  deflection  (y) .  Let  b\(x)  and  fc2(y)  be  the  corresponding  density 
function,  i.e., 


b\{x)dx. 


We  shall  assume  that  the  distributions  have  mean  zero  and  standard  deviations  a 


and  e  ,  respectively. 

y 

The  form  of  the  functions  b\  and  fc2t  the  values  of  b  and  a  ,  and  the  pattern 

x  y 

values  (a.,fc.)  must  be  determined  through  some  testing  procedure.  If  the  data 
i-  i 

are  sufficiently  complete  (individual  data  for  each  weapon  in  a  train),  the  indi¬ 
vidual  distributions  B\  and  B2  can  be  used  as  required  in  the  methods  of  Ref.  [1]. 

In  general,  Bj  and  #2  can  be  arbitrary  distributions.  For  normal  train  bombing, 

the  values  of  a.  and  b.  can  be  determined  from  the  delivery  conditions  and  inter- 
1  1 

valometer  settings.  For  the  methodology  of  Sec.  3,  it  is  necessary  that  B\  and  B 2 
be  gaussian  in  form,  so  that  only  0  and  0  are  needed.  The  spaci.-igs  a.  and  b. 

X  tj  "V  X- 

are,  of  course,  necessary  as  further  inputs.  However,  for  some  of  the  cluster- 
type  weapons,  the  complete  data  are  not  available.  Although  the  complete  data  are 
available  in  other  cases*  the  simplification  considered  here  is  usually  sufficient. 

Let  uo  assume  that  the  test  data  in  range  are  only  sufficient  to  estimate  the 
expected  number  of  weapons  landing  in  a  set  of  range  intervals.  Specifically,  for  a 
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set  of  K  Intervals  between  the  points  Xq,X\ ,*2, . , . ,X^t  we  obtain  the  corresponding 

y 

number  of  weapons  . where  the  total  number  is  M*l.  ,N..  We  form  the 

A  1*1  X 

percentage  distribution  within  intervals 


and  the  sample  density 


(A.i)  Vimx'3i — »  i*i»2 . K‘ 


We  shall  view  the  p ^  as  an  estimate  of  the  probability  density  function  p(x )  of  a 

J 

single  weapon,  where  ?(*)-(  p(«)dr  is  the  probability  that  an  individual  weapon 

J  —00 

landed  s  V-  Basically,  we  are  using  as  an  approximate  model  one  in  which  each 
weapon  is  randomly  drawn  from  the  positions  in  the  train.  We  therefore  assume 
that 


The  above  model,  in  effect,  assumes  that  each  weapon  in  the  pattern  is  aimed 
at  the  center  of  the  pattern.  For  symmetric  patterns,  the  center  is  estimated  as 
the  mean  point  X.  However,  the  data  are  often  quite  irregular,  no  that  some  smooth- 
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lng  must  be  done.  We  have  chosen  to  use  the  quartiles  X  ^  and  X  ^  as  character¬ 
istics  to  be  used  in  the  fitting  process ,  since  their  use  introduces  a  type  of 
smoothing.  The  center  X  then  is  defined  as  the  midpoint  between  the  quartiles 


For  symmetric  patterns,  XQmX .  Thus  X  25  and  X  75  are  defined  by  the  equations 


P(^25)-.25, 

K*.7s)-- ”• 

Define  I  and  J  by  the  conditions 


(A.  5) 


p(*i)<-25-  'tW'-25- 


Then  X  ^  and  X  75  are  given  by  the  equations 


where  P(X .)  is  given  in  (A. 2) . 

J 

Finally,  making  the  change  in  variable  x*X-X^  to  center  the  coordinate  sys¬ 
tem  on  the  pattern  center  and  using  (A. 4),  we  obtain  the  center  x and  the  quar¬ 
tiles  x  and  x 

a  2 J  • 75 
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x  «0, 
a 


^.75^.75-V 


X.25“*.25-V 


*.75_*.25 


^.25**^.  7  5 


L)-X.75* 


where  X  ^  and  X  ^  are  given  in  (A. 6).  Thus,  we  have  characterized  the  data  by 
three  characteristics:  x  ,  S,  and  x 

In  this  case,  we  do  not  have  enough  information  to  determine  che  form  of 

the  ballistic  distribution  Bi(X)  or  of  the  pattern  a..  Let  us  assume  that  the 

form  of  B i(X)  is  gaussian  with  variance  a2  and  that  the  a.  are  known.  Then,  for 

X  1 

a  delivery  of  N  weapons,  the  expected  number  falling  short  of  X  is 


N  /x-a.\  N  /x-a.\ 

.1  •.M-J'hr} 

VI  \  X  /  VI  \  x  / 


and  the  percentage  of  weapons  expected  to  fall  short  of  X  is 


if 

\  x  / 


The  function  P(x)  is  an  approximation  to  this  function,  fi(x) .  We  will  make  the 
assumption  that  the  N  weapons  are  evenly  spaced,  i,e., 

(A. 8)  a.m(2i-N-l)d, 

t 


where  2d  is  the  spacing.  We  have  characterized  F  by  three  parameters:  x  “0, 

a 

variance  S2,  and  upper  quartile  x  We  will  determine  values  for  d  and  a  under 

.  /o  x 

the  restriction  that  the  distributions  Q(x)  and  P(x)  have  the  same  variances  and 
quartiles.  Our  conditions  are  thus 


£-1 
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Since  x  ^  m  ~x  75,  the  latter  two  are  equivalent.  Subtracting  and  combining,  we 
obtain  the  equivalent  condition 


where  f(x,y)  is  defined  as 


f(x,y)-G(x+y)-G(x-y) , 


and  thus  /(-x,y)  ■  f(x,y).  In  summary,  our  conditions  are 


(A. 11) 


*2-%4  l 

x  i-1 


x  “0. 

Q 


i.e.,  d  and  e are  determined  from  these  conditions. 


Appendix  B 
STICK  DISTRIBUTION 


There  are  occasions  when  the  ballistic  dispersion  of  a  cluster  of  weapons 
may  be  approximated  by  what  will  be  called  the  "stick  distribution."  Consider  the 
problem  of  throwing  a  stick  at  a  line  perpendicular  to  the  stick  in  an  attempt  to 
cover  the  line.  Assume  a  stick  or  length  2L,  with  the  center  aimed  at  the  point 
0.  Assume  that  the  error  in  throwing  tne  stick  is  governed  by  a  gsussian  distri¬ 
bution  with  variance  e2.  The  probability  of  covering  the  offset  point  *  is  given 
by 


Prob (point  at  x  covered) 


where  g(x)m  exp(-x2/2)//2rT.  Define  the  function  h(.x,L )  as 


1 

(8. 1)  h(.x,l)~jjr f  g(y)dy. 

“x~L- 


The  function  h(.x,L)  can  be  viewed  as  a  probability  density  function,  since 


(B.  2) 


00 

f  h(x,L)dx~l. 

'•-or, 


We  will  call  the  distribution  function  H(x,L)  the  stick  distribution,  i.e.. 


(B.3) 


y+L 

g(z)d*dy. 

y-L 


More  generally,  we  have 


X/e 


(B.4) 


•IW  iW 
v  *•  00  “  00 


The  first  moment  of  H(X/a,L/e )  is 
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and  the  variance  ia 

(B.6) 


I>(^-»24 


Note  that  we  can  express  H  in  terms  of  G  and  g.  Using  Lemma  3  of  App.  D  and  (B.3), 
we  obtain 


(B.7)  #<*,£>- J  ■ftttHv+L)‘G(v-L)ldy 


[.(X+L)G(X+L)-(X-L)G(X-L)+g(X+L)-g(X-L)-) . 


The  stick  distribution  above  is  often  a  good  approximation  to  the  train  dis¬ 
tribution  P(x)  discussed  in  the  preceding  section.  As  in  App.  A,  assume  P(x)  to 
be  characterized  by  the  three  parameters  xfl*0,  variance  S2,  and  upper  quartlle 
x  (x  25*-£  75) .  Requiring  that  P  and  its  approximation  H  have  the  same  center 
coordinate,  variances,  and  quartlles,  we  obtain  the  conditions 


S2»s2+j-, 


where  S  and  x  ^  are  given  in  (A. 3)  and  (A. 7).  Subtracting  the  last  two  equations 
in  (B.8)  and  using  the  relation  x  25“"*  75*  we  ob£ain  411  alternate  condition 


Thus,  we  may  express  the  conditions  (B.8)  as 
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v°* 


where  the  function  F(x,y,L )  is  defined  as  in  (4.12)  as 


F(x,y  ,L)*H(x+y  ,L)-H(x-y  ,L)  . 
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Appendix  C 


THE  ELLIPTIC  COVERAGE  FUNCTION 


The  elliptic  coverage  function  P(R\ tR2\a,b)  is  defined  as  the  integral  of  a 
circular  gaussian  distribution  of  unit  variance  over  the  area  enclosed  by  an  off¬ 
set  ellipse  with  center  at  ( a,b )  and  axes  R\,R2  in  the  X  and  y  directions.  Thus, 


we  have 


where 


P(7?l  ,R2ia,b)=  f ]g(x)g(.y)dxdy, 
A  i 


-Kf fHth 


Making  the  change  of  variable  5**0  ,  r\mya  ,  we  obtain 

x  y 


i.e.,  an  equivalent  definition  is  the  integral  of  a  bivariate  gaussian  over  an 


offset  ellipse.  Similarly,  we  obtain 


\  x  y  x  y/  '  x '  '  y'  x  y 


'n  \d£,dn 


where 


Az- (x-a)2+(y-b)2£R2 . 


Another  definition  is  the  integral  of  a  bivariate  gaussian  over  an  offset  circle. 
Using  (C.l),  if  R\*R2mR  the  elliptic  coverage  function  becomes  the  circular  cover¬ 
age  function  P(fl,r),  i.e., 
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r 


K 


*  P(P,P;/2,2>)«p(p,  ^b5). 

From  Ref.  1,  alternate  expressions  for  P(P,r)  are 

P(P,r)-f  3exP-(~ -(|2+2-2-^  j0(ar) da 


(C.  A) 


mRfm£3±l&a&lIliSs)dMt 


where  Io(a)  flnd  are  Bessel  functions  of  the  first  kind  of  imaginary  argument. 

If  a*b’0,  the  elliptic  coverage  function  is  expressible  in  terms  of  the  cir¬ 
cular  coverage  function  (Ref.  7).  Thus,  from  Rtf.  2,  if  flj>/?2,  we  have 


(C.  5) 


(R\ +P2  R\~R?\  /R i~7?2  /? i+/?2\ 
P(P ! ,  P2 ;  0 , 0)  -P\-2"-*“l~/"P\“r"*~T“, /  * 


For  a*6,  7?i*/?2,  the  elliptic  coverage  function  has  been  evaluated  by  various 
computational  schemes.  Included  here  is  a  computational  scheme  well  adapted  for 
high-speed  computers  based  on  an  unpublished  paper  by  Oliver  Cross  at  Rand.  Using 
the  definition  in  (C.l)  and  converting  to  polar  coordinates,  we  obtain 


f 

(C.  6) 

i 

i 

R 

* 

is 

vth  :re 

V< 

I 

1 

(C.7) 

A  3 


Consider  first  the  case  where  the  origin  is  inside  the  ellipse  in  (C.7).  Then, 
we  obtain 


(C.8) 


-2t!(,7’(0)  ,  2. 

P(fl,,P2;a,fc)"jo  Jo  °.  -Wo 


'  ”  0 


where  t>(0)  is  the  radial  distance  at  angle  0  to  the  ellipse  E\  obtained  from 
(C.7)  as 


(C.9) 


i.e. , 


(C.10)  P(7?1,/?2;a,Z>)-~J 

E\ 

where  the  integration  is  around  the  boundary  of  the  ellipse  £j  given  by  (C.10). 
If  the  origin  is  not  inside  the  ellipse  in  (C.7),  P(7?j  ,/?2  ;<z,fc)  is  given  by 


»0j>  r2(0) 

(C.  11)  p'Rl  ?2ia,i/» j  j 


,.02^1(0) 

-  I  P  exp  I 

'Vo 


where  0j  and  02  are  the  angular  extremities  of  the  offset  ellipse,  r2(e)  is  the 
larger  radial  distance  to  the  ellipse  at  angle  0,  and  r*i(0)  is  the  smaller  radial 
distance,  as  shown  In  the  diagram  below. 
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Thus,  when  we  Integrate,  we  obtain 


(C.12)  ' 

where  r2(0)  is  the  radial  distance  to  the  ellipse  E\  of  (C.9)  on  the  far  side  from 
@1  to  02,  and  r2(C)  is  the  radial  distance  to  E\  on  the  rear  side  from  0j  to  02. 
Thus,  we  have 


(C. 13) 


exp('!^ri)]'i0. 


where  as  in  (C.10),  the  integration  is  around  the  boundaries  of  the  ellipse  E j  in 
(C.10).  Thus,  (C. 10)  holds  for  both  cases.  Let  us  make  the  transformation  on  0: 


(C. 14) 


r  cos6»f?i  cosyfa, 


r  sin0*P2  sin <p*-b, 


so  that 


r*2*^ (/?j  costfH-a)2+(P2  sincpffc)2, 


(C. 15) 


R?  sin«pfi 
SR i  coscphj’ 


/?l/?2+aj?2  cos«p+-fc/?i  simp 


Under  this  transformation,  the  offset  ellipse  E\  changes  Into  the  unit  circle 
with  center  at  the  origin.  Thus,  in  (C.10)  9  goes  from  0  to  2it,  so 


(C.  16)  P(/ri,H2;a,fc)=~J  ^^^&Q-{R\Rz+aRz  coscf*W?i  sincp)dcp, 


'0  r 


—14 1— 


where  r  is  given  in  (C.15).  Thus,  the  problem  of  evaluating  P(R\  ,/?2 \o.,b)  has 
beer  reduced  to  a  simple  quadrature  over  the  interval  (0,27i). 

Using  a  trapezoidal  type  of  numerical  integration  with  N  intervals  over  the 
range  (0,2n),  i.e,,  intervals  of  length  2n/N  and  the  midpoint  of  the  first  inter¬ 
val  at  cpg,  we  obtain 


(C. 17) 


P<»,  ,s2 M,i»-±  f  r y  Kv51^] ' 

0m  1 


where 


r2(«p)*(/?i  coscpf-a)2+(f?2  sincpffe)2, 

1— if  p2:.001, 

fHs  ifp<.ooi, 

V(<p)-f?lf?2+af?2  coscpfi/?!  sintp. 

The  series  approximation  for  Jf(p)  for  p<.001  avoids  roundoff  errors  for  p  small. 

The  accuracy  is,  of  course,  a  function  of  the  number  of  steps  N.  For  the  ratio 
of  f?i /f?2  reasonable,  i.e.,  between  ^  and  4,  a  value  of  tf*80  3eems  to  be  sufficient. 
For  extreme  cases,  i.e.,  a  ratio  of  1000  to  1,  N  must  be  taken  much  larger. 


! 


Appendix  D 


USEFUL  EXPRESSIONS 


We  will  be  using  Severn;,  relations  involving  the  gaussian  or  normal  distribu¬ 
tion  G(X)  and  its  density  function  g(X)*G' (X) ,  i.e., 


®W)-f  9(y)dy-\  SSEiadi^., 

^  —  os  v  i/v -n 


±cxp(-y2/2) 


For  completeness,  we  include  these  relations  as  lemmas,  together  with  short 
derivations, 

Lemxx  1.  Let  At  e ,  and  t  be  real t  e>0,t>0.  Then, 


Proof.  Using  (D.l),  (D.2),  and  a  change  in  variable,  we  obtain 


'■-01  - 


Ms, 

2*et' 


Changing  the  order  of  Integration,  expanding,  and  integrating,  we  obtain 


I,/  expf.-.-^^^^-^V 

\  2<«2+t2)//(27T)(a2+t2)  \/&-t2/ 


Similarly,  we  have 


r”  /  !  f(*+/l)2  s21\  dx  exp(-42/2(s2+t2))  /A//e2+t2\ 

L  «5»5r> 


We  notn  that  (D.4)  can  be  expressed  in  a  slightly  different  form: 


f  /  (*M)2\  ,  ^dx  g  (  A2  > 

exp  — 5-3 —  7==  exp-f - 


iernma  2.  Let  B,  and  e  be  real,  A^B  and  a> 0.  77!en, 


Proof.  Using  (D.l)  and  (D.2)  in  the  expression  for  Jj,  and  making  a  change 


in  variable,  we  obtain 


di/dx 
2ns  * 


Changing  the  order  of  integration,  we  have 


Is-f  dy-A-B. 
iB 


lemma  3.  Let  A,  B,  a,  b,  and  t  be  real,  A*B,  a*b,  t> 0.  Then, 


r 

G(.x)dx*AG(A)+g(A) , 

**  ~S» 


«(X+a)c(4^)+(fl+i)c(5 


-  )-  (B+a)6’(5l2.) 


Proof.  Integrating  by  parts,  we  obtain  (D.6)  immediately: 


f 


A 

G(x)dx~AG(A)ig(A) . 

—  CD 


Using  (D.6),  we  obtain 

(D.  8)  ?A  G(~^j dx*  (A+a)G +tg (4+a) . 


We  may  express  Ji,  in  the  form 


and  using  (D.8)  in  ear.h  of  the  four  integrals  in  (D.9),  we  obtain  (D.7). 
Lemma  4.  Let  A ,  L,  and  t  be  real  and  positive.  Then, 


(D. 10) 


Proof.  In  (D.7)  of  Lemma  3,  let  AmA,  Bm-A,  amL,  bm~L.  Then  (D.10)  follows, 
since  g(x)~g(-x) . 

Lernra  5.  Let  u  md  a.  be  real  and  s.  and  t  be  real  and  positive  for  n*l,  2, 

X>  v 

...,N.  Then, 


(D.li) 
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where 


i  n 

~rx 

a*  i- 


t-i  a j 


«  a. 


t*i  a? 


—  »  ai 

<2, 


a‘«a2  j  — . 
»2 
\ 


i*l  e? 


Proof.  When  we  use  the  definitions  of  8,5,  and  a2  above,  the  first  expo¬ 
nential  under  the  Integral  in  the  definition  of  I6  becomes 

exp eyp^[a?-(u-a)32^ 

Thus,  I 6  becomes 


(D. 12) 


U-  .xp^f  exp/'-tfc&aUi')  exp/ — —  Jfe-V 

\  2e2  /J-  V  2»2  /  \  2t2  227c/ 


Using  Lemma  1,  we  obtain 


(D. 13) 


*4  2+t2 


»xp(-2iai)  .xp(-i!i2»i-'). 

V  2a2  /  \  2 (e2+t2)/ 


Lemma  6.  Let  h(.xtL)’*[G(,x+L)-G(.x-l)  ]/2£.  27ien  i/’  t  and  e  are  real  and  posi¬ 


tive. 


f  J*  £V 

Ju-x\dx  ,_/u/t' ,L/t'\ 

J_„ 

?\  e  />et  At'  j 

where 


t  «/e2+t- . 


146- 


Proof. 


Making  appropriate  changes  In  variables  and  inverting  the  order  of  integration, 
we  obtain 


Using  Lemma  1,  we  have 
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